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Abstract 



In 1967, Faddeev and Popov were able to quantize the Yang-Mills theory by introducing new particles 
called ghost through the introduction of a gauge. Ever since, this quantization has become a standard 
textbook item. Some years later, Gribov discovered that the gauge fixing was not complete, gauge 
copies called Gribov copies were still present and could affect the infrared region of quantities like 
the gauge dependent gluon and ghost propagator. This feature was often in literature related to con- 
finement. Some years later, the semi-classical approach of Gribov was generalized to all orders and 
the GZ action was born. Ever since, many related articles were published. This review tends to give a 
pedagogic review of the ideas of Gribov and the subsequent construction of the GZ action, including 
many other toipics related to the Gribov region. It is shown how the GZ action can be viewed as a 
non-perturbative tool which has relations with other approaches towards confinement. Many different 
features related to the GZ action shall be discussed in detail, such as BRST breaking, the KO cri- 
terion, the propagators, etc. We shall also compare with the lattice data and other non-perturbative 
approaches, including stochastic quantization. 
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Chapter 1 

Introduction 



1.1 Patchwork quilt of QCD 

Quantum Chromodynamics (QCD) is the theory which describes the strong interaction, one of the four 
fundamental forces in our universe. This force describes the interactions between quarks and gluons, 
which are fundamental building blocks of our universe. At very high energies, QCD is asymptotically 
free, meaning that quarks and gluons behave like free particles^ However at low energies, i.e our 
daily world, due to the strong force, quarks and gluons interact and form bound states called hadrons. 
A well known example of these bound states are the proton and the neutron, but a whole zoo of 
hadrons has been observed in particle detectors. In fact, the only way to obtain information about the 
strong force is through bound states, as no free quark or gluon has even been detected. We call this 
phenomenon confinement. Although 40 years of intensive research have passed since the formulation 
of the standard model (which includes QCD), no good answer has been found to probably one of the 
most fundamental questions in QCD. Even the formulation of what confinement really is, is under 
discussion [1]. 

The difficulty for solving confinement lies in the fact that the standard techniques which have been 
so successful in QED, are not applicable in QCD. In QED the coupling constant is small enough-^, so 
one can apply perturbation theory, which amounts to writing down a series in the coupling constant. 
In QCD, the coupling constant increases with decreasing energy, a phenomenon know as "infrared 
slavery", and at low energy the coupling constant becomes too large and perturbation theory alone 
can never give a good description of the theory. Therefore, other techniques are required which we 
call non-perturbative methods. There exists a wide range of non-perturbative methods, which all try 
to approach QCD from one way or another. One should not see these different techniques as com- 
peting, but as patchwork trying to cover all aspects of QCD. Many different approaches have been 
developed to describe confinement, e.g. Abelian dominance [2, 3, 4], center-vortex dominance [5], 
light-cone dominance [6], the Kugo-Ojima confinement mechanism [7, 8], Wilson's lattice gauge the- 
ory approach [9], and the approach initiated by Gribov [10]. The last approach has been elaborated in 
widely scattered articles, but it has not been reviewed, and the relation of different approaches is not 
well known. 

Let us mention that even if one omits quarks in QCD, one can still call the remaining theory con- 

^E.g. in deep inelastic scattering (DIS) experiments, quarks can be treated as free particles. 
^This depends of course on the energy range of interest. 
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fining. Although there is no real experimental evidence, because the theory of gluons without quarks 
is a gedankentheorie, lattice simulations have shown that gluons form bound states which we call 
glueballs, and no free gluon can occur. Therefore, it is of interest to investigate pure QCD without 
quarks, and try to find out what happens. One could say that confinement is hidden in the behavior of 
the gluons. 

1.2 Gribov's gluon confinement scenario 

The gauge concept was introduced into physics by Hermann Weyl to describe electromagnetism, by 
analogy with Einstein's geometrical theory of gravity. According to Weyl, the vector potential A^(x) 
is the electromagnetic analog of the transporter or Christoffel symbol F^^ of general relativity, and 
serves to transport the phase factor of a charged field. It is remarkable that the theories of the weak 
and strong interaction which have been developed since Weyl's time and which, with electromag- 
netism, form the standard model, have also proven to be gauge theories. The group is extended from 
U(l) to SU(2) and, for QCD, SU(3), but these theories embody the geometrical gauge character intro- 
duced by Weyl. In particular they respect the principle of local gauge invariance in the same way that 
Einstein's theory of gravity respects coordinate invariance. Here we shall discuss heuristically how 
the gauge principle leads to Gribov's confinement scenario. 

According to the gauge principle, physical observables 0{A) and the action S{A) are invariant un- 
der local gauge transformations U {x) 

0{^A) = 0{A) (1.1) 
and S{^A) = S{A), where the gauge transformation is defined by 

= U{x)Aft{x)U\x) - -U{x)d^U^x). (1.2) 

(Conventions are given below.) We have written these statements for the continuum theory, but the 
corresponding statements hold in lattice gauge theory. 

Let us consider a toy example that illustrates how imposition of a gauge symmetry produces a mass 
gap. The hamiltonian for a free non-relativistic particle in one dimension is given by // = (— 1 /2m) 
With no further symmetry condition, the spectrum is continuous, with energy E = k^/2m, where k is 
any real number. Now suppose that the group of translations, x ^ x + pL, which is a symmetry of 
the hamiltonian, is taken to be a "gauge" symmetry, so the point x and the point x + pL are identified 
physically. Here p is any integer and L is a fixed length. The wave function \|/(x) is required to be 
invariant under this gauge transformation, \^{x + pL) = . Then the spectrum becomes discrete. 
En = {InnY/lmL^, where n is an integer. We describe this situation by saying that the gauge sym- 
metry has changed the physical configuration space — that is, the space of x's — from the real line 
to the circle of circumference L, and this changes the continuum spectrum into a discrete spectrum. 
This example illustrates a general phenomenon: a gauge theory must obey the constraint that the con- 
figurations A and ^A are identified physically, and such constraints have a dynamical consequences, 
making the spectrum more restrictive. According to Gribov's scenario, this changes the spectrum in 
QCD and explains why gluons are absent from the physical spectrum.^ 

^ A spectrum may of course develop a mass gap for other reasons such as the occurrence of a mass term in the action or 
by the operation of the Higgs mechanism. 
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The key step in the development of this idea was taken by Gribov [10] who showed that for a non- 
Abelian gauge theory, such as SU (2) or SU (3), the local gauge group is more powerful than for an 
Abelian theory, and imposes additional constraints. In an Abelian gauge theory the gauge condition 
3 • A = fixes the gauge (essentially) uniquely, whereas in a non-Abehan gauge theory, as Gribov 
showed, there are distinct transverse configurations, A ^A' with d - A' = d ■ A = 0, that are related by 
a 'large' gauge transformation A' = ^A. These are known as Gribov copies, and a non-Abehan gauge 
theory is subject to the additional constraint that these Gribov copies must be identified physically. 
Thus the physical configuration space is restricted to a region free of Gribov copies, sometimes called 
a fundamental modular region. Singer [1 1] showed that this situation is unavoidable in a non- Abelian 
gauge theory, and that the physical configuration space, is topologically non-trivial, in contrast to the 
situation in Abelian theories where the space of transverse configurations, satisfying 3 ■ A = 0, is a 
linear vector space. Gribov also understood the dynamical implications of the existence of Gribov 
copies, including a suppression of infrared modes due to the proximity of the Gribov horizon in in- 
frared directions. He made an approximate calculation of the gluon propagator and found that it did 
not have a physical pole. Thus the gluon is expelled from the physical spectrum as a result of the 
constraints of non- Abelian gauge invariance [10]. It was subsequently shown that the fundamental 
modular region (a region free of Gribov copies, defined below in sect. 2.2.2) in the minimal Landau 
gauge for a non- Abelian gauge theory is bounded in every direction [12]. This means that when you 
walk along any straight line emerging from the origin in the space of transverse configurations, you 
inevitably arrive at points that are (physically identified with points that are) back inside the bounded 
region, whereas in Abelian gauge theory you keep going forever. It was found in approximate calcu- 
lations by Cutkosky [13, 14, 15] and KoUer and van Baal [16], in which only a few low-energy modes 
were kept, that when gauge invariance of the wave functional is imposed, the spectrum changes from 
the perturbative spectrum of gluons to one that approximates the physical spectrum of glueballs. 

In the above remarks the physical configuration space was described in terms of a gauge fixing to 
the minimal Landau gauge. However to make sense, the physical configuration should be defined 
gauge-invariantly. We briefly indicate how this is done. The physical configuration Aphys correspond- 
ing to a configuration A is identified with the gauge orbit through A. This consists of all configurations 
A' = ^ A that are gauge-equivalent to A for all gauge transformations U = U {x),^ 

Aphys = {A':A' = ^A}. (1.3) 

The physical configuration is the space of gauge orbits, T = {J^hys}- This situation is described by the 
statement that the physical configuration space fP is the quotient space of the space of configurations 
J? = {A}, modulo the group of local gauge transformations Q = {U}, 

^ = ji/g. (1.4) 

A complete gauge fixing is a parametrization of the quotient space by choosing a single representative 
from each equivalence class. In practice the gauge must be defined in a convenient way. Singer's 

^It is worth noting tliat, in addition to gauge transformations U (x) that are continuously connected to the identity, the 
gauge orbit may also contain "large" gauge transformations that cannot be so connected. For example, if the gauge-structure 
group is the Lie group SU(2), whose group manifold is the 3-sphere S3, and if the Euclidean base-manifold (x-space) is also 
the 3-sphere 53, then a local gauge transformation U{x) is a mapping S3 — >^ S3. Such mappings, if they are continuous, 
are characterized by an integer winding number n, and the gauge orbit falls into disconnected pieces characterized by the 
winding number. 
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theorem is the statement that a choice of unique representative on each gauge orbit by a linear gauge 
condition, such as = 0, that is also continuous is impossible in a non-Abelian gauge theory. 



1.3 Approaches to the confinement problem 

There are many different approaches to the confinement problem. We shall not attempt to survey 
them here, beyond mentioning them. Lattice gauge theory, pioneered by Wilson, provides an unob- 
jectionable and numerically powerful approach, which for practical, non-perturbative calculations of 
the hadron spectrum has out-distanced the competition. However other approaches offer the possi- 
bility of analytic results, and possibly of exhibiting a simple intuitive confinement mechanism, the 
way electroweak interactions are understood by means of the Higgs mechanism. Among these is the 
approach in the maximal AbeUan gauge, with as a confinement mechanism the dual Meissner effect 
and condensation of color-magnetic monopoles, and a flux tube with a string tension. This approach 
is closely related to the approach in the center-vortex gauge. Other approaches make use of the light- 
cone gauge and the Coulomb gauge. Schwinger-Dyson calculations have been done in the minimal 
Landau gauge. Impressive results in 2-1-1 dimensions have been obtained by Nair and co-workers 
[17, 18]. 

1.4 Outline of the article 

The purpose of this review is to give a pedagogic overview of this approach by Gribov and the suc- 
ceeding construction of the so-called GZ action, which is scattered in the literature. This shall be the 
topic of chapter 2. We shall first explain the Gribov problem and discuss in detail the existence of the 
Gribov horizon with the corresponding Gribov region, which shall be of great importance for the solu- 
tion of the Gribov problem. Here, we can already show the relation of the Gribov horizon to Abelian 
and center- vortex dominance. We shall also cover many aspects of the even smaller fundamental 
modular region (FMR) and review all its properties. We shall then discuss the semi-classical solution 
of Gribov and its quantum generalization in which the cut-off at the Gribov horizon is implemented 
by a local and renormalizable action, which we shall call, as it has been called in the literature, the GZ 
action. We shall show the consequences on the gluon and ghost propagators, i.e. they are modified in 
the infrared region. 

Subsequently, in chapter 3, we shall elaborate on many different aspects of the GZ action. Firstly, 
we shall comment on the transversality of the gluon propagator. We shall also elaborate on the break- 
ing of the BRST symmetry of the GZ action, and revisit the Maggiore-Schaden construction [19], 
which attempts to interpret the BRST breaking as a kind of spontaneous symmetry breaking. We 
shall also comment on the fact that it is possible to restore BRST symmetry by the introduction of 
additional fields. Secondly, we shall investigate the right form of the horizon condition, which was 
under discussion in the literature [20, 21] and show that the solution can be found with the help of 
renormalizability. Thirdly, we shall dwell upon the famous Kujo-Ojima criterion [7], often used as an 
explanation of confinement and we shall show its relation to the GZ action. The KO criterion gives a 
relation between an enhanced ghost propagator and confinement and which is often discussed in the 
literature [22, 20, 23, 24, 25, 26]. We shall also show one should be careful with imposing a kind of 
boundary condition on the ghost propagator. Indeed one can show that imposing an enhanced ghost 
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propagator leads directly to the GZ action [27], and the KO criterion is not really fulfilled. In section 
3.7 the comparison with the lattice data shall be scrutinized in detail. We would also give a possible 
explanation for the discrepancy of the newest lattice data [28, 29, 30, 31, 32, 33, 34, 35, 36, 37] and the 
GZ action [38, 39, 40, 41, 42]. We note other analytical approaches, which are in agreement with the 
latest lattice data [43, 26, 44]. Other relevant analytic calculations of the gluon and ghost propagators 
are given in [45, 46, 109, 47, 48, 49]. In sect. 3.7.2 we present the lattice data which show clear evi- 
dence of violation of reflection positivity in the gluon propagator. It is noted that the Gribov form of 
the gluon propagator /[{k^Y + ^] also violates reflection positivity because of the unphysical poles 
at fc^ = ijy^, and it is proposed that the presence of unphysical singularities in the gluon propagator 
is the manifestation in the GZ approach of confinement of gluons. In sect. 3.8 a plane-wave source 
J^{x) = H^cos{k-x) is introduced, with free energy W{J) = Wk{H). A bound on the free energy per 
unit Euchdean volume, Wk{H) = Wk{H) /V, is reported that results from the proximity of the Gribov 
horizon in infrared directions. It implies ^mk^QWk{H) = 0, which is the statement that the system 
does not respond to a constant external color field H, no matter how strong. In this precise sense, 
the color degree of freedom is absent from the system and may be said to be confined. Finally in the 
concluding sect. 4, a critique of various approaches to the Gribov problem is presented, with a view 
toward lessons learned, open problems, and directions for future research. It includes in particular a 
review of the merits and difficulties of stochastic quantization which by-passes the Gribov problem. 

1.5 Limitations of the present approach 

In the best of worlds, an alternative perturbation series based on the GZ action, and starting with the 
Gribov propagator, would give numerically accurate results. Alas we do not live in the best of worlds 
and, as reported in sect. 3.7.1, numerical simulation on large lattices gives a finite result D{Q) > 
for the gluon propagator D{k) at k = in Euclidean dimension d = 3 and 4 [50, 51], whereas the 
alternative perturbation series gives D{0) = 0. (Numerical simulation does however give D{0) = in 
dimension d = 2.) We conclude that the GZ action successfully describes confinement of gluons by 
unphysical singularities in the gluon propagator, but the perturbation theory generated from this action 
does not give quantitatively accurate results. This deficiency of the alternative perturbation series may 
possibly be remedied by non-perturbative calculations with the GZ action [38, 39, 40, 41, 42]. 



1.6 The Yang-Mills theory: definitions and conventions 

Before starting all this, let us introduce the Yang-Mills action and establish our conventions for they 
can easily differ in different books and articles. QCD is a gauge theory, as all theories for the fun- 
damental forces. In fact, QCD is a special case of the more general SU (N) Yang-Mills theory, with 
N = 3. Therefore, we introduce the standard Yang-Mills theory and set the definitions and conventions 
used throughout this review. The derivation of the Yang-Mills action can be found in any standard 
textbook on quantum field theory [52]. 

We start with the compact group SU (N) of N xN unitary matrices U which have determinant one. 
We can write these matrices as 

U = e-'^^"^" , (1.5) 
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where X" represent the generators of the SU (N) group. For SU (2) we have Xa = Oa/2. The index 
a,b,c,...is called the color index and runs from {1, . ..,N^ — l}. These generators obey the following 
commutation rule 

[^",^*] = ifabcX', (1.6) 

and the SU{N) group corresponds to a simple Lie group. We can choose these generators to be 
hermitian, X'^ = X, and normalize them as follows 

Tr[XaXh] = ^. (1.7) 

The generators Xa belong to the adjoint representation of the group SU{N), i.e. 

UXaU'f = Xb{D^)i,a, (1.8) 

with {D^{Xa))bc = —ifabc- The structure constants of the SU{N) group have the following property, 

jrabcjrdbc ^ Nd'^ . (1.9) 

Now we can construct a Lagrangian, which is symmetric under this group. 
Firstly, we define the standard SU{N) Yang-Mills action as 

5ym = |d^Tr/vv/vv, (1.10) 

whereby 7)^, is the field strength 

Fuv = - ^y^M - ^8[AnAy] , (1-11) 
and the gluon fields which belongs to the adjoint representation of the SU{N) symmetry, i.e. 

= A^«. (1.12) 

The field strength can thus also be written as 

= F^X\ (1.13) 

whereby 

= - 3vA^ + g/a«AX ■ (1-14) 
Under the SU{N) symmetry, we define A^u to transform as 

a; = UA^U^--{d^U)U\ (1.15) 

s 

and one can check that this is compatible with the fact the A^ belongs to the adjoint representation. 
Consequently, from (1.11) we find 

F;^ = UF^U\ (1.16) 
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and therefore the Yang-Mills action is invariant under the SU{N) symmetry. Infinitesimally, the trans- 
formation (1.15) becomes 

5A° = -DfOf", (1.17) 
with the covariant derivative in the adjoint representation 

Df = a^-gr^^A^. (1.18) 
Secondly, we can also include the following matter part in the action, when considering full QCD 

Sr. = ld''x{wMapDll^) , (1.19) 

where flavor indices and possible mass terms are suppressed, and which contains the matter fields \j7,- 
and \|/,- belonging to the fundamental representation of the SU{N) group, i.e. 

% = UijWj, (1.20) 

or infinitesimally 

8^ = -igQ"X^'jWj. (1.21) 

The index / runs from {I,.. .,N}. Every \j7, and \|/, is in fact a spinor, which is indicated with the 
indices {a, P, . . .}. D'/i is the covariant derivative in the fundamental representation 

= d^d'J-igA"^{XT , (1.22) 

and j/j are the Dirac gamma matrices. One can again check that also the matter part is invariant under 
the SU (N) symmetry. The matter field \j7 represents the quarks of our model. As is known from the 
standard model, there is more than one type of quark, which we call flavors. For each flavor, we would 
need to add a term like ^m, but keeping the notation simple, we shall not introduce a flavor index here. 
The starting point of the Yang-Mills theory including quarks is thus given by 

5 = 5ym + V. (1.23) 
Most of the time, we shall however omit the matter part, and work with pure Yang-Mills theory. 
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Chapter 2 

From Gribov to the local action 



In this section we shall give an overview of the literature concerning the Gribov problem. First, we 
shall uncover the Gribov problem in detail by reviewing the Faddeev-Popov quantization [53]. Next, 
we shall treat the Gribov problem semi-classically, as done in [10]. Then, we shall translate the ideas 
of Gribov into a quantum field theory by formulating a local, renormaUzable action [54, 22]. This is 
often referred to in the literature as the GZ action and we shall use this designation. 

2.1 The Faddeev-Popov quantization 

In this section we recall the Faddeev-Popov method which solves the quantization problem of the 
Yang-Mills action at the perturbativev level [55, 56, 53]. Although this has now become a standard 
textbook item (see [57, 58, 59, 60, 61] for some examples), we shall go into the details of the calcula- 
tions to point out some subtleties. 

2.1.1 Zero modes 

We start from the Yang-Mills action given in the Introduction. Here we shall only consider the purely 
gluonic action, as the quantization problem arises in this sector. We recall that 



Unfortunately, this functional is not well defined. Indeed, taking only the quadratic part of the action. 




(2.1) 



Naively, we would assume the generating functional Z(7) (see [52]) to be defined by. 




(2.2) 




(2.3) 



and performing a Gaussian integration (A.I) 




(2.4) 
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with A/jy {x,y) = 5{x—y) (9^5^v — we see that this expression is ill-defined as the matrix A^y [x, y) 
is not invertible. This matrix has vectors with zero eigenvalues, e.g. the vector = 3pX(x), 

I dyA,y{x,y)Yy{y) = J dy[5{x -y){d\y -d,dy)]dyX{y) = 0. (2.5) 

Therefore, something is wrong with the expression of the generating function (2.2). Notice that this 
problem is present for SU (N) Yang-Mills action as well as for QED, i.e. the abelian version of the 
Yang-Mills action. 

The question is now where do these zero modes come from? Let us consider a gauge transforma- 
tion (1.15) of = 0, whereby we take U = exp(igX"x"), 

a; = --{d^U)U^ = Xad^,x\ (2.6) 

S 

or thus A^' = dfjX". This means that our examples of zero modes 7^ are in fact gauge transformations 
of Afj = 0. As we are integrating over the complete space of all possible gluon fields A^, we are also 
integrating over gauge equivalent fields. As these give rise to zero modes, we are taking too many 
configurations into account. 



2.1.2 A two dimensional example 

To fix our' thoughts, let us consider a two dimensional example. Consider an action, S{r), invariant 
under a rotation in a two-dimensional space, 

dre-*('') = / de / rdre-*^'') . (2.7) 
Jo Jo 

The "gauge orbits" of this example are concentric circles in a plane, see Figure 2. 1 . All the points 

on the same orbit, give rise to the same value of the action 5(r). Therefore to calculate W, we could 
also pick from each circle exactly one point, i.e. the representative of the "gauge orbit", and multiply 
with the number of points on the circle (see Figure 2.1). Now how exactly can we implement this? 
Mathematically, we know that for each real- valued function, we have that [62] 

with Xi the solutions of f{x) = and provided that / is a continuously differentiable function with /' 
nowhere zero. Integrating over x yields 

or thus we find the following identity 

/ dx5{fix)) = 1 . (2.10) 

^' \f'{^i)\ 

Applying this formula in our 2 dimensional plane, we can write 

— |d^5(:r(r,^)) = 1, (2.11) 



di|> 
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whereby f represents the line which intersects each orbit. Now assuming that our function J inter- 
sects each orbit only once, we can write 



3^ 



/ d^5(^(r,^)) 

J(r,(|))=0-^ 



1 



(2.12) 



However, to make an analogy with the Yang-Mills gauge theory in the next section, we rewrite this. 



ajF(r,e+^) 



3^ 



J(r,e+(|))=0 



I d^6(ir(r,e+^)) 



/ d^5(:r(r^)) = 1 , (2.13) 

r<l')=0"' 



and in this notation (|) represents the rotation angle of the vector r. This identity will allow us to pick 
on every orbit (a given r) only one representative, where ^(r^) = 0. Notice however that for every 
representative we have to multiply with a Jacobian, i.e. the derivative of the function at this point 
with respect to the symmetry parameter (|). As this Jacobian only depends on the distance r, we denote 
this measure as follows 



A^(r) 



3^ 



J(r*)=0 



Now inserting this identity into expression (2.7) gives 

W = JddJ rdrAyr{r) J d^6(5(r,e + ^))e*^W 
and transforming — > — (|), 

W = fd<\> fdQ [ rdrAj(r)5(:r(r,0))e'^W 



we are able to perform the integration over (|), which gives a factor of 27t. 

W = 2nJdQj rdrAj(r)6(5(r,0))e'^W . 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



This factor represents the "volume" of each orbit. We shall see that we obtain something similar for 
the Yang-Mills action. 



Finally, let us remark that it is of uttermost importance that jF intersects each orbit only once. Other- 
wise this derivation is not valid, and one should stick with formula (2.1 1). 
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Figure 2.1: Gauge orbits of a system with rotational symmetry in a plane and a function J which picks one 
representative from each gauge. 



2.1.3 The Yang-Mills case 

We can repeat an analogous story for the Yang-Mills action [53, 63, 64], by keeping in mind the 
pictorial view of the previous section. Due to the gauge invariance of the Yang-Mills action, we can 
also divide the configuration space A^^{x) into gauge orbits of equivalent classes. Two points of one 
equivalency class are always connected by a gauge transformation IJ = exp(— igX"6°), 

A^ = f/A^f/t--(9^[/)[/t, (2.18) 

s 

see equations (1.15) and (1.5). In analogy with the two dimensional example, we shall therefore also 

try to pick only one representative from each gauge orbit, which shall define a surface in gauge-field 
configuration space. As we are now working in a multivariable setting, i.e. infinite dimensional space 
time coordinate system x, and the —I dimensional color coordinate system a,b,.. ., the analogue 
of the identity (2.13) becomes. 

Ay I [df/]5(jr(A^)) = 1 , (2.19) 
whereby we have used a shorthand notation, 

§(iF(A^)) = nns(^'«w) 

X a 

[du] ~ nnd6"w- (2-20) 

X a 

Due to the multivariable system, the Jacobian (2.14) needs to be replaced by the absolute value of the 
determinant. 



A J (A) = I det !Kb {x,y)\ with (Mob {x, y) = — ^^^^ — 



(2.21) 



15 



This determinant is called the Faddeev-Popov determinant. Just as in the two dimensional example, 
this determinant is independent from the gauge parameter Q". 



Inserting this identity into the generating function (2.2) gives 

Z = J [dU] J [dA]Aj(A)8(jF(A^))e-^™ , (2.22) 

where we temporarily omit the source term JA. Analogous to the two dimensional example (see 
equation (2.16)), we perform a gauge transformation of the field A — > A^^ , so that transforms back 
toAju: 

A]^ = UA^U'' - - {d^U)U'' f/A^ - - {df,U)U^ = A^ . (2.23) 

§ s 

Expression (2.22) becomes 

Z = y [dU] I [dA]A^(A)8(^(A))e-^™ , (2.24) 

as the action, the measure [dA] and the Faddeev-Popov determinant are invariant under gauge trans- 
formations. Now we have isolated the integration over the gauge group U, so we find 

Z = V J [dA]A3r(A)5(5(A))e-^™ , (2.25) 

with V an infinite constant. As is common in QFT, one can always omit constant factors. It is exactly 
this infinite constant which made the path integral (2.2) ill-defined. 

Let us now work out the Faddeev-Popov determinant. This determinant is gauge invariant, and does 
not depend on 0", therefore we can choose A so that if satisfies the gauge condition !F{A) = 0. In this 
case, we can set 9" = 0, 



my) 

Applying the chain rule yields. 



(2.26) 



e=o&j(A)=o 



H.bixy) - [^,^_nM^^_^) 



(2.27) 



e=o&j(A)=o 



First working out expression (2.18) for small 9 

< = Af,-{Df,eyx"+o{e^) , (2.2s) 

and thus 



(2.29) 



J(A)=0 



This is the most general expression one can obtain without actually choosing the gauge condition, jT. 
Let us now continue to work out the Faddeev-Popov determinant for the hnear covariant gauges. For 
this, we start from the Lorentz condition, i.e. 

!F%A^{x)) = d^^%x)-B-{x), (2.30) 
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with B''{x) an arbitrary scalar field. With this condition, fW^fc(x,y) becomes, 



<Mab{x,y) = -d,Dl'h{y-x)\^^^^^^. (2.31) 

Because in the delta function in expression (2.25), the condition ir(A) = is automatically fulfilled, 
so we find, 

Z = y"[dA][det[-a^£);^*5(y-x)]]6(aA-S)e-^™ . (2.32) 

Still, we cannot calculate with this form. However, luckily, there is a way to lift this determinant into 
the action. For this, we need to introduce Grasmann variables, c and c, known as Faddeev-Popov 
ghosts. As described in the appendix in expression (A. 3), we have that (by setting t] = f] = 0) 

&QiMab{x,y) = J[dc][dc]expJdxdy^{x)!Kb{x,y)c\y), (2.33) 



and thus 

Z = J[dA] [dc] [dc]d{dA - B) exp -Sym- J dxc" {x)d^Dfc'' 
and we have been able to express the determinant by means of a local term in the action. 



(2.34) 



Finally, we would like to get rid of the dirac delta function. For this, we can perform a Uttle trick: 
since gauge-invariant quantities should not be sensitive to changes of auxiliary conditions, we average 
over the arbitrary field B" (x) by multiplying with a Gaussian factor, 

l[dB]5{dA-B)exp(J^JdxB^^ = exp(^^|dx(a,Ap2j , (2.35) 

whereby a corresponds to the width of the Gaussian distribution. Taking all the results together, we 
obtain the following gauge fixed action: 

- j dx[^d^Df"^ ^ '"-"^ 



S = 5YM+/dx(c"9^Dfc^-— (a,Ap2) , (2.36) 



This gauge is called the linear covariant gauge. Taking the limit a — )■ returns the Landau gauge. In 
this case the width a vanishes and thus the Landau gauge is equivalent to the Lorentz gauge (2.30) 
with 5 = 0. Another widely known gauge is the Feynman gauge whereby a = 1 . The Landau gauge 
has the advantage of being a fixed point under renormaUzation, while in the Feynman gauge, the form 
of the gluon propagator has the most simple form. 



In conclusion, we have obtained the following well defined generating functional: 

Z(7) = j [dA][dc][dc]exp -S + j d^^^^A^ 
with the action S given in equation (2.36). 



(2.37) 
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2.1.4 Two important remarks concerning the Faddeev-Popov derivation 



We need to make two important remarks. 

• First, notice that in fact, we need to take the absolute value of the determinant, see expression 
(2.21). In many textbooks this absolute value is omitted, without mentioning that mathemati- 
cally, it should be there. Subsequently, in equation (2.33), we have neglected this absolute value 
in order to introduce the ghosts. It was thus implicitly assumed that this determinant is always 
positive. However, in the next section, we shall prove that this is not always the case. Only 
when considering infinitesimal fluctuations around A^u = 0, i.e. in perturbation theory, is this 
determinant a positive quantity (see section 2.2). 

• Secondly, closely related to the first remark, this derivation is done under the assumption of 
having a gauge condition which intersects with each orbit only once. We call this an ideal 
condition. If this is not the case, one should in fact stay with an analogous formula such as 
(2.11), namely. 



where N{A) is the number of Gribov copies for a given orbit^ Again, in the next section, we 
shall show that the condition (2.30) is not ideal by demonstrating that the orbit can be intersected 
more than once. In fact, a mistake is made here. 

2.1.5 Other gauges 

Here we have worked out the Faddeev-Popov quantization for the Linear covariant gauges which 
encloses the Landau and Feynman gauge as a special case. However, many other gauges are possible. 
We can divide the gauges in several classes. The first class are the covariant gauges, which besides the 
linear covariant gauges also includes e.g. the 't Hooft gauge [65] and the background fields gauge [66] . 

A second class of gauges are the noncovariant gauges, i.e. gauges which break Lorentz invariance. 
The most famous example is probably the Coulomb gauge, whereby !F" = ViA'f (see e.g. [67] for a 
derivation). This gauge condition has the same form as the Landau gauge condition but ins = d—l 
dimensions. Consequently our results concerning the Gribov region, the fundamental modular re- 
gion etc. that hold for configurations A^(x) in Landau gauge, hold also in Coulomb gauge for the 
space components of configurations A,(/,x) in Coulomb gauge at a fixed time /, with the substitution 
d ^ s = d — l. A local action in Coulomb gauge that implements a cut-off at the Gribov horizon has 
been given in [68], similar to the action we shall present here shortly in Landau gauge. 

Some other examples are the axial gauge, the planar gauge, light-cone gauge and the temporal gauge. 
A nice overview on the second class can be found in [69]. Finally, there are some other gauges like 
the Maximal Abelian gauge [70], which breaks color symmetry, and some more exotic gauges which 
break translation invariance. For a nice overview on different gauges, we refer to [71]. 

For this review, we shall mainly work in the Landau gauge. 

'For each copy, the Faddeev-Popov determinant is the same, therefore, the sum in equation (2.11) can be replaced by 
\+N{A). 




(2.38) 
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2.1.6 The BRST symmetry 

Now that we have fixed the gauge, the local gauge symmetry is obviously broken. Notice however that 
the global gauge symmetry is still present as one can check by performing a global gauge transforma- 
tion on the action (2.36). Fortunately, after fixing the gauge, a new symmetry appears that involves the 
ghosts, namely the BRST symmetry. This symmetry is most conveniently expressed by introducing 
the ^7-field, 

S = 5YM + /d'^x(^Z7«a^« + a^+c%Dfc^^ , (2.39) 

whereby the path integral is now given by 

Z(y) = y^[dA][dc][dc][d^7]exp -5 + y^dx/°A^ , (2.40) 

with S here the action (2.39), and bis a. bosonic field, b" is in fact an auxiliary field, sometimes referred 
to as the Nakanishi-Lautrup field [72]. Since it has no interaction vertices, one can easily integrate out 
this field to get back the action (2.36). Therefore, the actions (2.39) and (2.36) are equivalent. It was 
found by Becchi, Rouet and Stora [73] and independently by Tyutin [74], that the action 5 of equation 
(2.39) enjoys a new symmetry, called the BRST symmetry [75], 



sS = 0, (2.41) 



with 



sc^ = b", sb'' = 0, 

= -igc"{XT^i s% = -WaC%XT ■ (2.42) 

One can check that s is nilpotent, 

= 0. (2.43) 

a property which will turn out to be very important.^ This BRST symmetry is of utmost importance, 
as it is useful for several properties. Firstly, it is the key to the proof of the renormalizability of 
the Yang-Mills action. Secondly, the BRST symmetry is also the key to the proof that the Yang-Mills 
action is unitary in perturbation theory . Let us explain what (perturbative) unitarity means. We define 
the physical state space i?^ubs, which is a subspace of the total Hilbert space, as the set of all physical 
states |^K)phys- A physical state is defined by the cohomology of the free BRST symmetry^ [76, 77] 

^o|¥)phys = and \\\f)p^y^^ so{...) , (2.44) 

where sq is the free BRST symmetry. Now a theory is unitary if 

1. Starting from physical states belonging to ^ubs, after these states have interacted, one ends up 
again with physical states G i^ubs- 



motion 

3' 



^Without the introduction of tiie b-field, the BRST operator would be nilpotent only on-shell, i.e. using the equation of 
on. 

This means, switch off interactions or set g = 0. 
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2. All physical states have a positive norm. 



It is precisely the BRST symmetry which allows one to prove these properties . Also notice that 
by fixing the gauge, we have introduced extra particles, the ghost particles c and c. Because these 
particles are scalar and anticommuting, they violate the spin statistics theorem. For the theory to 
be physically acceptable, these ghost particles must not appear in the physical spectrum. This is of 
course related to issue of unitarity and one can show that the ghosts are indeed excluded from i^ubs 
by invoking the BRST symmetry. 

Finally, let us remark that in some approaches, BRST symmetry is regarded as a first principle of 
a gauge fixed Lagrangian [8 1 , 82], rather than instead using Faddeev Popov quantization. 

2.2 The Gribov problem 

Let us now explicitly show that in the Landau gauge^, the gauge condition is not ideal. Gribov 
demonstrated this first in his famous article [10] in 1977, which has been reworked pedagogically in 
[83]. For the gauge condition (2.30), he explained that one can have three possibilities. A gauge orbit 
can intersect with the gauge condition only once (L), more than once (L') or it can have no intersection 
(L"). In [10] Gribov explains that no examples of the type L" wee known, however that many examples 
of the type L' are possible. 




Figure 2.2: The three possibilities for a gauge orbit w.r.t. a gauge condition. Original figure from [10]. 

Let us quantify this. Take two equivalent fields, and A'^^ which are connected by a gauge transfor- 
mation (1.15). If they both satisfy the same gauge condition, e.g. the Landau gauge, we call and 
A^ Gribov copies. We can work out this condition a bit further, 

A; = f/A^f/t-i(a^f/)i7t, d^^ = o & a^; = o, 

d.UA.U'f + UA.d.U'' - \{dp)U'^ - '^{d,U){d,U'^) = . (2.45) 

Taking an infinitesimal transformation, U = \ + a.,U^ = \ — OL, with a = a"X", this expression can 
be expanded to first order, 

-a^(a^a+ig[a,A^]) = 0, (2.46) 

^See [78, 79] for the original proofs, or [80] for a more recent version of the proof. 
^From now on, we shall work in the Landau gauge, unless explicitly mentioned. 
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or from equation (1.18) we see that this is equivalent to 

-df,Df,a = 0. (2.47) 
When is transverse, = 0, the Faddeev-Popov operator that appears here is hermitian, 

-a^D^ = -D^d^, (2.48) 

It has a trivial null space consisting of constant eigenvectors, a^vj/ = 0, that generate global gauge 
transformations which are unfixed by our gauge fixing. The relevant infinitesimal Gribov copies are 
in the space orthogonal to the trivial null space. 

In conclusion, the existence of (infinitesimal) Gribov copies is connected to the existence of zero 
eigenvalues of the Faddeev-Popov operator. This is a very important insight as now we can under- 
stand the two remarks made in the previous section. 

Firstly, for small A^, this equation reduces to — a^a = 0. However, it is obvious that the eigenvalue 
equation 

-a2x|/ = s\|/, (2.49) 

only has only positive eigenvalues e = > (apart from the trivial null space). This means that 
also for small values of A^ we can expect the eigenvalues e(A) to be larger than zero. However, for 
larger A^, this cannot be guaranteed anymore, so negative eigenvalues can appear (and will appear) 
for sufficiently large A, and thus the Faddeev-Popov operator shall also have zero eigenvalues. This 
means that our gauge condition is not ideal. Secondly, if the Faddeev-Popov operator has negative 
eigenvalues, the determinant of this operator can switch sign and the positivity of this determinant is 
no longer ensured. An explicit construction of a zero mode of the Faddeev Popov operator has been 
worked out in [84, 85, 83]. 

Finally, let us also mention that in QED no Gribov copies are present. We can show this with a 
simple argument. In QED, the gauge transformations are given by 

a; = A^-a^x, (2.50) 

and thus, for the Landau gauge, a^A^ = 0, the condition for A^ to be a gauge copy of A^ becomes 

a^; = => dlx = 0, (2.51) 

which does not have any solutions besides plane waves. As a plane wave does not vanish at infinity, 
they cannot be used for constructing a gauge copy A^. 

2.2.1 The Gribov region: a possible solution to the Gribov problem? 
Definition of the Gribov region 

Now that we have shown that the Faddeev-Popov quantization is incomplete, we need to improve the 
gauge fixing. Gribov was the first to propose in 1977 [10] to further restrict to a region of integration, 
the so-called Gribov region H, which is defined as follows: 

a = {a;;, a^; = 0, > 0} , (2.52) 
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whereby the Faddeev -Popov operator fM is given by equation (2.31) 

!M"'{x,y) = -d,Df8{x-y)=[-dX'+d,fai,cA';)8{x-y). (2.53) 

This is the region of gauge fields obeying the Landau gauge and for which the Faddeev-Popov operator 
is positive definite. We recall that a matrix is positive definite if for all vectors CO, 

J dxdyOi"{x)!M''''{x,y)(o''{y) > . (2.54) 

In this way the problem of the absolute value of the determinant would already be solved (first remark 
on p. 18). The border of this region 5Q is called the first Gribov horizon and at this border the first 
(non-trivial) eigenvalue of the Faddeev-Popov operator becomes zero. Crossing this horizon, this 
eigenvalue becomes negative. This is depicted in Figure 2.3. Consecutively, one can define the other 
horizons similarly, as drawn on the picture, where the second (6^2), the third (6^3), . . . eigenvalue 
becomes zero. However, keep in mind that this picture is a very simplified pictorial view. In reality, 
the space of gauge fields is much more comphcated. 




Figure 2.3: The different regions in the hyperspace dA ~ 0. 



An alternative formulation of the Gribov region 

We can also define the Gribov region as the set of relative minima of the following functional^ 

FA{U)^\\A^f = TvjdxA^^{x)A"^{x) = \jdxA^^^\x)A^^\x), (2.55) 

which corresponds to selecting on each gauge orbit the gauge configuration which minimizes A^. 
Notice that there can be more than one minimum. It can be seen relatively easily that this definition 
agrees with the Gribov region. Assume we have a gluon field for which the functional (2.55) as a 

^The original derivation can be found in [86, 87, 12], a more recent version in the appendix of [88]. 
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local minimum atU = 1. Firstly, in order to have an extremum, varying ||A|p w.r.t. an infinitesimal 
gauge transformation (1.17) must be zero, 

= 5 (^/ dxA;{x)A;{x)^ = I dx5A;ix)A;ix) = -j dxDfQ\x)A^^{x) 

= -J dxa^e"(x)A^(x) = J d;ce°(x)a^^(x) = . (2.56) 

As this equation must be zero for all 9", we must have that dijA'^{x) = 0. Secondly, this extremum 
must be a minimum, therefore, differentiating again, 

62||A||2 = - Jdxd^Q"{x)dA^^{x)= Jdxe"{x){-d^D'^'')e\x)>0 VG. (2.57) 
This implies that the operator —d^D^ = M"^ must be positive definite, see equation (2.54). 
Properties of the Gribov region 

How exactly one can implement this restriction, is the topic of the next sections. First, we shall address 
some profound geometrical questions and discuss some properties of the Gribov region. 

• Firstly, does each orbit of gauge equivalent fields intersect with the Gribov region? This prop- 
erty is of course of paramount importance as it would not make any sense to integrate over 
an incomplete region of gauge fields. A first step towards establishing this property was made 
in [10] where is was proved that for every field infinitesimally close to the horizon bQ., there 
exists an gauge copy at the other side of the horizon, infinitesimally close again. Later, it was 
then actually proven that every gauge orbit indeed intersects with the Gribov region [12, 89]. In 
[89] is was mathematically rigorously proven that for every gauge orbit, the functional (2.55) 
achieves its absolute minimum. Moreover, since every minimum belongs to the Gribov region, 
every gauge orbit intersects with the Gribov region. 

• Does = belong to the Gribov region? This is important as this means that the perturbative 
region is also incorporated in the Gribov region. In fact, we can prove this very easily. Taking 
Ap = 0, the Faddeev-Popov operator becomes = — a^S"*, which is positive definite. 

• We can also prove that the Gribov region is convex [12]. This means that for two gluon fields A^ 
and A^ belonging to the Gribov region, also the gluon field A^ = ccA^ + PA^ with a, P > and 
a + P = 1, is inside the Gribov region. To demonstrate this, we need to show that 3lf°*(ccA^ + 
PA^) is positive definite. However, from expression (2.53), we immediately see that 

5lf«^(cxAi + pA2) = aM^'^iAl^ + ^M^'^Al) . 

As a, P > 0, this sum of two positive definite matrices is again a positive definite matrix. 

• Finally we can show rather easily that the Gribov region is bounded in every direction [12]. 
This is the essential property that leads to a mass gap, as explained above in sec. 1.2. Assume 
we have a gluon field A^ / located inside the Gribov region Q.. Then we can show that, 
for sufficiently large X > 0, the gluon field XA^u is located outside of H. Firstly, as the matrix 
Mf^i^n) = ^ijfabc^^ is traceless (it is already traceless on the color indices), the sum of all the 
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eigenvalues of f^2^^ is zero. Therefore, for A^u 7^ 0, there should exist at least one eigenvector^ 
CO with negative eigenvalue k, i.e. 



J dxdy(o" (x) {x,y)(o'' (j) = K < . (2.58) 

Secondly, as ^^''^(A^) is linear in A^u, f^°^(XA^) = X!M2^{Aij) has the same eigenvector co with 
eigenvalue Xk. Therefore, 

JdxdyO)%x)!M"''{'kAf,){x,y)o)''{y) = J dx 0)%x){-dl)o)''{x) + Xk , (2.59) 

which shall become negative for large enough X. Consequently, fAf^^(AAp) is no longer positive 
definite and AA^ is located outside the horizon. Therefore, Q. is bounded in every direction. 
Moreover, in [90] it has been proven that the Gribov region is contained within a certain el- 
Upsoid. One may also find the exact location of the Gribov horizon for A which has a single 
fourier componant [91], 

A^^{x) = c''^cos{k-x), (2.60) 

where is orthogonal to for A transverse, kfjC^ = 0. Let f^v (c) be the Lorentz tensor defined 
by t/jvic) = c^Cy, and let n^{c) be its largest eigenvalue, t^Vy = n^Vft. Configurations of the 
form (2.60) intersect the Gribov horizon at 

n^{c) = 2k^. (2.61) 

Because n^{c) depends quadratically on c^, if we write c'^ = kC^, then the Gribov horizon occurs 
at n^(C) = 2. This exhibits the suppression of infrared modes by the Gribov horizon. 

With all these properties, restricting the integration of gluon fields to the Gribov region looks Uke 
a very attractive option to improve the gauge fixing. Unfortunately, the Gribov region still contains 
Gribov copies. This was first discussed in [86]. Let us repeat their reasoning. Assume a gluon field 
A^ belonging to the boundary of the Gribov region, then we have that 

6||A|p = 0, 

S2||A||2 ^ ^ 30, y"dxe''(x)fW«^(x)e^(x) = 0. (2.62) 

As the Faddeev-Popov determinant has zero modes, this means that it is inconclusive whether ||A|p 
is a minimum^. We have to consider the third variation 5^||A|p, 

5'||A|P = gfabcldxd^e%x)e\x)ir/{x)e^{x), (2.63) 

which is, generally speaking, not zero^. Therefore, a gluon field on the boundary of the Gribov region 
is not a relative minimum of the functional (2.55), and thus there must exist a transformation U: 

||A||2>||A^||2. (2.64) 

^We assume this eigenvector to have norm 1. 

^This is a consequence of the second derivative test as can be found in any textbook on basic mathematics. 

^One can compare this with whose first and second derivatives are zero atx = 0, while the third derivative is positive. 
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Since 1 1 [(1 — | and 1 1 (1 — e)A| | are continuous in £, and the difference B= | |A| p — | |A^| ^ is a 

finite number, it follows that for £ > sufficiently small, the configuration (1 —z)A satisfies the same 
inequality, — e)A|p > — e)A]^|p. Thus the configuration (1 — 8)A is not an absolute mini- 
mum and moreover because the configuration A = lies in the interior of Q. and because H is convex, 
the configuration (1 — e)A lies inside H. Thus H contains Gribov copies. 

Also numerical results have confirmed this, see e.g. [92]. In fact, it is not surprising that the Gri- 
bov region still contains copies. By looking at the functional (2.55), it seems obvious that on a gauge 
orbit, the functional (2.55) can have more than one relative minima. Two relative minima of (2.55) 
on the same orbit are Gribov copies which both belong to Q.. Also Gribov was already aware of this 
possibihty [10]. Interesting explicit examples of Gribov copies are given in [93] and references found 
there. 

2.2.2 The fundamental modular region (FMR) 
Definition of the FMR 

What is then the configuration space free from Gribov copies? It is obvious that from the functional 
(2.55) which defines the Gribov region, we can also define a more strict region, i.e. the set of absolute 
minima of the functional (2.55). As we take for each gauge orbit, the absolute minimum, we shall 
select, on a given orbit, only one gluon field, namely the gauge configuration closest to the origin. 
This region is then called the fundamental modular region A. Restricting to this region of integration 
is also called the minimal Landau gauge^^. A is then a proper subset of H, A C H. Notice that the 
absolute minimum of the functional (2.55) can only determine the minimum up to a global gauge 
transformation. Indeed, as mentioned on p. 19, fixing the gauge does not break the global gauge sym- 
metry and by performing a global gauge transformation H independent from the space time coordinate 
X, expression (2.55) does not change, 

IIA^II^ = TrJ dxHA'j^ {x)H'' HA]^ {x)H^ =Tt J dxA'^ {x)A^ (x) = \\A^\\^ . (2.65) 

Therefore, saying that we picked out from a gauge orbit exactly one configuration always means mod- 
ulo global gauge transformations. 

In fact, the FMR A would be the exact gauge fixing if the global minima of the functional (2.55) 
are non-degenerate. However, it is proven that degenerate minima can and do only occur on the 
boundary of the FMR, 6A [85]. Therefore, if one would integrate over 

Z = /[dA]e-^™, (2.66) 

J A 

these degenerate minima do not play any role, as they have zero measure. This agrees with endpoints 
of a function which do not play a role when integrating over a function. 

Properties of the FMR 

Let us discuss again some properties of the FMR, which resemble those of the Gribov region. 

• Firstly, all gauge orbits intersect with the FMR. This is in fact akeady demonstrated in the first 
bullet point on p.23. 

'"Sometimes this is also called the absolute Landau gauge, while the minimal Landau gauge can refer to taking one 
arbitrary minimum of the functional (2.55), depending on the author or article. 
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• = belongs to the FMR as is the smallest possible norm. 

• A is convex [86]. This is a bit more involved to prove than for the case of the Gribov region H. 
We have to show that if Aj,,A^ G A, alsoB^ = ?Aj, + (1 — ?)A^, with? G [0, 1]. For this we work 
out the functional (2.55), 

||A^||2 = Tr|dxA,^WA,^(x) 

= Tr dx (uAftU"^ - ^ (a^,f/)i7^^ (^UA^U"^ - ^ {d,uU)U^^ . 
= ||A|p-2-Tr / dx{A^U%U)-\Tv f dx {{d^U)U^ {d^U)U'^) . 
As Aj, and A^ both belong to the FMR, we have that 
||^i,[/||2_||^i||2>Q ^ -2-Tr j Ax{Alu''d^U)-^Tr j Ax{{d^U)U\d^,U)U^)>Q 
\\A^'"\\^-\\aY>0 ^ -2-Tr /"dx(A2f/ta^f/) -Ixr /"dx((a^[/)[/^(a^[/)f/f) >0. 
These two inequahties are linear in A, and they yield 

-2-Tr / dx ((/A^ + (1 -t)Al)U%U) - ^Tr f dx {{d^U)U\d^U)U'^) > , (2.67) 
from which follows 

\\B^\\^-\\B\f>0. (2.68) 
Therefore, B belongs to the FMR and the FMR is convex. 

• A is bounded in every direction. This is obvious as A C H with H bounded in every direction. 

• The boundary of A, §A has some points in common with the Gribov horizon [85]. 

• Some points on the boundary 8A are Gribov copies of each other. 

2.2.3 Other solutions to the Gribov problem 

Firstly, a very important result has been proven by Singer in [11], whereby it was shown that with 
suitable regularity conditions at infinity there is no gauge choice that is continuous, that is, there is no 

choice of (unique) representative of each gauge orbit that is continuous in the space of gauge orbits. 
Thus a gauge free of Gribov copies is a singular gauge that is therefore very difficult to handle in 
calculations. These gauges do exist, e.g. the space-like planar gauge [94] which has no Gribov copies, 
but which breaks Lorentz invariance. 

Many other attempts have been made, such as improving the Faddeev-Popov gauge fixing in [95] 
whereby the absolute value of the Faddeev-Popov determinant was lifted into the action. However, 
the number of copies is not properly accounted for, and therefore, as far as we know, no further cal- 
culations have been done in their framework. Also in [96, 97] an attempt to improve the gauge fixing 
has been done. However, as far as we know, the meaning of this model remains unclear in the infrared. 
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It has been pointed out [98, 99] that if one takes the Faddeev-Popov determinant, det{M{A), with- 
out the absolute value, and integrates over all configurations, including all Gribov copies, then one 
gets the correct result. This happens because one is evaluating the signed intersection number of the 
intersection of the gauge fixing surface, for example 3 • A = 0, with each gauge orbit, and this is a 
topological invariant, the same for each orbit. This has the feature that the Euclidean Landau-gauge 
or Coulomb-gauge functional weight changes sign, and the correct result is obtained by cancellation 
between different Gribov copies which give equal but opposite contribution. This could lead to large 
errors if approximations are made. 

Finally we mention that stochastic quantization [102] with stochastic gauge fixing [103] provides 
a geometrically exact solution to the quantization of gauge fields whereby a gauge-fixing "force" 
is introduced that is tangent to the gauge orbit. The stochastic process can be represented by a lo- 
cal, perturbatively renormalizable [104, 105, 106, 107] action in d -\- 1 dimensions. The idea here 
is a continuum analog of the Monte Carlo method in lattice gauge theory whereby a stochastic pro- 
cess is invented whose "time" average approaches the desired equilibrium average. Here "time" is 
machine time or the number of sweeps over the ff-dimensional EucUdean lattice. There is also a 
"time"-independent J-dimensional version of stochastic quantization. This method has perhaps been 
insufficiently explored although, in fact, the infrared critical exponent of the gluon propagator has been 
calculated in the time-independent stochastic quantization [46, 109]. The merits of this approach are 
discussed in sect. 4.5. 

2.2.4 Summary 

In conclusion, to be absolutely sure that one has a correct quantization of the Yang-Mills theory, one 
should really restrict to the FMR in order to have a completely correct gauge fixing whereby only one 
gauge configuration is chosen per orbit. However, no practical implementation of this region has been 
found so far in the continuum. Some other attempts of improving gauge fixing are interesting, but 
not very convenient or too difficult to handle. However, if we restrict ourself to the Gribov region, it 
is possible to perform practical calculations. Gribov has done this semi-classically, and as we shall 
describe in detail below, it is possible to build an action which automatically restricts to the Gribov 
region. One can still object that the Gribov region still contains Gribov copies, but there has been a 
conjecture [108, 110, 111], that the important configurations lie on the common boundary 5An 5D. of 
the Gribov region Q. and the FMR A. Therefore, the extra copies inside the Gribov region would not 
play a significant role, and it would be sufficient to restrict to the Gribov region. 

2.3 The relation of the Gribov horizon to Abelian and center-vortex 
dominance 

Let us go a bit more into the details of the hypothesis that the important configurations lie on the 
common boundary 8An 8n, as a unified discussion is not available in literature. 

We start by considering whether or not, for a given configuration A, there exists a non-zero solution CO 
to the equation, 

Z)p(A)a) = 0. (2.69) 
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This equation is gauge-invariant, as shown in the appendix A.2. So the property for a configuration 
A to allow a non-zero solution to (2.69) holds also for every gauge copy of A, Dij{^A)^(i) = 0. This 
property, which distinguishes certain gauge orbits, is thus of geometrical significance, and we may 
suspect that these gauge orbits play a distinguished role in the dynamics of QCD. 

A gauge orbit with this property has a peculiar feature. In (at least) one direction, namely the di- 
rection defined by co, the gauge orbit is degenerate because the infinitesimal gauge transformation, 

a; = A^ +£)^(A)8C0 = A^ , (2.70) 

with 8 infinitesimal, leaves A invariant. Therefore, the gauge orbit through A has (at least) one dimen- 
sion less that a generic gauge orbit. We call (2.69) the "degeneracy property". 

We now relate degenerate gauge orbits to the confinement scenario in maximal Abelian gauge. Ac- 
cording this scenario, the functional integral in the maximal abelian gauge is dominated by Abelian 
configurations (or more precisely those nearby) [4]. In SU{2) gauge theory, any Abelian configuration 
can be written as^^ 

A«(x) = 5''VW> (2-71) 

where is an arbitrary Abelian configuration. This configuration possesses the degeneracy prop- 
erty (2.69). Indeed for the the x-independent global generator, 

a" = cd"\ (2.72) 

where c is a constant, we have 3^0) = 0, and we easily verify 

D^(A)ro = 3^0) - ig[A^, ro] = , (2.73) 

which is the condition for a gauge orbit to be degenerate. Thus, the hypothesis that in the maximal 
Abelian gauge the functional integral is dominated by Abelian configurations is compatible with the 
gauge-invariant hypothesis that the functional integral is dominated by degenerate gauge orbits. 

The abelian form (2.71) is preserved by the group of local U{1) gauge transformations about the 3- 
axis, that induces the local Abelian gauge transformations, — )• a^u -t- d/jX, and we may suppose that, 
by such Abehan gauge transformations, the Abelian configurations are made transverse, 3^a^ = 0. 
Each transverse Abelian configuration corresponds to a unique distinct Abehan field tensor /^(x) = 
3^av — 3va^, with inversion a^ = (3^)~^3p/^v automatically satisfying d^a^ = 0, so different transverse 
Abelian configurations are gauge inequivalent. On the other hand the FMR, A, in minimal Landau 
gauge is bounded in every direction. So some transverse Abelian configurations lie inside A and some 
lie outside. This situation is pictured in Figure 2.4. 



In the mathematical Uterature such an orbit is called reducible. 
This can be generaUzed to other SU {N). 
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Figure 2.4: The plane of the figure represents transverse configurations 3 • A = 0. A is the fundamental modular 
region, with boundary in green, and Q. is the Gribov region with boundary in blue. The straight line 
DCOAB, in red, represents the linear vector space of transverse, Abelian configurations. A* = S'^^a^. 
The segments AB and CD that lie outside A are gauge equivalent to the segments A'B' and C'D' , 
also in red, that lie on the common boundary of A and H. 



We shall now show that when a transverse Abelian configuration A that lies outside the fundamental 
modular region A ^ A is gauge transformed A' = ^A (by a non-Abelian gauge transformation U) to 
the fundamental modular region, A' £ A, where A is the set of absolute minimum of\\^A\\^ on every 
gauge orbit, then A' lies on the common boundary of the fundamental modular region A and the Gri- 
bov region Q., A' G dKr\dQ. [110]. This is also illustrated in Fig. 2.4. 

To prove the statement we observe first that the equation, D^,(A')(o' = 0, follows from (2.69) by gauge 
invariance, where (a' = U (aU^ is now an A;-dependent gauge transformation 3^(o' / 0. From this we 
have immediately 

a^D^(A')(o' = 0, (2.74) 

where the operator on the left is recognized as the Faddeev-Popov operator. The existence of an x- 
dependent solution co' to (2.74) is the defining condition^-^ for a configuration A' £ Q. to lie on the 
Gribov horizon A' G dQ., and moreover we have A' G A C so we conclude that A' G dO.. Further- 
more since A' G dO. and A' G A and A is included in A C Q., it follows that A' necessarily also lies 
on the boundary of A, A G dA. Thus it lies on the common boundary A' G 3An QED.'"* 

This result is interesting because the Gribov horizon dD. arises as an artifact of gauge fixing in the 
minimal Landau gauge, whereas the degenerate gauge orbits have a geometrical significance. On the 
lattice there are also center vortex configurations where, for SU (2), all link variables have the value 
±1. There is a confinement scenario in the maximal center gauge (as there is for the maximal Abelian 
gauge) according to which the dominant configurations in the maximal center gauge are center vor- 
tex configurations. Center vortex and Abelian configuration both lie on degenerate gauge orbits: for 
Abelian configurations the number of missing dimensions is the rank of the gauge group, whereas 
for center vortex configurations in SU (N) there are N^-l. When the missing dimension is greater 

'•'The condition that CO be x-dependent, 3^0) / 0, is necessary because the minimal Landau gauge condition does not 
fix global gauge transformations U = expco, which have as infinitesimal generator x-independent co, with 3^0) = 0. These 
satisfy (2.74) for every transverse configuration A (including those in the interior of A), when co is x-independent, 3^co = 0, 
for we have 3^£)p(A)co = Z)^(A)3^co = 0. 

'^Note that relative minima of the minimizing functional Fj^(U) = ||^A|p for degenerate gauge orbits occur on the 
boundary 3fl by the argument given above for absolute minima. 
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than 1, then degenerate configurations are singular points of the Gribov horizon of wedge or conical 
type [111]. The proof in [111] was presented for lattice gauge theory, but the argument carries over 
to continuum gauge theory. The argument given above that, in minimal Landau gauge, degenerate 
gauge orbits intersect the common boundary, 3A fl dD., applies to both Abelian and center- vortex con- 
figurations. Thus the hypothesis that abelian configurations (center vortex configurations) dominate 
the functional integral in the maximal Abelian gauge (maximal center gauge) is compatible with the 
hypothesis that configurations on the common boundary of the fundamental modular region and of the 
Gribov region, 3An3n, dominate the functional integral in the minimal or absolute Landau gauge. 
Since the present state of QCD is a patchwork of different confinement scenarios, it is gratifying that 
the confinement scenario in center vortex gauge or maximal Abelian gauge is compatible with the 
scenario of dominance of configurations on the Gribov horizon. All three scenarios are compatible 
with a gauge-invariant scenario of dominance by degenerate gauge orbits. 

2.4 Semi classical solution of Gribov 
2.4.1 The no-pole condition 

Gribov was the first one to try to restrict the region of integration to the Gribov region [10, 83], which 
was done in a semi-classical way. He restricted the generating functional to the Gribov region by 
introducing a factor V{Q.) in expression (2.34), 



Z{J) = f [dA]exp[-5YM] 

JQ. 



J [dA] [dc] [dc]V{D.)5(dA) exp -Sym - J dxc\x)df,Dfc\. 



(2.75) 



whereby we are working in the Landau gauge, 5(3A). Now the question is how to determine this 
factor V{Q.). One can see that there is a close relationship between the ghost sector and the Faddeev- 
Popov determinant, which is clear from calculating the exact ghost propagator. For this, we start from 
expression (A.3) 

/ = j [dc] [dc] exp j d'^xd'^y Ca{x)Aab{x,y)cb{y) + j d'^x +Ca{x)Jl{x)) 



-jd'xd'yr,{x)A-J{x,y)4{y) 



(2.76) 



= CdetAexp 
where in our case: 

Aab{x,y) = -d^,Dfd{x-y). (2.77) 
From this we can calculate the ghost propagator, 

{ca{x)cb{y)), = afO05^^ 

= J [dA] V(a)8(a^p det(-a,Df )A,,^ (^, j)e-^™ . (2.78) 
Taking the Fourier transform and keeping in mind that we have conservation of momentum 

{ca{p)cbi-p))c = I [dA]\/(a)5(a^Apdet(-9^Df ) (j d{x-y)eP^^-y^Aj{x,y)] e^^^- , (2.79) 
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we can compare this expression with the one loop renormalization improved ghost propagator starting 
from the Faddeev-Popov action, 



{Ca{p)cb{k))^ 



From this expression, we can make some interesting observation. Firstly, for large momentum 
we are within the Gribov region Q., as perturbation theory should work there. Indeed, for large k^, 
g{k^) ^ \/{k^ln^), which is the perturbative result. Secondly, we notice this expression to have two 

poles: one pole at = and one pole at k^ = A^exp (-J^^f )■ Th. first pole indicates that for 

k^ w 0, we are approaching a horizon, see expression (2.80). As for all k^, "Pi is always positive, we 
stay inside the Gribov region. The second part of the ghost propagator ¥2 is not always positive for all 
k\ For k^ < A^exp (-j,^-^) , ^2 becomes complex, indicating that we have left the Gribov region. 
Therefore, V{D.) should make it impossible for a singularity to exist except at k^ = 0. 

From these observations, we can construct the no-pole condition. For this, we shall calculate g{k^,A)ah 
where the gluon field is considered as an external field. This comes down to calculating det(— 3^D^^) 
xA^^{x,y) from expression (2.79), i.e. we shall calculate the following diagrams: 



dip+k)5"''g{e) 
1 1 



487t2 



(2.80) 



IP, 



Gik,AU = 




Figure 2.5: The ghost propagator with external field to second order. 



In momentum space, these three diagrams are given by 
h = 5"'{2%y'5{k-q)^ 

In fact, this is all we can say about these diagrams, unless we take into account that after determining 
y(n), which shall be a function of the external gluon field, we shall always need to integrate over A. 
This means that the gluon lines are connected, rendering the second diagram to be equal to zero. For 

'^The Feynman rule for the ghost-ghost gluon vertex is given by ifc^/n^i, with akb resp. from c, A, and c whereby the 
outgoing momentum fe^ stems from c. 
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the third diagram, the incoming momentum k shall equal the outcoming momentum q. Formally, we 
can therefore rewrite the third diagram as 

whereby we have introduced the infinite volume factor V, to maintain the right dimensionality. More- 
over, we also know that the color indices ^ = ^. In order to calculate the correct prefactor, we take the 
sum over the color factors, see formula (1.9) 

r-^a 1 s A\ 1 , 1 1 A^g^ f d^'q ^ I ^ {k-q)f,qv 

^ {l+o{k,A)) , (2.83) 



whereby 



Now we can rewrite this 

as in this way we are considering the inverse, or the IPI diagram. This inverse contains more infor- 
mation as we are in fact resumming an infinite tower of Feynmandiagrams. The condition that the 
Faddeev-Popov operator has no zero modes, reduces to the requirement that 

G{k,A) < 1 . (2.86) 

We look into this requirement a bit more. As we are working in the Landau gauge, q^^{q) = 0, and 
thus is transverse, 

4{-q)Aiiq) = a)(A)(^5,v-^) =a)(A)P^, (2.87) 

moreover, multiplying with d/jy, we find that CO(A) = ^A^^, with d the number of dimensions. 
Therefore, we can simplify a. 

As it is possible to prove that a{k,A) decreases with increasing k^, see appendix A.3, where we have 
used the fact that A^(— ^)A^(^) is positive, condition (2.86) becomes, 

a(0,A) < 1 . (2.89) 

Taking the limit A:^ ^ in cs{k,A) yields. 
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whereby we used the fact that / ^^qf{q^)qnq\/q^ = 1/J 6^ / d^qf{q^). 
In summary, the no-pole condition is given by 

y(n) = e(i-a(o,A)), (2.91) 

with a(0,A) given by expression (2.90), or thus, using the Heaviside function, 

ViQ) = r'"'^' ^S^-oiOA)) ^ (2.92) 
we can insert this into the path integral (2.75). 

2.4.2 The gluon and the ghost propagator 
The gluon propagator 

Our goal is calculate the gluon propagator in Fourier space 

(a«(^)A^(p)) , (2.93) 

at lowest order including the restricting to the Gribov region. We start from the path integral (2.75), 
while introducing appropriate sources for the gluons, 

with N = Z^\j = 0). We do not need to take into account the integration over [dc][dc] as we are 
only calculating the free gluon propagator, also we only need the free part of Sym- Translating this in 
Fourier space, we have that 

or thus. 



whereby 



, (2.94) 

7=0 



<(^) = 8^'(P^^^§^^ + 8.vfc'+Q-l)^A) , (2.95) 

also includes the part stemming from a(0,A) in expression (2.90). Now invoking the Fourier transform 
of (A.l), we find 

A;{k)Kip)) = 5{k + p)9i I g(detO-V2(^«^^)-i(fc) . (2.96) 
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This determinant has been worked out in appendix A.4, resulting in 

^-1.^.2 ,w f A 2 , 2 1 



and thus 



(detO-i/2 ^ 



(2.97) 



As we assume {K^^^y^{k) not to be oscillating too much, we apply the method of steepest descent^^ 
to evaluate the integral over P, 



(A«(;t)A^(p)) = 8(-t + p)fA^V(W«)-i(-t)|p^p^, (2.99) 
whereby we have absorbed 27ii into fA^. Po is the minimum of /(P), i.e. 

/(Po)=0 

We define the Gribov mass, 



(2.100) 



4 ^ _Po^2 2 (2 101) 

which serves as an infrared regulating parameter in the integral. As in fact, V is equal to infinity, in 
order to have a finite y, Po ~ V. Therefore, 1/Po can be neglected and we obtain the following gap 
equation, 

which shall determine Y^. Now, we only have to calculate the inverse of 

«)(/:) = B"'(y\y^ + d,yk^+(^-l]k,k,) , (2.103) 



whereby we have set P = Po which yields, 

as one can check by calculating {K^)-\k){K^^){k) = S'^^S^. For a = 0, the inverse becomes tran- 
verse and the gluon propagator is given by 

(a;(^)A^(p)) = 6{k + pW^J^P,,ik), (2.105) 

as !AC shall cancel e^^l^"^ due to normalization. 

'^The infinite parameter to apply the method of steepest descent is the Euclidean volume V, as is illustrated by the expUcit 
calculation (2.129). 
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The ghost propagator 

Now that we have found the gluon propagator, we can calculate the ghost propagator. In fact, this 
comes down to connecting the gluon legs in expression (2.81). We easily find, as in (2.88), 



k^l-G{k) ' 



m = ^S'^ j . (2.106) 



To calculate 1 — o{k), we rewrite the gap equation (2.102) as 



k^ ^ J {2nY{q^+'f) 
and thus we write unity in a complicated way, 

To investigate the infrared behavior, we expand this integral for small k^, whereby up to order k^ 

y {k-qf 4-2^ + 1 q^ q^ W'/ ' 

and thus we can split in three parts. The first part is given by 



108) 



whereby Iy = f q2(^J_^_.^) is a number depending on y. The second part is zero, at is it odd in q, 
and the third part is given by 

R' (k) - -45 k k. [ ^"^P I Ak k. f^^^^Mn^ (2 110) 

The first term of this expression is given by, 

gap /. ^ 4 

while the second term is given by 

1 f d'^q 1 1 

Therefore, 

<W = 4(^-i5^^2 + (^25^+2^^fcv)^2^)/y (2.111) 
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Taking all results together, we obtain, 



l-G{k) 



Ng 



2kfjky 



k^ 



NgV 



d-l 



1 3 

+ ■ 



d d^ + 2d 



or thus, the ghost propagator is enhanced, 

ab /u2\ 



g^'ik 



1 d^ + 2d 



k*d^-3d + 2Ngyy 



(2.112) 



(2.113) 



As an example, for J = 4, we find easily that ly = 1/ (327t'^7^) and thus 



1 1287iV 



(2.114) 



Also in three dimension we find enhancement of the ghost. In two dimensions, the calculations are 
not so straightforward as there is a problem with switching the limit and the integration. We refer to 
[41] for more details on this. The conclusion however remains the same, also in 2 dimensions the 
ghost propagator is enhanced. 

In fact, looking at the calculations, a = 1 means that the 6-function has become a 5-function. This is 
due to the fact that we could neglect l/Po in expression (2.100) as we are working in an infinite vol- 
ume y — oo. In other words, by limiting to the Gribov region, the ghost propagator has an extra pole, 
which indicates that the region close to the boundary has an important effect on the ghost propagator. 



2.5 The local renormalizable action 

After the publication of Gribov, his result was generalized to all orders by constructing a local renor- 
malizable action [54, 22] which implements the restriction to the Gribov region, and which, following 
custom, we shall call the GZ action. In this section we shall first analyze a toy model to demonstrate 
how the GZ action was obtained [112]. 

2.5.1 A toy model 

We start with the simple quadratic action for a real scalar field in Euclidian dimension d, 

S = |d^x^(aA(x))2, (2.115) 

whereby we omit color and Lorentz indices. We assume the Gribov region is contained within the 
"ellipsoid" in A-space, 

H{A) = J d^x^A(x)[(-a2)-iA](x) = cV. (2.116) 
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It will be essential that the action S{A) and the "horizon function" H{A) are both integrals over a 
density and are thus bulk quantities, of order of the Euclidean volume V. Consequently the constant c 
remains finite in the limit V ^o°. We work at finite but large volume V, and use the fourier transform 

A(x) = y-i£Afeexp(jA:-x), (2.117) 

k 

J{x)=V-^Y,Jkexp{ik-x), (2.118) 

k 

where / J{x)A{x) = V~^Y,k-^-kAk will be a source term. This gives 

S{A) = ^^Y.^kk^^-k, (2.119) 



k 



and the Gribov region is bounded by 



' 

H{A) 

which is seen to be an ellipsoid in the infinite dimensional space of the A^. To restrict to the Gribov 
region, we need to consider the following generating functional 

Z = J[dA]Q{cV -H{A))e-^ , (2.121) 

as the 6-function assures that H{A) < cV. 



If we introduce yk=Ak/k the ellipsoid becomes a hypersphere and, as is known for hyperspheres, 
the volume gets more and more concentrated on the surface as the dimension grows. (Indeed in a 
space of dimension the volume element in radius is given by r^~^dr, and the (normalized) inte- 
gral over the interior of a sphere r < R, given by f^^B(R — r)r^^^dr, approaches Jq5(R — r)r^^^dr, 
where the 6-function is replaced by the 5-function as the number of dimensions N becomes infinite, 
oo.) Therefore, we can replace the 9-function with a 6-function and the generating functional 
(2.121) becomes 

Z = j[dA]b{cV -H{A))&-^ . (2.122) 
Let us remark that also Gribov akeady noticed this, see the end of the previous section. 



We use the formula 

^^ei(^-^>, (2.123) 

so we find 

Z = r - I [dA]e-^e(^+P)(^^-^(^)) 

J-oo 271 J 

/oo 
d/e-^(«+P) , (2.124) 
-oo 
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where G(P) = — In /[dAjeP^^'^^^^^^^e"'^. Here we have distorted the contour of integration it ^ it + 
(3 for real P > 0, which improves the convergence because H{A) > is positive. A saddle point 
approximation for the /-integration is now justified because = 0{V) is a large quantity 



-G(P*) 



whereby P* is the solution of 



G'(|3) 
cV 

cV 
cV 



/[dA]//(A)eP('^^-^W)e-^ 
/[dA]eP(^^-'f^(^))e-^ 

/[dA]g(A)e-P^W-^ 
/[dA]e-P^(^)-s 

(^(A))b ■ 



We thus obtain a generating functional of Boltzmann type, 

Z = |[dA]e-[P*^(^)+^l 



(2.125) 



(2.126) 



(2.127) 



with P* determined by (2.126). As in statistical mechanics, the micro-canonical ensemble defined by 
the "horizon condition" H{A) = cV has been replaced by the canonical ensemble [1 13] of Boltzmann 
type.^' 



We may verify by explicit calculation that with the partition function of canonical type (2.127), the 
horizon function H{A) has zero (relative) variance in the infinite volume limit, so that it is equivalent 
to the micro-canonical ensemble. More precisely we now show that {H) = 0{V), (H^) = 0{V^), and 
the variance = {H^) - {Hf = 0{V), is of order V, which is smaller than the mean-square by a 
volume factor, so the relative variance dH^/ {H'^) = 0{V~^) vanishes for V ^°o. This is the behavior 
of a generic bulk quantity in statistical mechanics. We have 



CW) = ^E^/[<lAlAtA_texp 



_Lv— — 



|[dA]exp -i-£A^(^/ + A^A(-p) + i£A;,y_ 



7=0 



1 

-r— 

2Y^ + P 
V r d'^k 1 
2 J (27ipFTp ' 



p 



2^ P 



-1 



7=0 



(2.128) 



where we have written P for P*. Here we used the Fourier transform of expression (A.l) in the 
appendix, and means to leading order in (large) V, so the sum over k may be replaced by an 

'^Here the horizon function H{A) and p* are mathematically analogous to a Hamiltonian and inverse temperature in 
statistical mechanics, and should not be confused with a mechanical Hamiltonian and physical inverse temperature. 
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integral. This shows that (H) = 0{V), as stated. We now evaluate {H'^{A)), 



4;j4 



4y2 ^ p2 dJkdJ-k^jpdJ-p 



exp 



V 



+ 



JnJ 



exp 



2V 



1^ 1 1^ 1 1^ 1 / 



d'^p 



2 J {2ny{p^ + ^)^ ' 



7=0 



(2.129) 



Thus we have {H^iA)) = {H{A)f + 0{V), and the relative variance 5H^/ {H^{A)) = 0{V-^), van- 
ishes like as asserted. We have verified that in the thermodynamic hmit the Boltzmann distribu- 
tion (2.127) does behave like a 5-function distribution. 



2.5.2 The non-local GZ action 

To restrict the region of integration of the Yang-Mills action to the Gribov region, we need to consider 
the path integral/^ 

Z = y"[dA]detfW(A)e^™e(^A)) , (2.130) 
where X{A) is the lowest eigenvalue of the Faddeev-Popov operator, 

g^-b = + = -d'd"' + gfabcA'^d, , (2.131) 

whereby we are working on-shell, d^A^ = 0. By introducing the 6-function, Q(k{A)), we insure that 
the lowest eigenvalue X{A) is always greater than zero. Note that all constant vectors 3^00 = are 
eigenvectors of the Faddeev-Popov operator, with zero eigenvalue fMoi = 0. As these eigenvalues 
never become negative, we shall not consider these trivial eigenvectors and work in the space orthog- 
onal to this trivial null-space. 



Degenerate perturbation theory 

In order to find the lowest lying (non-trivial) eigenvalue ^(A), we shall apply perturbation theory 
whereby fM^^ is the unperturbed operator. For the moment, we work in a finite periodic box of 

edge L, which shall approach infinity in the infinite volume limit. The momentum eigenstates 1*1'^'^^): 
n G Z"^/ {0}, while s runs over all the colors, s = I, . . . — I, are given in configuration space by 



x,a 



d 



18 



In this section we give a more complete derivation of the results of [54]. 
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They form a convenient basis of eigenvectors of the operator 5%, 

(0), _ /^2jty 2,vi,(0)x 



^^ra=(^^J (2-133) 

We designate by Hq the lowest non-zero momenta, i.e. no = (0) • • ■ jO, ±1,0, . . . ,0), with nQ = l. Thus 
\^H^s) represent the vectors belonging to the lowest eigenvalue^^ which is given by 

4? P.134) 

Notice that the space spanned by the vectors is 2d{N^ — l) = T dimensional, so we must apply 

degenerate perturbation theory. We call this space ^Hq 

^ = spw^i\¥ll),s=l,...,T). (2.135) 

The projector onto this space !Ho is given by 



Po = LK'M<i\- (2.136) 

nos 

The other eigenvectors of 5%, have corresponding eigenvalues X.?' given by 



2 

^£^ = (t) ^i^"'>l- (2.137) 

(The entire space can be decomposed into ^ + !Hi + fHz + . . ., whereby the Hi are defined as the 
spaces spanned by the vectors belonging to the same eigenvalue^*^). 

Let us now switch on the perturbation fMi. The T degenerate eigenvalues Xq of fWo split up into 
T different eigenvalues, Xi,...Xt, of 51/ which are its T lowest (non-trivial) eigenvalues. Within the 
degenerate subspace H, with eigenvalue X.o, any hnear combination |\|/a^^) = Y.b^abl^b''^) of the 
eigenvectors, Iv^*^^), is also an eigenvector, where Uab is an arbitrary T-dimensional unitary matrix, 
and a = {no,s). Consequently, a very small perturbation ^Mi causes a (large) finite change from the 
(arbitrarily chosen) eigenvectors Ixi/i"^), given in (2.132), into some new r-dimensional basis of eigen- 
vectors, and this change is non-perturbative. Fortunately however a small perturbation causes only a 
small change to the T-dimensional degenerate subspace !Ho itself, because the projector Pq onto !Ho, 
given in (2.136) is basis independent, and is thus adapted to any basis. Indeed, if we change basis in 
= LbUab ), we have 

= I WP) I = If/a. Ivf) (xKi"^) I = £ Wf) {^vP I = Po. (2.138) 

a ahc b 

For this reason degenerate perturbation theory is done in 2 steps. In the first step we make a trans- 
formation S, which will be calculated perturbatively, that provides a diagonalization of the Faddeev- 
Popov operator 9t{ to within a finite T xT dimensional matrix K, 

S-^^S = K. (2.139) 



'^The lowest lying eigenvalue is of course zero, belonging to the constant vectors, but we are not considering these 
constant vectors anymore. 



20 



We can order the eigenvalues by size, ^ belongs to the (« + l)th eigenvalue. 
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Here K acts within the subspace ^Hq and satisfies 



PoK = kPo = k, (2.140) 

and S satisfies 

SPq = S. (2.141) 
In the second step the diagonalization of M is completed by a simple T xT matrix diagonalization 

fZ-^Kf/ = Xdiag, (2.142) 

where X.diag is a T x T-dimensional diagonal matrix that consists of the T lowest (non-trivial) eigen- 
values, X\ to Xt of M, A.diag = diag(X.i ...X.r), and Uab is a finite T xT unitary matrix that is calculated 
non-perturbatively. For our purposes it will not be necessary to perform step 2 because it will be 
sufficient to make use of the property 

T 

Y,K = Tr(C/"'KC/) = TrK. (2.143) 



To make the first step, we seek a transformation 5, and a T x T-dimensional matrix K that satisfy 
(2.139), or 

fW5 = 5K, (2.144) 

where S and K also satisfy (2.141) and (2.140). To determine 5 and K we write them as a perturbation 
series: 

oo oo 

5=£5„ K=J^K„. (2.145) 

By substituting them into (2. 144), and identifying equal orders, we find 

9^So = 5oKo, (2.146a) 

^Si + MiSq = 5iKo+5oKi, (2.146b) 

^S2 + MiSi = 52Ko+5iKi+5oK2, (2.146c) 

9^S^ + MlS2 = 53Ko-F52Ko-F5iK2-F5oK3 , (2.146d) 



The first equation is the free equation, which is solved by 

5o = Po, Ko = ?ioPo, (2.147) 

where for simplicity we have written Xq = = (27t/L)^. To solve the higher order equations, we 
use the normalization condition, that 5^ for n > 1 maps into the space orthogonal to 

PoSn = Vn>l, (2.148) 
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and we also have 

SnPo = Sn Vn. (2.149) 
We now multiply the remaining equations by Pq, and use 

Po^Mo = ^MoPo = h)Po- (2.150) 

Firstly, we find 

PoiMiPo = Ki (2.151) 
In an analogous fashion, we find for the other equations 

K2 = PoMiSi K3 = P0M1S2 ... , (2.152) 
etc. To find 5i, we start from equation (2.146b) which may be written, upon using equation (2.151), 



5%5i-5iKo = PoO^iPo-O^iPo. 

As Kq = XqPo, we find 

(fA^, - Xo)5i = (Po - I)^iPo ^Si = i!Mo- V)"' (Po - I)MiPo 
In an analogous fashion, we can deduce S2 from equation (2.146c), 

^{^-%oI)S2 = {Po-I)^Mi{9^-XQiy\Po-I)^MiPo 
+ (fWo - V)"^ {Po -IWiPo^iPq , 

or thus 



(2.153) 



(2.154) 



52 = 



(5%-V) \Po-lWi] Po + (^^-V) ^{Po-lWiPo^iPo, 



(2.155) 



(2.156) 



whereby we made use of equations (2.151), (2.152) and (2.154). With the expressions for and ^2, 
we obtain 

Ko = 'koPo, 

Ki = Pq^MiPo, 

K2 = Po^Ui{^MG-'koiy\Po-I)fMiPo, 



K3 = Po^Mi[{9%-XoI) \Pq-I)Mi 

+PqMi (fWo - V)"^ (^0 - iWiPofMiPo . 
We can write these expressions in terms of the matrix elements. 



(2.157) 



K 



no "Oj^o ' 



' not' 



{¥^l\!Mi{!Mo-'koI) '{I-Po)iMi\% 



(0)^ 
Hot' 



4"''"' = K°ii^^i[(^^-v)"'(/-Po)fWil'i>i'g!) 



(0)\ 



(2.158) 
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The infinite volume limit 

We shall now show that in the large-volume limit, we can do some simplifications. 



In this limit we let L — > oo, while keeping a typical momentum k = 2nn/L finite. So it is advanta- 
geous to change notation, which we also simplify, and write 

\k,s)^\¥^J); {k,s\ = {¥^J\, (2.159) 

where k = {2n/L)n. Compared to a typical momentum k, the lowest non-zero momentum, ko = 
2%n()/L, with |Ho| = 1 goes to zero. Icq — 0, and the simplification comes from neglecting |^o| = 27i/L 

compared to \k\. We shall use the notation koja^ etc. for momentum vectors with lowest non-zero 
magnitude |^o| = Fol = 27t/L. 

We start with the expression for K2. We rewrite (/ — Pq) in terms of a bra-ket expansion, 

/-Po= I I |X,^)(^,^|= £ ''Y.\k,s){k,s\- (2.160) 

|i||>0 *=1 \\k\\=2%lL \\k\\>2n/L 

If we now let the operator (fWo — XqI) ^ act on this term, we obtain 

(5%-V)"'(/-A))= L ''Y^{k'-kl)-'\k,s){M. (2.161) 

\\k\\>2%lL ■'=1 

In the infinite volume limit we neglect k^ = (27t/L)^ compared to k'^ 

{k^-kl)-^ ^{k^)-^ (2.162) 
and the restriction on the summation becomes vacuous, 

^ r d'^k 

E -v/ (2.163) 

\\k\\>2Tl/L ^ ' 

where V = L'' is the EucHdean volume. Consequently we may make the substitution 

{^-'kQl)-\l-Po)^Y. L {k')~^%s){k,s\ = ^^- (2.164) 

k ■'=1 

Inserting this in k2°"'"''' yields, 

^uM, ^ ^^u\ \ks) {J^ {h\ \kot) . (2.165) 
We can work out the matrix elements with the help of equation (2.132) 

(^s\M,\kQt) = i^:o,;.^/d''xe-'(^-^»)V.rcA^W 
iko,,^ J d^^e-'^V.fcA;i(^) 

{%u\fl^i\ks) = iA(,,^^|d''xg/„,cA^(x)e'(^-^o)- 

ifco,^^|d^^g/„,,A^(x)e'^-^ (2.166) 
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where we have made use of partial integration and the Landau gauge condition d/jA/j = 0, and neglected 
^0 and which are of order 2n/L, compared to k. We thus obtain at large L, 

k 

= I d^^d^g/«.cA;(x) (^^^^j 8fstdAi{y){iko,v) , (2.167) 

which gives, for large V, the simple expression, 

^,uU ^ I g/„,,A^(x) {0^-')^ gfstdAi{y){ikoy) . (2.168) 



We now turn to K3, given in (2.158). It consists of two terms. The first term contains only pro- 
jectors I — Pq orthogonal to the subspace in intermediate positions and may be treated like K2. The 
second term in (2.158) contains the projector Pq in an intermediate position, instead of / — Pq which 
appears in the expression, 

...iI-Po)fMiPoiMi\kot) = ... \k,r){k,r\9^i\k'^,s){k'^,s\Mi\kot) 

= |^,r)i4,i|d''^e-«Vr.cA;^W 

xi^o.v:^ / d'y e-'(^'-^) -^g/,,,A^(y) , (2.169) 

where we have used (2.166). We observe that with Pq in an intermediate position (instead of / — Pq), 
there is (1) an extra factor k^ of small magnitude |^o| = In/L (instead of a finite momentum k), 
and (2) an extra factor of l/V that is associated with the finite sum L|yt"|=23i/L (instead of the sum 
1^-^ — > y / d^k/iliiY which cancels the \/V). Consequently, in the large-volume Umit, we may ne- 
glect the second term in (2. 158) compared to the first term. The first term in this equation is evaluated 
at large V by the argument used for K2, with the result for K3 given in the next equation below. 

In conclusion, in the large volume Umit, we obtain the following matrices 

^Y^t = (ik'Q^^)^jd'xd'ygUcAl{x){^-'M,M^-'),^^^^ 

etc. (2.170) 

In the large volume limit, the higher order terms K„ may be evaluated like K3. For each K„, all terms 
in which the projector Pq appears in an intermediate position are neghgible compared to the one 
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remaining term which is evaluated Uke K2. We notice that 

= %^^ -9vtQ^Mi%:^^ + Mq^9v[i^^^9\{i9v[^^^ - ... , (2.171) 
so we can sum the whole series K„ starting from n = 2, 

^o-lot ^ ^lyMt ^ 1 I d''xd^3;g/„,,A^(x) [M-') gUX{y){iko,,) • (2.172) 

Taking the trace of k 

At large volume, the spectrum of (nontrivial) eigenvalues = {littf' /Lf- > {2n/L)^ of fWo = —9^ 
becomes dense on the positive real line < < 0°. As the perturbation (A) gets turned on, the 
spectrum at first remains dense on the positive real Une. However when A crosses the Gribov horizon 
(which is bounded in every direction) the lowest eigenvalue becomes negative. Let us consider how 
this can happen. 

One possibility is that the spectrum behaves as in ordinary non-relativistic potential theory. One 
bound state develops with finite negative binding energy Xi < 0, while the rest of the spectrum X.„ for 
n> I remains dense on the positive real axis. This is illustrated in Figure 2.6, case (b). Among the T 
lowest eigenvalues Xi...Xt (which are the eigenvalues of the matrix k we have just calculated), Xi will 
be finitely negative, while X.2...Xr are at the bottom of the almost dense positive continuum starting at 
0, X„ for n = 2,...T. In this case the sum of the first T eigenvalues becomes negative when the 
first one becomes negative 

T 

Xi<0^£?i„<0. (2.173) 

Another possibility is that the spectrum remains (almost) dense, but the bottom of the spectrum moves 
a finite distance into the negative region. This is also illustrated in Figure 2.6, case (c). In this case, the 
lowest T eigenvalues all become negative (almost) together. These are the eigenvalues of the matrix 
K we have just evaluated and X„ = TrK. 
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A(0) 



HA) 





(a) 


(b) 


(c) 

A = 











Figure 2.6: Eigenvalues A,,, (A) of the Faddeev-Popov operator !M{A) at large Euclidean volume V. (a) The 
unperturbed spectrum A — Q. (b) A configuration A slightly outside the Gribov horizon, with one 
bound state, at a finite negative value, (c) A configuration A slightly outside the Gribov horizon, 
for which the almost continuous spectrum begins at a negative value. 



Next, we next evaluate TrK = TrKo + TrKi + TrK^. Firstly, the trace of Kq is given by 



TrKo= £ £ kI"'"" = XoT. 

|fl|| = l u=l 



(2.174) 



Secondly, the trace of Ki is zero. 



TrKi 



I I 

|H|| = 1 i(=l 

2n 1 



/ - / - \ JIM I i I 



/(0)x 
mi I 



\n\\ = \ u=\ 



i«o,^,Yy y gfuucAlix) = , 



(2.175) 



which vanishes because fuuc = 0, and we have used (2.166). Finally, by (2.172), the trace of at 
large volume V is given by 



TrK, = £ £ 

|h|| = 1 u=\ 



2^0 



J dWj gfurcA'^ix) (iAf-^,,., 5AvrfA^(3') • (2.176) 



Other, perhaps more pathological, cases may be considered, but we shall assume that for the configu- 
rations that dominate the Euclidean functional integral in the infinite-volume limit, it is justified, as it 
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is in the two cases just considered, to replace the condition that all eigenvalues of K be positive by the 
condition that the sum of its eigenvalues be positive, TrK > 0. In any case, this condition is weaker 
than the one it replaces, just as the restriction to the Gribov region is weaker than the restriction to the 
FMR. Consequently results obtained by this procedure will be an underestimate of the effect on QCD 
dynamics of the restriction to the FMR. 



In conclusion, the following quantity should be positive 

2 



TrK = 2(^V 



>0. 

Therefore, we can set 

Z = J[dA]e-^^''Q{d{N^-l)V -H{A)) , (2.177) 

where the integral extends over transverse configurations, S^ff = Sym — IndetfW^, where the "horizon 
function" is given by 

HiA) = jd''xjd'ygfbaiAl{x){'M-'Y^?>{x~y)ghknrAl{y). (2.178) 

It takes its name because, at large volume, the Gribov horizon is given by 

H{A) = d{N^-\)V. (2.179) 
Note that the coefficient d{N'^ — 1) is the number of components of the gluon field, Aj^. 

The non-local GZ action 

We start from (2.177), and we represent the 6-function, Q[d{N^ — 1) —H{A)\ by its fourier decompo- 



sition,-^^ 



where, for convenience, we have written 

K{A)=d{N^-\)-H{A) . (2.181) 
Upon interchanging order of integration we obtain from (2.177) 



where the "extended free energy" W(z), depending on the complex variable z, is defined by 

exp W{z) = J [dA] exp[z*:(A) - 5eff(A)]. (2.183) 

■^^As before, we have suppressed the source term {J, A) in the action. This is justified in the present section, as long as 
the source J{x) remains of compact support in the infinite-volume limit, because the relevant terms will turn out to be bulk 
or extensive quantities, of the order of the Euclidean volume V. 
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The position of tlie ie allows us to continue the path of the (O integration into the lower half plane, 
CO — >(0 — iX, with X > 0, where we will look for a saddle point. The CO integration now reads 



1 r°° 

— (ico exp[-In(X + /co) + '^1/'(X + /co)], 

'71 7-00 



271 

We look for a saddle point in X at CO = 0. If so, it is located at a stationary point of 
namely at the solution X* of 



0. 



We have 



so 



3X. 3X. X. 
exp WiX) = J [dA] exp[XK{A) - 5eff (A)] , 
dW{X) 



(2.184) 

(2.185) 

(2.186) 
(2.187) 
(2.188) 



where the subscript indicates that the expectation- value is in the ensemble (2.187), and the saddle- 
point condition reads 

^ (2.189) 



or 



{H)x = d{N^-\)V- 



1 



(2.190) 



It determines a solution X = X*. Since the Gribov region is defined by H{A) < d{N^ — the sign 
is correct for X. > 0. 

We expand <I>(?i + /co) about the saddle point X*, so the co integration now reads, to leading order 

1 



271 

X 



exp[-lnr + 'W^(r)] 
d'^^{X) 



(2.191) 



/oo 
Jco exp 
-oo 



3X2 



,coV2) 



Whether or not the saddle point method is justified and accurate depends on the sign and magnitude 
of the second derivative. 



d'^^jX) 
~3X^ 



_1_ d^-WjX) 
3X2 

1 

X2 
1 

X2 



+ {K\-{K)l 
+ {H\-{H)l>0. 



(2.192) 



It is positive, because 1/X2 is positive, as is the variance of any quantity, so the CO integration is splen- 
didly convergent. 
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We now come to an important point. The horizon function H{A) is a bulk or extensive quantity. 
Indeed, from (2.178) it may be written as the integral over a density, 



H{A) = j A^xgfMAl{x)Nbi{x). (2.193) 



where 



Nm{x) = [{M-'Y"'gfbkmA'^ {x) = I d'y {M-'Y'"b{x-y)gfbUn.Al{y). (2.194) 

Shortly we will write it as an integral over a density which is a product of local fields. It is a general 
property in statistical mechanics that the variance of a bulk quantity is of order V, 

{H\-{H)l = 0{V). (2.195) 

as is illustrated by the explicit calculation in the toy model, (2.129). This makes the coefficient, of 
(ip-, in the CO integration (2.191) also of order V, ^-^^ = 0{V). Thus the width of the peak in co is of 

order Aco = — > 0, so the saddle-point approximation becomes accurate in the limit V — > oo, 
and the co integration contributes a factor of order V^^/^. Thus, we obtain for the partition function 

Z^V-^^^expW{X*). (2.196) 

Again, it is a general property that the (extended) free energy W(k) is a bulk quantity of order V, so 
at large volume the coefficient may be neglected and from (2.187) we obtain the main result, 

Z = J [dA] exp{ X[d{N^ - \)V-H{A)] - S^s }, (2.197) 

where 7i has the value 7i* determined by (2.190). 

We have just seen that // is a bulk quantity H = 0{V). On the other hand X* is a parameter, analogous 
to the Boltzmann factor P = 1 /kT, that remains finite in the limit V — > oo, as we shall verify by explicit 
calculation. Therefore the term 1 /X,^ in the saddle-point equation (2.190) becomes negUgible at large 
V, so at large V we may write the saddle-point equation as 

^^&^=d{N^-\)-0+, (2.198) 

where 0+ is an arbitrarily small positive number. Since the Gribov region is defined by H{A)< d{N'^ — 
l)V, with H{A) = d{N^ — l)V on the boundary, this equation shows that the probability distribution 
of H{A) gets concentrated on the boundary 3n in this limit. We call this equation the "horizon 
condition." It may also be written 

In the following, we change notation slightly, by setting X = '^, and we write the last result as 

Z = J[dA] [dc] [^]e-lSYM+s^+Id''xh,{x)-fJd''x d{N^-i)] ^ (2.200) 

where hi (x) is given by 

hi{x) = j d'ygfbalA"^{x){M-'Y'"b{x-y)gfbkmAl{y) . (2.201) 
The parameter y is fixed by the gap equation 

= 'fd{N^-\)-0+. (2.202) 
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Remarks 



One thing which needs to be pointed out, is what happens with the path integral at the boundary 
of the Gribov region. At the boundary, one of the eigenvalues X of the Faddeev-Popov operator 0\{ 
approaches zero. Because hi contains the inverse of the Faddeev-Popov operator, the probability in 
the path integral vanishes rapidly; roughly speaking, a factor e^^ enters the path integral. On the other 
hand, we have argued that only the boundary of the Gribov region gives contributions for V — oo. This 
could sound contradictory. However, we can give an example which demonstrates what is going on. 
The path integral shall be a result of two competing functions. Firstly, we have a factor r^^^ stemming 
from the integration, where N approaches infinity in the thermodynamic limit, and we simple take r 
to represents the fields, while secondly, a factor e~«^ , with R the size of the boundary, represents the 
horizon function. The following function shows us what is going on in the path integral 

lim r^-'^e-T^- . (2.203) 
In the figures below, one can see how for larger N, this function evolves into a delta function 




0.5 1.0 1.5 2.0 



Figure 2.7: Evolution of the expression (2.203) for growing whereby we chose R — 3. 

Another element which looks contradictory is the following. In perturbation theory, or for large 
momenta, only the small area perturbing around A = is important, while we have just shown that the 
configurations get concentrated on the boundary of the Gribov region. Perhaps this can be explained 
by noticing that the Gribov parameter 7^ cannot be accessed in perturbation theory. Indeed, as <x 

^QCD °^ ^ ' perturbatively, y = 0. 
The correct horizon function 

In order to establish renormalizability [22], and as an infrared regulator [108], the horizon function 
(2.201) was refined into the function 

S, = limj d'x h2{x) = lim I d'x J d'^y (D^ (x)) {0^-'r{x,y) [d'; {y)y' (y)) , (2.204) 

whereby y(z) is defined through 

y2(z) = e'V- (2.205) 
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The lime_).o operation corresponds to replacing the space time dependent y^(z) with the constant Gri- 
bov parameter 'f-, and at 6 = we get D'^'^{x)y^ = fahcgA^i^, and so (2.201) is indeed recovered. We 
note here that the limit, lime_j.o, in expression (2.204) will be taken after an appropriate localization 
of the horizon function, a point which we shall outline in detail in what follows. 

In conclusion, by (2.200), the non-local action is given by 

5ni = Sym + % + 5h - / d''x d{N^ - 1) , (2.206) 

with 5ym the Yang-Mills action and 5gf the gauge fixing term, and the horizon condition reads, by 
(2.199), 

^=0, (2.207, 
or, fromZ = J[d^]e~^'^, where J[d^] stands for integration over all fields, 

{h2{x)) = yU{N^-l)-0-, (2.208) 
where h2 (x) is defined in (2.204), and we have made use of translation invariance to eliminate a factor 

y = /d^x 

2.5.3 The local GZ action 

In this section, we shall localize the action Sni by introducing some extra fields. Looking at the 
following standard formula for Gaussian integration for a pair of bosonic fields, see expression (A.2) 



Cdet-Uexp J d''xd''yf{A-^y''{x,y)j\y) 

= I [dcp][dcp]exp[| d'xd''y(^-rix)A"\x,y)^>\y))+l d'x (cp^7^(^) W^" W) 



, (2.209) 



we observe that we can get rid of the inverse of the Faddeev-Popov operator in 5/,, eq. (2.204), by 
introducing new fields. For every index /, defined hy ...f = .. we can write,^^ with A —M, and 

Jf = T,=D"^f{x)=D";f{x), 



exp {-Sh) = Y{ deiOvC j [d(p] [d(p] exp ( lim d''x J d^j cp^* {x)^M'''' {x,y)(p\ 
+ 1 d'x {D^{x)y'{x)) (pf (^) + {DK^mx)) 9- 



(2.211) 



■^^The complex conjugate notation for (p and tp is purely formal, because, if taken literally, the integration (2. 211) over ip 
and (p would explode exponentially for positive eigenvalues of 5Vf. To give a meaning to this integration, we write = i'Xp 
and we take (pj" and to be a pair of real independent bose fields, so (2.209) correctly represents the identity 

j d(pdx exp{ix^(p + ixJ + <?J) = j 5(fW(t) + 7)exp(97) = det"^fWexp(-7fW"V), (2.210) 

for every real, non-singular matrix M and real J . We note that M need not be hermitian, as happens for 9i{{A) when the 
gauge condition is taken off-shell, d-A^Q. However there is no harm in writing ^ for 1%. 
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whereby we have introduced a pair of conjugate bosonic fields {^1^-,^f^) = ((p°,(pf). We can then also 
lift the determinants det!M into the exponential by introducing a pair of Grassmann fields {(^^'^,0^^^) 
= (rof , (o"). Making use of the standard Gaussian formula for Grassmann variables 



C(detA)exp (^-| d'xd'y rM){A-'r\x,y)4{y) 

= J [dco] [dm] exp J d^xd^y m"{x)A"%\y) + J d^x {J^{x)03'' {x) + ro"(x)4(x)) 



(2.212) 



whereby we set the sources and equal to zero, we obtain 



exp(-5,)= n /[da)][da)][d(p][d(p]exp[|d^x|d^((p^(x);Af«^x,y)(pf(y) 

-coKxW"{x,y)co^iy)) +lim^ / d^'x (D^Wy^W) cp« (x) + (D« (x)y2(x)) (p;.'(x 



By (2.206), the new localized action thus becomes 

Sgz = Sq + Sj , 



with 



and 



S'o = Sym + Sg[ + I d'^x {f^dyDfi^l' - m^d^Dfa^;^ , 

- lim| d^x [(D-(x)7^(x)) (p^'X^) + (^rWY'W) +Mn' - 1)] 
limj d'^x [7^(x)D-((p«^(x)+(prW)-Y'^(A^'-l)] 

I d'x [t^d- {i?;'{x)+r;{x)) -y'd{N' - 1)] . 



(2.213) 
(2.214) 



(2.215) 



Notice that, as already remarked, the limit 9 — )• 0, in equation (2.215) has been performed after lo- 
calization. As one can see from (2.205), taking this limit is equivalent to setting y^(x) equal to the 
constant y^. At the level of the classical action, total derivatives may be neglected, Sy becomes 



5, = |d^x[/gr''^A«((p^^ + (p^^)-YMA^'-l) 



(2.216) 



Notice here that starting from the first horizon function /ii(x) given in (2.201) and undertaking the 
same procedure, we would end up with exactly the same action Sy. This can be understood as we have 
neglected the total derivatives. Although at the classical level the actions derived from hi (x) and h2{x) 
are the same, at the quantum level they may be different. 



We use the relation between the local action Sgz and the nonlocal action 5ni 

Z = J[dA] [db] [dc] [dcje"*"' = J [dA] [db] [dc] [dc] [d(p] [d(p] [dco] [droji 



-Sgz 



to obtain from the horizon condition = 0+ the local form 

- {8f''A;,{ip'; + ip';))+2y'd{N^ - 1) = 0+ , 



(2.217) 



(2.218) 
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where we have again used translation invariance to ehminate a volume factor V = J df^x. Alternately 
we can write the horizon condition as 



ar 

8f 



,2 - 0+, (2.219) 



where the "free energy," — F, is defined by 

e"^ = J [d<I)]e-^Gz , (2.220) 
and /[d^] stands for the integration over all the fields. 



There exists a freedom of redefinition of fields and we shall perform a shift of the field (o", 

<(x) ^ (Ofix) + 1 d''z(fW-i)«^(x,z)g/,,,a^[«c^)(z)(pf (z)] , (2.221) 
so that the action becomes 



Sgz = So + Sy, (2.222) 

where S'q has been replaced by So 

5o = 5[, + |d'^x(-gr^^a^ro?(D^'^c'')(pf) , (2.223) 

and an integration by parts on 3^ has been performed. For reasons having to do with BRST symmetry, 
this shift makes the action Sgz renormalizable [22], and in section 2.6, we shall present the proof 
that it is. We would like to stress that this is far from trivial, especially since no new parameter is 
needed to take into account vacuum divergences, which would lead to a modification of the vacuum 
term — / d^x d{N'^ — !)• In addition, the algebraic formalism employed in the next section also gives 
a clean and simple argument why the extra term appearing in equation (2.223) is necessary, without 
invoking the non-local shift (2.221). 

Another way to see that the term — g/"^''3^(0?(D^^c'^)(pf is allowed is to observe that it gives the 
action, considered in its dependence on the variables (c,(o) and (c,(o), the structure of a triangular 
matrix, with terms in cc, roco, and roc, but no term in cro. It is easy to see by drawing diagrams that, 
because of this triangular structure, any term in the action that involves roc — whatever its precise 
form may be — contributes only to diagrams with an entering c and an exiting ro. Thus the roc terms 
in the action do not affect correlators of pure gluons (and quarks), and we have the pleasant freedom 
to choose them at our convenience. 



2.5.4 The gluon and the ghost propagator 

Now that we have the local GZ action at our disposal, we can easily calculate the gluon and ghost 
propagator, at lowest order. We shall show that we obtain the same results as Gribov obtained, see 
section 2.4.2. 
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The gluon propagator 

To calculate the tree level gluon propagator, we only need the free part of the action Sgz, 



^GZ 



(2.224) 



where the limit a — >■ is understood in order to recover the Landau gauge. The . . . stands for the 
constant term —d{N^ — 1 )y* and other terms in the ghost- and CO, W- fields irrelevant for the calculation 
of the gluon propagator. We take y^g to be of order g^, which is equivalent to the rescaling — >■ 
Next, we integrate out the cp- and cp-fields. As we are only interested in the gluon propagator, we 

GZ ^ ^ 

w 



simply use the equations of motion, = and = 0, which give 



be 



1 



-32 



y^gfabc^a 



(2.225) 



We use this result to rewrite 5q2> 



"^GZ 



+ ... 



/ 



(2.226) 



where the last step relies on the relation (1.9). We continue rewriting so we can easily read off 
the gluon propagator 



GZ 



Ig^Nf 



8^ - 3^3v 



a 



-1 



lab 



(2.227) 



The gluon propagator can be determined by taking the inverse of A^ and converting it to momentum 
space. Doing so, we find the following expression 



?>{p + k){2Tl) 



p'^ + 2g^Nf 
" V ' 



PuPy 



:ab 



(2.228) 



which is exactly the same expression as Gribov found, see equation (3.129). We can already observe 
that this expression is suppressed in the infrared region, while displaying complex poles at p'^ = ±iy^, 
where ^ = Ig^N'^ is taken to be of order g^. This structure does not allow us to attach the usual 
particle meaning to the gluon propagator, invalidating the interpretation of gluons as excitations of 
the physical spectrum. In other words, gluons cannot be considered as part of the physical spectrum. 
In this sense, they are confined by the Gribov horizon, whose presence is encoded in the expUcit 
dependence of the propagator on the Gribov parameter y. 
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The ghost propagator 



We shall calculate the ghost propagator to one-loop order, see figure 2.8, taking into account the 
horizon condition (2.219). 




Figure 2.8: The one loop corrected ghost propagator. 



In momentum space, the ghost propagator is given by 



(2.229) 



whereby 



1 1 
1 



k^ 



N 



P,v{q) 



k^ 



'l+a{k^)) + 0{g'), 



(2.230) 



with 



o{k^ 



2 1 f d'^q ik-q)^ky 



k^J {2%y {k-qY q^ + 2g^Nf y"^" q^ 



Ng' 



2^M^v f d^q 1 



k^ J {iTlY {k-qY q'^ + lg'^Nf 



'/jv ' 



q^qy 



(2.231) 



which is identical to expression (2.106), apart from the redefinition of y. Resumming the one-particle 
reducible diagrams gives 



g{k') 



1 



1 



k^l-a{k^)' 

We are again interested in the low momentum behavior and therefore calculate cj(0). 



(2.232) 



a(0) 



Ng 



2 ^fJ^V 5? ^ 



d''q 1 



q 



Ng' 



,d-l 



d J {2%Y q^ q^ + 2g'^Nf 
d''q 1 



d J {2%Y q^ +^8'^Nf ■ 

Notice the dimensional regularization that is necessary to make this integral well defined. 



(2.233) 



To calculate it, we invoke the gap equation (2.219). Firstly, we calculate the effective action. The 
one loop effective action Fy'^ is obtained from the quadratic part of our action 



e 



"d<I)le ■^z 



(2.234) 
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with 



Notice that we need to keep the constant term —d{N^ — l)y* as it enters the horizon condition. 
A straightforward calculation gives the one loop effective action in d dimensions, 

= _.(^^_,,/+(^^^(._,)/|^,„?!±|W. (2.235, 

Now we can apply the gap equation (2.219), 

or thus 

which exactly expresses 

a(0) = 1 , (2.238) 

see expression (2.233). By (2.232) this means that the ghost propagator is enhanced, just as we 
expected from the semi-classical calculation of Gribov. Moreover, this result has been explicitly 
checked up to two loops, see [1 14, 1 15]. A more detailed calculation to this order yields 



gik')^^, (2.239) 



as originally found by Gribov. 



2.6 Algebraic renormalization of the GZ action 

We shall now prove that the GZ action is renormalizable to all orders by using algebraic renormaliza- 
tion, see [116, 117]. In [54] and [22], it was first shown that the GZ action was renormalizable. The 
Ward identities in the form given below were first given in [19]. A first algebraic proof was given in 
[83], and made complete in [21]. A recent alternative proof can be found in [1 18]. While the proof of 
renormalizability is technical, the underlying reason is simple: the term Sy, given below, that breaks 
BRST symmetry is soft, of dimension 2. 

2.6.1 The starting action and BRST 

We start with the action 

5gz = SYM + S^i + So + Sy, (2.240) 
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with 



(2.241) 



We recall that we have simplified the notation of the additional fields ((p^'',(p|^'',(jo^'',co^'') in Sq as 5o 
displays a symmetry with respect to the composite index / = (/y,c). Therefore, we have set 



sc 



((p-,(pr'«r'O=(W>9f,«r,(0«). 
The BRST variations (2.42) can be logically extended for all the fields, 

sAl = -{D,cf, " ^ 

and remains nilpotent, = 0. We note that 5o is BRST-exact, 



(2.242) 



2 

sb'' = 0, 



(2.243) 



(2.244) 



and the action of s on the field A hidden in D^^ = D"^ (A) explains the origin of the third term of 5o, 
eq. (2.241). Consequently the term 5o is an allowed gauge-fixing term within the BRST approach [22] 
that has no effect whatsover on the expectation-value of gauge-invariant quantities. However in the 
GZ action, Sqz = + Sy, the y-dependent term, Sy, breaks this BRST symmetry [54, 39], 



.5gz = ^5, = -gi I d'xf'' [a>^^- (D^c-) ((pf + <) 



(2.245) 



However the breaking is soft because the mass dimension of the integrand here, and also in Sy, is 
2, instead of 4 (in 4 dimensions). In section 3.4, we shall elaborate on the meaning of this BRST 
breaking. 

[As shown below (2.223), we have the freedom to modify the action by adding any term in roc. 
This allows us to write the starting action in the alternate form, 



where 



For we have 



Sgz = Sym + Sgf + SQ + Sy, 
Sy = -'fg J d^x ^r^^A^cp^*^ + jr^^A°ro^^ + -{N^-l)'A. 



(2.246) 
(2.247) 



-'fg J d^x (^r^"A°(p^^ + r^^(A°9j^ +Dfc'^m^;) + ^{N^-l)'f 



(2.248) 
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which indeed differs from Sy, eq. (2.241), by a term in tOc. Note that the second term in Sy is 5-exact. 
This will simplify the renormaUzation somewhat, as we shall note parenthetically along the way.] 

In order to discuss the renormalizability of Sqz [or Scz], we should treat the breaking as a com- 
posite operator to be introduced into the action by means of a suitable set of external sources. This 
procedure can be done in a BRST invariant way, by embedding 5gz into a larger action, namely 



^Gz = 5YM + 5gf + 5o + 5s, (2.249) 



where 



(2.250) 



We have introduced 3 new doublets {U^\ M"^), {V^\ Np and {T^', R'p with the following BRST 
transformations, and 

sU^'=M^\ sM^' = 0, 

sV^'=N^', sN^' = 0, 

sT^' = R^' , sR^' = . (2.25 1) 

We have therefore restored the broken BRST at the expense of introducing new sources. However, we 
do not want to alter our original theory 5gz- Therefore, at the end, we have to set the sources equal to 
the following "physical" values: 

jj<^^ I I 'j^^i I Q 

A* Iphys A* Iphys 1^ Iphys ' 

<lphys = <lphys = -<lphys = ^^.252) 

which yields 

Sgz = 2^Gz Iphys ■ (2.253) 
The doublet (T^', R^^) was introduced in [21] and provides for the first time a correct algebraic renor- 
maUzation of the action ^gz given in (2.240). This happens because the terms gf^'V^'D^^^c'^Ta'i and 
R'^^gf'^^'^D^^c'^TSi cancel at the physical value. 

[Alternatively we may embed 5gz into a larger action, 

Sgz = 5ym + % + 5o + 5s, (2.254) 

where 

5s = 5s|i?=7'=o (2.255) 

depends on only the two doublets, {U,M) and {V,N). The original action 5gz is regained at the 
physical values of the sources U, M, V, N, 



Sgz = ^GZ 

We shall provide an algebraic renormaUzation of 5gz-] 



(2.256) 

Uphysf ^physi ^phys) ^phys 
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Table 2.1: Quantum numbers of the fields. 





A" 


c" 


c" 




(Pf 




(Of 


(Of 


dimension 


1 





2 


2 


1 


1 


1 


1 


gliost number 





1 


-1 











1 


-1 


Qf -charge 














1 


-1 


1 


-1 



Table 2.2: Quantum numbers of the sources. 



jjai 


M^' 


V 


yai 
M 








L" 


2 


2 


2 


2 


2 


2 


3 


4 


-1 





1 








-1 


-1 


-2 


-1 


-1 


1 


1 


1 


1 









2.6.2 The Ward identities 

Following the procedure of algebraic renormalization outlined in [116], we should try to find all 
possible Ward identities. Before doing this, in order to be able to write the Slavnov-Taylor identity, 
we first have to couple all nonlinear BRST transformations to a new source. Looking at (2.243), we 
see that only and c" transform nonlinearly under the BRST s. Therefore, we add the following 
term to the action Egz, 

5ext = |d'^x(^-*:°(D^cr + igL°r'<^cV^ , (2.257) 
with and L" two new sources which shall be put to zero at the end, 

^^lphys = ^"lphys = 0- (2-258) 
These sources are invariant under the BRST transformation, 

sK^ = 0, sL" = 0. (2.259) 

The new action is therefore given by 

^GZ = ^GZ + 5ext- (2.260) 

The next step is now to find all the Ward identities obeyed by the action Eq^- We have listed all the 
identities below: 

1 . The Slavnov-Taylor identity is given by 

5(E^z) = 0, (2.261) 



with 



S^GZ S^GZ I ^^GZ ^^GZ 



, Lfl ^^GZ , ji^a ^^GZ , „a ^^GZ , »^a/ ^^GZ , j^rai ^^GZ , nai ^^GZ 

+^ -^^"^'W 
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2. The U{f) invariance is given by 

f/,-,E^z = 0, (2.262) 

2 5v 5v 5v 5\ 5\ 

H f-^ 

By means of the diagonal operator Qf = Uu, the /-valued fields and sources turn out to possess 
an additional quantum number. One can find all quantum numbers in Table 2.1 and Table 2.2. 



3. The Landau gauge condition reads 



4. The antighost equation yields 



-gf = d,A;. (2.263) 



^+3.^ = 0. (2.264) 



5. The linearly broken local constraints yield 



^^GZ I a ^^GZ I „f rj.di ^^GZ _ „fabcAb^7ci 

\ I fA 

S^GZ , 3 ^^GZ „rabcy^b^^GZ „rabcAbT,ci 



-^+^^W^-^f "^'^ = ^r'W- (2.265) 

6. The exact "Rij symmetry reads 

%E^z = 0, (2.266) 

with 

7. The integrated Ward identity is given by 



!Fi L'oz ^j<fx\c^~S + = • (2.268) 



8. Due to the presence of the sources T^^ and R"^, the Ghost-Ward identity is broken, an identity 
which is important in ordinary Yang-Mills theory. However, it shall turn out that this is not a 
problem for the renormalization procedure being undertaken. 

[The alternate action, 

r^z = rGz + 5ex,, (2.269) 
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satisfies tiie above 7 Ward identities with R = T = 0. In addition it satisfies tlie Gliost-Ward 
identity which is characteristic of Faddeev-Popov theory in the Landau gauge, 



-'GZ — ^cl 5 



with 



/A 



+ 8f 



■abc I —b 



rib 



ib 



"""W So? Wi ' SA^if 



and 



(2.270) 



(2.271) 



(2.272) 



where all fermionic derivatives are left derivatives. This Ghost- Ward identity, which is not 
satisfied by Sgz, simplifies the proof of renormalizability of Sgz-] 



2.6.3 The counter-term 

The next step in the algebraic renormalization is to translate all these symmetries into constraints on 
the counter-term £^2' which is an integrated polynomial in the fields and sources of dimension four 
and with ghost number zero. The classical action Eq^ changes under quantum corrections according 
to 

^GZ ~^ ^GZ + hl.hz, (2.273) 

whereby h is the perturbation parameter. Demanding that the perturbed action (Eq^ + h'Loz) fulfills 
the same set of Ward identities obeyed by see [116], it follows that the counterterm Eq^ is 
constrained by the following identities. 

[In the alternate formulation we have 

^GZ ~^ ^GZ + ^S^z- (2.274) 

The 7 identities listed below that are satisfied by the counter-term Eq2 are also satisfied by Eq^, with 
R = T = 0. As in Faddeev-Popov theory in the Landau gauge, it also satisfies the integrated Ghost- 
Ward identity, 

^"£&z = 0, (2.275) 

which is not satisfied by Eq^-I 



1 . The Unearized Slavnov-Taylor identity yields 

«E^z = 0, (2.276) 
with !B the nilpotent Unearized Slavnov-Taylor operator. 
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and 



2. The ?/(/) invariance reads 



= . (2.278) 



Uijl^'oz = 0- (2-279) 



3. The Landau gauge condition 



4. The antighost equation 



— ^ = 0. (2.280) 



'^^+3/'^ = 0- (2.281) 



5. The Unearly broken local constraints yield 



5 5 5 5 \ 



6. The exact symmetry reads 

%E^z = 0, (2.283) 

with % given in (2.267). 

7. Finally, the integrated Ward identity becomes 



/ 



d"x(."|y+co?^ + <^j=0. (2.284) 



Now we can write down the most general counterterm of d = 4, which obeys the linearized 
Slavnov-Taylor identity, has ghost number zero, and vanishing Qf number, 

5:^z = «o5ym + « I d'^x |ai^:«A« + aiB^c^A" + as L«c« + a4f/;' 3^^ + «5 V^' B^ro" 

+b^gfabcRt^\Al + b^gfabcTH'tiAl + ^5<9^a)^ + ^67;r'9^W } > (2-285) 

with ao, . . . , an and bi,...,be arbitrary parameters. Now we can unleash the constraints on the coun- 
terterm. Firstly, although the ghost Ward identity (2.270) is broken, we know that this is not so in 
the standard Yang-Mills case. Therefore, we can already set 03 = as this term is not allowed in the 
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counterterm of the standard Yang-Mills action, which is a special case of the action we are studying 



23 



[In the alternate formulation we have an identical expansion for ^q^, but with hi = . . . = b(, = 0. 
The integrated Ghost- Ward identity, Q Zq^ = 0, yields 03 = as in Faddeev-Popov theory in the 
Landau gauge. This simplifies the above argument that 03 = 0.] 

Secondly, due to the Landau gauge condition (2.280) and the antighost equation (2.281) we find, 

ai = ai. (2.286) 
Next, the linearly broken constraints (5) give the following relations 

ai = -ag = -ag = aio = an = -h = ^4 , 

a4 = as = —as = aj , bi =b2 = b$ =b^ = Q . (2.287) 

The Rij symmetry does not give any new information, while the integrated Ward identity relates the 
two previous strings of parameters: 

a\ = — ag = — ag = a\Q = an = —^3 = = 04 = 05 = — = aj . (2.288) 
Taking all this information together, we obtain the following counterterm 



c 



- gfabcV^^^c' - 8fabc^^,m'li^%c' - gfabX^d^c'^ + gfabcT^%c'^'l^ . (2.289) 

2.6.4 The renormalization factors 

As a final step, we have to show that the counterterm (2.289) can be reabsorbed by means of a mul- 
tiplicative renormalization of the fields and sources. We set ^ = {A^^.c" ,0" ^b" ,(^^ and 
= {K^,U',M^',N^\V^',U^\R'^', Tp, and we define the renormalization constants 

80 = Zf,g 
<t>o = Z^^^cj) 

0)0 = Z4,a>, (2.290) 
If we try to absorb the counterterm into the original action, we easily find. 



= 1 + 



(2.291) 



■^^In particular, since we will always assume the use of a mass independent renormalization scheme, we may compute 03 
with all external mass scales (= sources) equal to zero. Said otherwise, 03 is completely determined by the dynamics of the 
original Yang-Mills action, in which case it is known to vanish to all orders. 
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and 

7I/2 _ 7I/2 _ 7-1/4 -1/2 _ ,ai 

— Z.C ^8 ~ ^ 2 ' 

Zk = zll\ 

Zl = Zj/^ (2.292) 

The results (2.291) and (2.292) are already known from the renormalization of the original Yang-Mills 
action in the Landau gauge. Further, we also obtain 



71/2 




1/2 -1/4 _ , 

^A — ^ 


71/2 


_ 7-1/2 




71/2 


- 




Zm 


= l-/^y = 


7-1/2 -1/4 


Zn 


— ^A ' 




Zu 


1 ^U^^ 

= l+^y = 




Zv 


1 


.7-1/2-1/4 

'^g ^A ' 


Zt 


= l+/^y = 




Zr 


1 

= l-^y = 


, 7-1/27-1/4 

^g ^A 



(2.293) 



This concludes the proof of the renormalizability of the action (2.240) which is the physical limit of 
Eq2- Notice that in the physical limit (2.252), we have that 

Z^ = Zg^/^Z^^/^. (2.294) 



2.7 Relation between Gribov no-pole condition and the GZ action 

To end this chapter, we would like to point out that the equivalence between the no-pole condition and 
Zwanziger's horizon condition has been checked up to third order in the gauge fields. All the details 
of the calculation can be found in [1 19]. 
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Chapter 3 

Features of the GZ action 



In this section we shall elaborate on the GZ action by discussing a variety of topics. Firstly, we shall 
calculate all the propagators of this action and discuss the transversality of the gluon propagator. Sec- 
ondly, we devote some effort to scrutinizing the BRST breaking of the GZ action and its consequences. 
We shall also demonstrate how one can restore this broken BRST by the introduction of additional 
fields. Next, we shall briefly touch the hermiticity of the GZ action and we shall also devote some 
words on the form of the horizon function in relation to renormalizability. Finally, we shall consider 
the Kugo-Ojima criterion in relation to the GZ action. 



3.1 The propagators of the GZ action 



Firstly, we shall calculate all the propagators of the GZ action, which is useful when doing loop 
calculations with the GZ action. We start by taking only the quadratic part Srgz of the action ^gz into 
account 



5rgz 



/ 



d^x 



be 



We see three different parts appear: 

5rgz = / d4x[^(-a28^vaA)8"'A^ + ^^"3/.A^-^a.^"A;+(pfa2(p; 

The roto propagator 

The goal is to calculate the propagator 

As we are working in Euclidean space, the path integral is given by 



ah 



+ d^X 



-Sqz 



(3.1) 



(3.2) 



(3.3) 



(3.4) 
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with [d^] the integration over all the fields. In order to calculate the propagator in momentum space 
we can now employ formula (A.3) 

J d''xd'^y(Oi,ix)A^^{x,y)(xK{y) + J d^'x {JUx)(0^ix)+m^{x)J^{x)) 

j d'^xd'^yj!^{x)A-^\x,y)Jl{y) 



[dco] [dro] exp 
= CdetAexp 
where in our case: 



(3.5) 



A{x,y) = b{x-y){d''W , 
A-\x,y) = b{x-y)^^^' . 



(3.6) 



Now going to Fourierspace, 



1 = j [da)][da)]exp j ±^(o{-p){-p^Mp) + J ^ (U-pMp) + oi{p)M-p)) 



f d"p ~ I ~ 
= CdetAexp -J j^M-p)^Mp) 

we can calculate the propagator in Fourierspace: 
6 6 
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_ab , ~~cdnA\ — Siac Sibd si^N ^ Sir „ i ;,^^'^^^4• 



{p)K{k)) = 8^^8''^8^^^8(p + /:)(27i 

7=0 • ' P 



(3.7) 



(3.8) 



The ghost propagator 

Completely analogously, we find 



c''{k)d'{p)) = B^>^-{2n)'5{p + k). 



(3.9) 



The mixed operators 

For the final part, we shall use the following general formula. 



I{A,J) = J [d(p]exp d'^xd'^y(p{x)Aix,yMy) + J d'^x (p(x)/(^ 



= C(detA)-i/2exp 



yd'xd'yJ{x)A-\x,y)J{y) 



(3.10) 



However, as we see in the action (3.2) we need to rewrite the complex conjugate Bose fields (p and (p 
in terms of real fields. Therefore, we shall introduce the real fields U and V 



2 

— i 

T 



yafc _ iJjab 



(3.11) 
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We can thus rewrite the relevant part of the action as 



We rewrite this in matrixform as 

exp[-5;,Gz] = ^M-ll d'x [A",{x) b'^jx) V^'{x) Uj^^x)] 



(3.12) 













■.an 







13^ 











gpxg^pgf^ 



b"{x) 
Vi^ix) 



where 



P -g .^A-^v 



32 



3^3v 

32 



(3.13) 



are the respective transverse and longitudinal projectors. Notice that we have rewritten (3.12) in a 
symmetric way. Now we can apply the general formula (3.10), meaning that we have to find the 
inverse of A. 

■ -32p^5'^^ 

3v8'"^ 

-fgf^^^a:^ 





3^8^^" 


-fgr%yi 


" 




" Ab c 
T 




f~b xy 

*-'V CO 















pn c 


pn 


<-f CO 


Tjn gh 
^ X 





gaKgfa'g£y 







jij c 

■■K T 


jij 

•'K 


CO 


T ij gh 
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gPXg^pgt9_ 
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^Xco 
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^^x -1 
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g^ggrtg 
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After some calculation we find for A ^ (x) : 



9^ 



p fiic -2gf ;s fiio -2gf fijrfxyrp ^g^f 





For the propagators, we need to go to Fourierspace, 



-92 










/[dtp] 



exp 



\jd'xX{x)A{x)X^{x) + J ^,n-p)J{p) 

C{daA)-'l'^^v\f ^j'{-p)A-\p)J{p), (3.14) 



with 



and A ^ (p) given by 



-2gf 



-iPv 



p fi/c -^ST a fijozlsf fijrfxyrp ^g'T' 

^kt/ p4+2g^Nf p* J J p2{p-*+2g2Nf) 



[Ja h Jv Ju] , 

fbxyp -'^rf 
4gY 



(3.15) 



















llpgyi^i. 



We now have all the ingredients to calculate the propagators. 

AA-propagator 

We have for example, 



1 



hJA^yp)hJ^,{-k) (271) 

1 



8 \ 



or equivalently 



where we have defined 



_p4 + 2A^gY 



.2.4 ^A'V 



(3.16) 



(3.17) 
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A^7-propagator 

Next, 



or thus 



= T^siWnk)) = -i^6«^%^, (3.18) 



4{p)b''{k)) = -i^5"''d{p + k){2n)\ (3.19) 



^7^7-propagator 



b"{p)b\k)) = d-''^d{p + k){2n)\ (3.20) 



The propagators with U and V 

In an analogue fashion, we find 

A",{pK'^{k)) = r''^^^,P^{p){2n)'d{p + k) 

-2gf 
P^{P^ 



b^pK^-ik)) = fb^ip,-^{2nf5{p + k), 

vf{p)vf{k)) = {r''f''P.^ p2j'X'^^) ^y^^^^'%^ {2nrHp+k) , 



Uf{p)Uf{k)) = -^5-5''^5^(27i)45(;7 + ^:), 



;A;(;,)t/f (^)) = (b^p^'^ik)) = (vf{p)Ufik)) = 0, (3.21) 
which can be rewritten in terms of (p and (p again, 

A^(p)9v^('t)) = (A^(p)^;'(-t)) =r^^^=|^P^(p)(27c)^8(p + fc), 
b%p)<?','{k)) = lb%p)t\k)) = r'^ip,^{2i()%p + k), 



(3.22) 
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In summary, we have the following large set of propagators in the theory: 



Al,{p)b\k)) = -i^5"'5{p + k){2n)' 



b'^{p)b\k)) = 5^'^5{p + k){2n)\ 



p' 

b\p)^l'{k)) = (b'^ipnl'ik)) = r'^ip,^i2nr5ip + k), 
^»%{k)) = (r^T-^-P^ ^,(^4+Tg2^^) + ^5-5^^5,v) {2nrd{p + k) , 



(3.23) 

3.2 The transversality of the gluon propagator 

One might wonder whether the gluon propagator still remains transverse in the presence of the Gribov 
horizon, this can however be proven easily as done in [39]. We start with the following Ward identity, 
as given in expression (2.263) 

whereby F is the quantum effective action. Introducing sources I''{J^) for the fields fe°(A°) and per- 
forming the Legendre transformation, this identity translates into 

'• = S,§. (3.25) 

whereby is the generator of the connected Green functions. Subsequently, acting with on this 
expression, and by setting all sources equal to zero, we retrieve 

52Z^ 



= 3; 



^ 5y«(x)5y^^(y) 



= 9^(A^(x)A^(j)). (3.26) 
/,/=o 



As the gluon propagator is in fact the connected two-point function, we have proven the transversality 
of the gluon propagator. 
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3.3 The hermiticity of the GZ action 



Let us also shortly comment on the hermiticity of the GZ action. We recall that the action is given by 
equation (2.240) 



with 



5gz = So + Sy, (3.27) 

So = Sym + S,f + I d'x (cp^3, (d^ (pf ) - ^d, (D;f cof ) - gr'%m^D';'c'i^) , 

Sy = -Y2gy'd'^x(^r^''A«(p^'^ + r^'^A«9f + ^(A^2_i^^^ (328) 



If we define 



(p^ = 9 , 9^ = 9 , co^ = CO , co^ = CO , (3.29) 



we see that the GZ action is almost Hermitian, provided that the gauge condition 3 • A = is treated 
on-shell,^ up to the term — g/''*'^3^(o"D^^c'^(p?. However, we recall that this term was introduced for 
renormalization reasons by the shift (2.221). Therefore, returning to the action before the shift, the 
GZ action is Hermitian. Moreover, for practical purposes, the term — g/"^''3^a)?Dj;'^c^(p- is almost 
redundant as it cannot couple to any Feynman diagram without external c and CO legs, as discussed in 
the paragraph below (2.223). Therefore, we can conclude that in practice, the GZ action is Hermitian. 
[Similar remarks hold also for 5gz-] One comment which one can make is that the kinetic term (p3^(p 
in the GZ action, seems to have to wrong sign, as does the term ^^d/j {p^^(^^), when the Faddeev- 
Popov operator is positive. This has already been noted above, and explained in (2.210). 



3.4 The soft breaking of the BRST symmetry 
3.4.1 The breaking 

We recall here that the GZ action (2.240) is not invariant under the BRST transformation (2.243). 
Indeed, in equation (2.245) it was shown that 



Ay — SS — SSy 

or alternatively 



(3.30) 



A, = sS = sSy =-gfl d'xr'' (Ay; - (D^c-) ip';) . (3.31) 

One sees that it is exactly the presence of the Gribov parameter y which prevents the action from being 
invariant under the BRST symmetry. Indeed, when y = 0, we can integrate out the fields (p,(p,(0, CO 
and we are left with the original Yang-Mills theory. Therefore, this breaking is clearly due to the 
introduction of the horizon into the Yang-Mills action. Nevertheless, this fact does not prevent the 



'The gauge condition may be treated on-shell provided the Z)-field is shifted appropriately under hermitian conjugation. 
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use of the Slavnov-Taylor identity to prove the renormalizability of the theory. Since the breaking Ay 
is soft, i.e. it is of dimension two in the fields, it can be neglected in the deep ultraviolet, where we 
recover the usual notion of exact BRST invariance as well as of BRST cohomology for defining the 
physical subspace, see section 2.1.6. However, in the nonperturbative infrared region, the breaking 
term cannot be neglected and the BRST invariance is lost. Breaking of BRST invariance may be 
regarded as a consequence of the fact that the non-perturbative gauge fixing introduced above has 
been done in the Landau gauge, and has not been carried out in other covariant gauges. In fact it has 
been shown that if the GZ action is modified by Sqz Sgz + ^b^, then the expectation-value of gauge 
invariant quantities becomes dependent on the would-be gauge parameter a [120, 121]. According 
to the Maggiore-Schaden construction [19] that will be described shortly, the soft BRST breaking 
introduced here may alternatively be regarded as a spontaneous symmetty breaking. 



3.4.2 The BRST breaking as a tool to prove that the Gribov parameter is a physical 
parameter 

The breaking term (3.31) has an interesting consequence as it allows one to give a simple algebraic 
proof of the fact that the Gribov parameter y is a physical parameter of the theory [39], and that 
as such it can enter the explicit expression of gauge invariant objects: e.g. the correlation functions 

{F^{x)F^{y)) or the vacuum condensate (^^). 

We can demonstrate this as follows, see [39]. Taking the derivative of both sides of equation (3.31) 
with respect to one gets. 



= ^ A, = -§1 dV"'^' (a>^'- - (Z^rc-) {^'; + ) . (3.32) 



35 



As we have seen, the BRST operator s which is defined in equation (2.243) is nilpotent, thus we have 
that 1^ cannot be cast in the form of a BRST exact variation, namely 



dS 
3f 



a/^Ay, (3.33) 



for some local integrated dimension two quantity Ay. Otherwise would have to be equal to zero. 
To show indeed that is a physical parameter, we assume for a moment the contrary 

sSy = 0, (3.34) 

whereby there is no breaking term Ay. Since Sy depends on the auxiliary fields ((p^*^, 00°*^, 00°'^) 
which constitute a set of BRST doublets, we know that, 

5y = j5y. (3.35) 

due to the doublet theorem, see appendix B.2,i.e. Sj is 5-exact. Subsequently, taking the derivative of 
both sides of expression (3.35) with respect to one obtains 

4' '4- 
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which implies that 7^ is an unphysical parameter. Indeed, for a gauge invariant quantity Q, we would 
have that 

^ = I m g e-^°^ ~ j m g sSyC-'°- ~ J ms Sy) e-^- ~ . .)) = (3.37) 

and correlation functions of gauge invariant operators would be completely independent from 7^. 

However, when sSy 7^ 0, due to the presence of the soft breaking term Ay, we have that ^ / 0. 
Therefore, the existence of the breaking Ay could be an ingredient to introduce a nonperturbative mass 
gap in a local and renormalizable way. 



3.4.3 The Maggiore-Schaden construction and spontaneous breaking of BRST sym- 
metry 

Let us also discuss the following paper [19]. The authors of this paper proposed to interpret the BRST 
breaking as a kind of spontaneous symmetry breaking. For the benefit of the reader, we shall firstly 
explain the Maggiore-Schaden construction, and then clarify its status. 



The Maggiore-Schaden construction explained 

In [19], the first step was to add the following BRST exact term to the Yang-Mills action 5ym: 

5i = jyd^x(c«a^^ + ro^^av£>f(pj;') , (3.38) 

with s, the usual nilpotent BRST operator as defined in (2.243). The first term is the usual gauge fixing 
5gf, while, the second term introduces the new fields W^, (H^, and (p^^*^. As s is nilpotent, the action 
Sym + does not break the BRST symmetry. Moreover, we recall that, see expression (2.243) 

sm^ = (pf, 59f = 0. (3.39) 

Therefore, as the nilpotent symmetry defines two doublets (9, co) and (9, ro), and by appendix B.2, we 
can exclude these fields from the physical subspace. The second step was to shift a number a fields in 
the following way 

= b'^^+gi'f'^'x,, (3.40) 

where the primed fields {v^'^;^ , cp^"^ , c'", b"") are new fields which all have vanishing vacuum expectation 
value and which respect translation invariance. The BRST operator s acts on these new variables 
according to 

sc"' = b'\ sb'" = 0, 

s(p'f = (of , s&f = . (3.41) 
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Notice now that we can write s(}f^ as follows, 



(3.42) 



and thus, this BRST variation becomes x-dependent. Moreover, also the vacuum expectation value 
becomes x-dependent, i.e. 

(3.43) 



SCO 



ah 



Because the expectation value of an s'-exact quantity is broken, and because the action is s-exact, it 
is concluded in [19] that the BRST symmetry s is spontaneously broken. As a third step, the action 
5ym + 51 is rewritten in the primed variables. One finds for 5i 



5i = s|d4x(c'«a^« + ra-avDf<-gY'</ai>cA^) 
After working out the BRST variation s. 



(3.44) 



si=j d^x [z7'«a^^+c'''a^ (^f c^)] + / 



+ J d\[-gij^'^'fabcAl-gi^ffabcDl'c' 



(3.45) 

one finds that the action Sym + is very similar to the GZ action, up to three parts. Firstly, the 
constant part A'f{N'^ — 1) is missing, therefore a BRST invariant part ^2 = —'f's J d^xdfjiS^" is added 
to the action. This gives. 



(3.46) 



where we have neglected the first term as it is a total derivative of a field that respects translation 
invariance. Secondly, the last term on the first line of expression (3.45), i.e. y^X/jd^D'^^^)'^'' needs to 
be deformed. If we naively assume that we may perform a partial integration and neglect the surface 
terms we find 

J dWxudyD^Wt = - j dW£>f = -i^gf^' j d^xA«(p;f^ (3.47) 

Thirdly, comparing with expression (2.240), we see that there still is an extra term present which 
precisely yields the alternate action 5gz, eq. (2.248): 



= Sgz ■ 



fd^d 



(3.48) 



In conclusion, we can say that the action 5gz (which differs from 5gz by a physically irrelevant term) 
is obtained from an exact variation 
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3.4.4 Local version of Maggiore-Schaden construction 



In [39] this construction was re-examined and some comments were made. The main point of criticism 
is assumptions that were made concerning the partial integration with neglect of surface terms at 
infinity. The main point of criticism are the assumptions which have been made concerning the partial 
integration. Looking at the partial integration (3.47), one sees that the surface terms are neglected. 
Normally, this is not a problem as we always deal with fields which vanish at infinity. However, here, 
the surface term has the following form 



/ba 



(3.49) 



The presence of the spoils the argument that one can discard this surface term. In fact, to be correct, 
one would have to impose extra conditions on the fields to justify the dropping of the surface terms, 
i.e. the fields need to drop faster to zero than usually assumed. The partial integration is really needed 
to obtain the similarity with the GZ action. 



Keeping these criticisms in mind, we present a version of the Maggiore-Schaden construction that 
is based on local quantities, so partial integration may be done without surface terms at infinity, fol- 
lowing a method proposed in [157]. We shall consider the BRST symmetry of a local lagrangian 
density and derive a conserved Noether BRST-current and corresponding Ward identities. 

Consider the local Lagrangian density 

= Lym - 3^^'M° - a^c'^Df - avcpr^f 9f - Y'£>f 

-gi'foffabcD'/c^ - y'Df^'^" -fd{N^-\). (3.50) 

where s acts formally as in (3.41) and (3.42). We shall discuss the status of s shortly. Although the 
BRST operator s produces an x-dependent term slsP^ = -I- 'f'x^d"'', the Lagrangian density Z is 
nevertheless translation-invariant because C5°^ appears in the first line only in the combination d^^a^^ 
or &''^fabc- Because s is nilpotent, = 0, this Lagrangian density is ^-invariant, 

sZ = 0. (3.51) 



We now consider a local infinitesimal operator that acts on all fields <I>a according to 

d^^a{x)=e{x)s^a{x). (3.52) 

To control all partial integrations, we stipulate — and this is the main point — that the otherwise 
arbitrary function e(x) vanishes outside a small region. To be concrete, suppose that the fields satisfy 
periodic boundary conditions, <I>„(x+L) = <I>a(x), where are the edges of the volume, as happens at 
finite temperature. Then the global BRST operator s, satsifying ^©^^(x) = (p^^(x) +Y^x^8"'^ is not well 
defined because x^ is not periodic. This substantiates the criticism of [39] regarding the status of the 
global BRST operator s for, with periodic boundary conditions, the global operator s does not exist. 
However the local operator 8e = e{x)s is well defined, provided we stipulate that the support of e(x) is 
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restricted to a region contained within <x^<L^/2, and this will be sufficient for our purposes. 

(In practice we will choose z{x) = r|8(x — z), which is concentrated at the point x = z- Here r| is an 
infinitesimal constant.) With this restriction, the right hand side of (3.52) extends unambiguously to a 
function that is periodic in x^i with period L^^. Moreover Noether's theorem holds, 

= Z{x)sL{x) + a^8(x) S^a 

= a^8(x);^(x), (3.53) 
and yields the Noether current of the would-be BRST symmetry. 

The BRST-current j^,{x) inherits an explicit x-dependence from the transformation jc5°^ = (p^°^ + 
y^X/jd"^ given by 

= U/j~\~^/jy ^\ ; (3.55) 

where 

_yajyad^d _ ^ ^g^gai-^ + ^'f (D^a^b - gf^^'o'^fc'^ip'f^ (3.56) 

tpv = 'f{Df(4" - gr''<c^9f) (3.57) 

have no expUcit x-dependence. With periodic boundary conditions jij{x) is, strictly speaking, not well 
defined, because of this explicit x-dependence. However the product d^e{x) jij{x) is well-defined, and 
in our manipulations we shall always keep the well-defined product involving £(x) until, at the end of 
the day, 8(x) will drop out and we obtain identities satisfied by the quantities Mju(x) and tijv{x) that are 
manifestly consistent with translation invariance and periodicity. 

The renormaUzable action 5gz we shall consider is the integral of the Lagrangian density 

5gz = jd'^xl, (3.58) 

which is formally equivalent to Sgz- Because 8(x) vanishes outside a small region, a partial integration 
is justified in calculating the variation of the action and we obtain 

Se^GZ = ~ J^''^ ^pJp 

= -J d'^x 8 (x) (3pMp + + d^t/M Xv ) . (3.59) 



and 
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In classical physics the action is stationary under any infinitesimal variation along the classical path, 
so Se^GZ = for any e{x) and, choosing e{x) = ri8(;ic — z), we obtain, classically d^j^{z) = for 
—L/j/l < Zfj < Lfj/2. One may also show that, classically, the tensor t,jy is separately conserved 
^^tt^^v = 0, for it is the conserved Noether current of a symmetry generator /jy of the Lagrangian 
density 

a^py,! = (3.60) 

where a^p^ is a translation of by a constant, a^py^^f = and p^^a = 0, for / (5^^- Thus, 

classically, we obtain two identities that are translationally invariant and consistent with periodicity, 
^n^iM = 0, and B^m^ = —t/jfj. We note that is 5-exact, 

f^v=^Y'<cPv", (3.61) 
so U/j is classically a conserved current, modulo an j-exact piece, d/jU/j = —sy^D^^(p^". 

In quantum field theory a conserved current leads to a Ward identity. To obtain it we start from 
the expectation value 

{0)=N jd^O exp(-5Gz) , (3.62) 
and make the infinitesimal change of variable = <I>a + 5e<I>q, which gives the identity 

0=(5eO-5e5GzO) (3.63) 

or 

(8eO) = - j d'^xz{x) {d^j^{x)0) 

= -Jd'^x e(x) ( (a^Mp + + ^^,t^^ xv) O) . (3.64) 

To summarize so far: although the global BRST symmetry is not well defined when periodic bound- 
ary conditions are imposed, the corresponding local transformation involving the factor e{x) is well 
defined, and the Noether BRST-current of the Lagrangian density L is conserved classically, and the 
usual Ward identity holds for the Noether BRST-current provided that £{x) vanishes outside a suffi- 
ciently small region. 

As an application, we apply the Ward identity (3.64) for the case O = ro^'^(y). On the left-hand 
side of (3.64), the 5-exact quantity acquires a non-zero expectation value 



= fe(y)yv§'^^ (3.65) 

and we obtain 

£{y)y'yyd'' = -j d'' x z{x)[F^\x - y) + G'^^{x - y)xx] , (3.66) 

where 

Ff{x-y) = ({d^u^ + t^){x) <(y)\ (3.67) 
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G^^(x-y) = (d^t^xix) (^''{y)) . (3.68) 
We specialize to 8(_y) = ri5(_y — z) and find 

d{z-y)y'z,5'' = -Ff{z-y)-G'^i{z-y)zx , (3.69) 

or, with x = z — y, 

b{x)iz,h'='' = -F;^(x) - G'^^{x)zx . (3.70) 

This holds for all x and z. The expUcit dependence on z is consistent on both sides of this equation, 
and we obtain two identities, 

= -5(x)y'6^6^^ 

Ff{x) = 0, (3.71) 

both of which are consistent with translation invariance and periodic boundary conditions. The iden- 
tity involving G^(x) could have been derived using the Goldstone theorem for the symmetry (3.60), 
for which t/jy is (—7^ times) the Noether current, and which is spontaneously broken 

(p^rof ) = 8^v8''^ . (3.72) 



To solve the identity satisfied by G^^(x), we introduce the quantity 

= ^^((^-gr^^Dfc^cpfjCxXly)) , (3.73) 

so 

qvW = a.ri(..), (3.74) 

and the identity satisfied by G^{x) reads 

9;.ri(x) = -Y2S(x)5xv5^^ (3.75) 
We take fourier transforms and obtain 

^■^/i(^) = -7'5xv5^^ (3.76) 
We decompose F^^^ (k) into invariants, 

fi W = S'^'feSxvri +Mpvr2 +A:v5x^r3 + Wvr4) , (3.77) 

and the last identity yields 

ik^Fi = 

Y2 + T3+k^TA = . (3.78) 
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This is easily solved for Fi and r4 with the result, 



(3.79) 



The pole term in ^ is a Goldstone pole, and arises from the symmetry-breaking (^cOy^) / 0. Since 
the BRST current is fermionic, the Goldstone particle is a fermion. However because of the expUcit 
x-dependence, the pattern of symmetry-breaking is not the usual one, and from the Ward identity for 
the BRST-current j/j we obtain two identities corresponding to the two quantities U/^ and 



The formula we have obtained holds to all orders in perturbation theory. However we may check 
it at tree-level. To zeroth order in g we have 



which has the fourier transform 



k^ 



(3.80) 



(3.81) 



This is consistent with the exact result, and shows that there are no radiative corrections to the first 
term in (3.79). 



The identity F^'' = may be treated similarly. It reads 



(3.82) 



To solve it we observe that the right hand side is the quantity — r^^(^ — y)^ which we have just 
considered, and we define the propagator 



The the last identity now reads 



or, in terms of fourier transforms. 



a,A^(x-y) = -n^(x-y), 



(3.83) 



(3.84) 



ik^'^iik) 



= -h"^ky 



1 



where we have used our previous solution (3.79). We decompose A^(A:) into invariants 



A^(^:) = 6^^(6^Ai+^^^vA2) 



and obtain 



/(Ai+fc2^2) 



if- + {d-\)T2 



(3.85) 

(3.86) 
(3.87) 
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This is easily solved for Ai with the result 



cd iu\ get/ 



(3.88) 



This formula holds to all orders of perturbation theory and, as before, we may verify that the pole 
term is obtained in zeroth order of perturbation theory using 



A^^{x-y)\o = -^U;{x)W,'{y) 



(3.89) 



To conclude, we have seen that with periodic boundary conditions, the would-be BRST operator s 
is not globally well defined because of an explicit ;ic-dependence. However the local infinitesimal 
variation Sg = z{x)s is well-defined provided that s(x) is non-zero only in a small region, and this 
is sufficient for the existence of a classically conserved Noether BRST-current and corresponding 
Ward identity. A Goldstone pole results from the non-zero expectation-value of the 5-exact quantity 
5(5^^. These are characteristic features normally associated with spontaneous symmetry breaking. 
However the pattern of symmetry breaking is modified, and there are two independent Ward identities, 
corresponding to the two quantities and t^, both of which are consistent with translation invariance 
and periodicity in x. 

Other symmetries 

The Lagrangian density L, eq. (3.50), formally possess all the symmetries of 

Lq = 2|^o, (3.90) 
that it inherits from the formal change of variable (3.40), for we have 

2o((p,9,c,Z7) = 2((p',(p',c',Z7'). (3.91) 

Among these symmetries is the BRST symmetry that was discussed in the preceding paragraph. The 
other symmetries of Lq may be treated in the same way: the symmetries of Lq become would-be 
symmetries of L that are not defined globally. But for each such symmetry there is a conserved 
Noether current of L and corresponding Ward identity, each with its Goldstone boson or fermion.^ 
These will be discussed elsewhere. 



3.5 Restoring the BRST 

3.5.1 Adapting the BRST symmetry s is not possible 

One can ask whether it might be possible to modify the BRST operator, i.e. s — > Sm, in such a way that 
the new operator Sm would be still nilpotent, while defining an exact symmetry of the action, SmS = 0. 

In [39] a simple argument was presented, discarding such a possibility. Firstly, we recall that the BRST 
transformation (2.243) defines an exact symmetry of the action when y = as this corresponds to the 

^In [157] Goldstone bosons and fermions were shown to exist at the critical point ^ = selected by the horizon 
condition. In the present discussion we have not used this condition explicitly, and it appears that we have found Goldstone 
bosons and fermions for any non-zero value of Y- However it may be that a solution with non-zero y may exist only if the 
horizon condition is imposed. We are grateful to Martin Schaden for a discussion on this point. 
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physical situation in which the restriction to the Gribov region has not been implemented. Therefore, 
one should search for possible modifications of the BRST operator which depends on y, namely 

Sm = S + Sy, (3.92) 

whereby 

Sy = y-dependent terms , (3.93) 

so as to guarantee a smooth limit when y is set to zero. However, taking into account the fact that y 
has mass dimension one, that all auxihary fields ((p^'^^,(p][5'^, ro^'^, (0^^) have dimension one too, and that 
the BRST operator s does not alter the dimension of the fields^, it does not seem possible to introduce 
extra y-dependent terms in the BRST transformation of the fields {'^^ j^)^*^ ,(0^ , M]^'^) while preserving 
locality, Lorentz covariance as well as color group structure. However this argument does not exclude 
the local Maggiore-Schaden construction just introduced. 

3.5.2 Restoring the BRST by the introduction of new fields 

However there is another way to restore the BRST [122] following the idea that was initiated in [123] 
and [124]. In these papers it was shown that the broken BRST symmetry s can be rewritten as a 
non-local symmetry, being not nilpotent in [123] and nilpotent in [124]. In a subsequent paper, one 
succeeded in localizing the non-local BRST symmetry of [123] by the introduction of new fields. Let 
us summarize here these results. 

The non-local BRST symmetry s' 

First we review the results of [123], and write down a non-local BRST symmetry which is obeyed by 
the GZ action. We start from the standard GZ action (2.240), set gy^ = 6^, and drop the vacuum term 
for brevity because it does not influence any variation of the action, 

5gz = 5ym + J d^x [b"d,A'' + c%Dfc'^ + J d'x (cp^avDf cp^^ 

-m^^dyDffal' - gf''%(^'D^/c''<p^' - Q^f^'A^ (i^l' + 9^^) ) . (3.94) 

Following [123], we first drop the final term of the above expression, i.e. gf^'^d/jid^^D^'^c'^ip^^, so that 
there is no coc term in the action. This has no influence on the soft breaking of the action because the 
^-variation of the action we now consider, 

S'gz = SYM + ld'x{b"d^; + c%Dfc') + 

+ J d'x (^r^d^Dfip^' - m";d^D^^(o''; - e V'^a^ (< + cpj^) ) , (3.95) 

can be written as 

sS'cz = Ay- ld''xs(^gfakb{Dyc)%(iS";^>'';) (3.96) 

= |d^x(-eVD^T''(<+9f)-eV"''A>^'+gr''-(^^^^ 

^It is understood that the dimensions are assigned to the fields A^,b",c",c" as in Table 2. 1 of sec. 2.6. It is apparent that 
the BRST operator s does not alter the dimension of the fields. 
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where Ay can be found in equation (3.31). Therefore, we can indeed safely leave out the term 
^yafecg^fjjae^j^rf^djpce After the localization however, this term shall be picked up again. Accord- 
ing to [123], the positivity of the Faddeev -Popov operator inside the Gribov region allows to rewrite 
(3.96) as 

= I d'x |^(z)rAe)[(avDv)-i]- J^^Gz+eV'''^'A;[(av£»v)-Y'"g|^5Gzj , 0.97) 

with 

From (3.97), we can now read off a new nonlocal BRST symmetry, s'Sqz = 0, generated by 

s'A"^ = -Dfc\ /c" = ^gr^'cV, s'(f = ^7«-(DfA^)[(av£)v)-Y", ^'b" = 0, 
s'(p^ = (sf;, s'(a";' = 0, s'm";' = + Q^[{D^d^)-^YPf''^'Af,, s'^"^' = 0. (3.99) 

We repeat that this symmetry s' is not nilpotent, j'^ ^ 0, see also [123]. 
The local BRST symmetry 59 

In [122], this non-local BRST symmetry s' was localized into sq. This was done by adding an auxiliary 
part to the GZ action, 

S'^z = Sgz + S,u. (3.100) 
with 5gz given in (2.240) and 5aux given by 

5aux = / d^'x {a"d^Dflt + a"d^Dfal'+gf^'{d^(f){Dl^c'^)a' - a«£»f aJ aJ) 

+ J d^x(^R"d^D;'t + Q"d^D"J'Q' + gf'%Q"D'/c'r -q'k^ , (3.101) 

where 

A" = /«^^((p^'- + (p^^') 

= D;;^(gr^^'a^p-o . (3.102) 

Here we have introduced of a new set of doublets (a«,a«), (n,a°), (Pf .^Pf), (2",/?°), 

(^^r)i.e. 





= a", 


sQ." 


= 0, 






ja° = 0, 








= 0, 










= R", 


sR" 


= 0, 






.r=o. 
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a«, a°, , , and 2° are bosonic fields, while the other fields H", O", ^f, , and R" are 
anti-commuting fields. It can be checked that the action 5q| is invariant under the following BRST 
transformation sq 





= - {D,cr , 


SBC'' 




SBc" 






= 0, 








= 0, 








= 0, 




= a", 




= 0, 




= a^ 


sea" 


= eV, 




— vuab 




= 0, 












= R\ 


seR'^ 


= 0, 








= eV. (3.104) 



Again, this symmetry sq is not nilpotent, Sq ^ 0. However, one can check that Sq = 0. 



Let us summarize some properties of the action Sq^,- 

• Firstly, we need to ask what happens if we set 9^ = 0, or equivalently, = 0. In this case, 
we are not considering the restriction to the Gribov horizon anymore and we are back in the 
standard Yang-Mills framework. In order for the theory to be meaningful, we would expect that 
5Q2|e2=o equivalent with the ordinary Yang-Mills action 5ym + 5gf. This is indeed the case, 
as is proven in [122]. 

• Secondly, one should also prove that the GZ action and the new action Sq^ ^e equivalent. In 
other words, for (|) G {A^, b", c", c", (^f, ^f, (af, (af], we require the following identification 

(^(xi) . . .^(^„))ioe = j [d<I>]ioc^(^i) • . .^(;c„)e-^Gi = (^(;ci) . .Mxn))Gz ■ (3-105) 

This was also proven in detail in [122]. 

• Thirdly, in order for the action 5^1 to be meaningful at the quantum level, one should prove 
that this action is renormalizable. However this has not yet been proven. In [122] many other 
symmetries were given which could be useful for the renormalization of the action. It is worth 
mentioning that s% not being nilpotent, is not necessarily a problem to prove renormalizability. 
From the result Sq = 0, it should be possible to construct a nilpotent symmetry, see [125]. 

3.5.3 Two different phases of the restored BRST symmetry 

We investigate the BRST operators s' and sq just defined. They are closely related. Indeed, in (3.99), 

we have ^'ro^^ = (p'^^ + B'^[(D^d^)-^]"PfPi'Af„ while in (3.104), sq acts according to sqW"^' = (p^ + 
O^P^*^, where P^'^ satisfies the equation of motion 

Dfa^P^<^ = r*^A^, (3.106) 

— CIC 

obtained by (left) differentiating (3.101) with respect to P^ . This equation has the formal solution 

p- = (D^a^)-i]«''r^^A?, (3.107) 
and with it, the action of sq agrees with the action of s'. 
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We shall arrive at the unexpected conclusion that there are two different solutions Pv^'**^ and P^,^^**^ 
to equation (3.106) which provide two different phases of s' and .ve. 

Solution 1: For simplicity we suppose that is transverse, dfjA^j = 0, so the Faddeev-Popov 
operator Dp(A)9p is hermitian. We quantize in a finite Euchdean periodic box of edge L, so the 
spectrum of £)^(A)3^ is real and discrete, and there is a complete set of eigenvectors, 

<(A)a^V^: = ^.<- (3.108) 

For generic A not precisely on the Gribov horizon, as we shall suppose, the null space of D^(A)3^ 
consists of constant functions 3^\J/ = 0. To obtain solution 1, we use the eigenfunction expansion of 
the inverse operator, 

[iD,d,rr\^,y) = (3.109) 
where the prime means sum over non-zero eigenvalues. Solution 1 is thus 

= l^f^'y^i{y)f'"'Al{y). (3.110) 

As an explicit example of solution 1 in SU (2) gauge theory, take A^ = ccos(fcci)5v25^^, where 
k = Inm/L, m is a non-zero integer, and xi is the first componant of x^, which is transverse, = 

c5^^32 co?,{kx\ ) = 0. It is easy to verify that Mf" = cos(^xi )m" , where u" is any constant color vector, is 
an eigenfunction. Indeed we have A^3^ cos (/:.«i) = ccos(/:.«i)8''^32Cos(fcxi) = 0, soD^^3^\j/^ = 3^\j/' = 

—k^'\^. For this choice of Ay, the eigenfunction expansion of P^'' is given by Py'^ = (— /:^)~^ccos(fcci)/"^'^ 
(-A:2)-i/«i«^A^. 

Solution 2: The simple identity, 

D^^(A)3,5''% = gr^^A^, (3.111) 
reveals unexpectedly that a second solution of (3.106), for any Ay (x), is given by'^ 

Remarkably, it is independent of A^. Solutions 1 and 2 are different, as one sees by comparison 
with the explicit example of solution 1. Note that a perturbative expansion yields solution 1, whereas 
solution 2 is non-perturbative. 

With solution 2, the BRST transformations (3.99), s^J^^' = + e2[(Dvav)"']"^/^'^'A^ sW^' and 
s^Tsfl^ = cp^'' + 0^P|^'' coincide with the BRST transformation of Maggiore and Schaden, = (p^"'' + 
'f-X/jd"'-', after obvious relabeling, and with 6^ = gy^. One can reject solution 2 by imposing periodic 
boundary conditions, in which case is not well defined. However even with periodic boundary con- 
ditions, the local infinitesimal variation 5 = e{x)s', or 5 = e{x)sQ, remains well-defined for solution 2 
provided that e{x) vanishes outside a sufficiently small region. As shown in sec. 3.4.3, consideration 
of this variation is sufficient to establish the existence of a conserved Noether current and correspond- 
ing Ward identity, and the presence of the Goldstone pole normally associated with a spontaneously 

"The horizon function H{A) = ~f d''xgf'"''A'/, (x) [{D^d^)- ' )'"'g/**"'A*] {x) was derived in (2. 178) by solving the eigen- 
value problem (by summing the perturbation series). Consequently solution 1, the eigenfunction expansion, should be used 
to evaluate if (A). 
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broken symmetry. On an infinite volume, the system may choose the second solution (3. 112), in which 
case we have 

= e^g-^x^S"*. (3.113) 

In this case sq is spontaneously broken because the vacuum state is not ^e-invariant, 

se'Vo^O. (3.114) 

In conclusion we note that the action Sq^ = Sqz + "^aux, is local and possesses a localized BRST 
operator sq. Its field equation possess two solutions. Solution 1 preserves BRST symmetry, whereas 
solution 2 spontaneously breaks this symmetry. At this point we do not know which solution the 
system will choose. If the correlators of the system described by the action 5q| coincide with the 
correlators described by Sqz (for the fields appearing in ^gz), then according to the calculation of sec. 
3.4.3, BRST will in fact be spontaneously broken. However Sq^ appears to have two phases: one in 
which BRST is preserved and one in which it is spontaneously broken. 

Interpreting the restored BRST symmetry 

BRST symmetry is thought to encode the geometrical character of a quantum gauge field theory. It 
is therefore encouraging that the local action Sq^ possesses a BRST symmetry — as does Sqz in 
the local sense discussed in sec. 3.4.3 — even though this theory is presumably only approximate 
(because it was derived by integrating over the Gribov region instead of the fundamental modular 
region).^ If BRST is spontaneously broken, as indicated by the calculation of sec. 3.4.3, the vacuum 
itself is not BRST-invariant, s^'q 7^ 0, eq. (3.1 14), and we cannot identify the positive metric subspace 
by the usual condition = 0. Here we enter into an unfamiliar realm in which it is not known 
how to define a positive-metric physical subspace. However, to the authors' knowledge, proof of 
positivity of the subspace i^*? = relies heavily on the existence of an asymptotic field for every 
fundamental field, including the gluon and quark fields. But the very meaning of confinement is that 
the asymptotic gluon and quark fields do not exist in the confined phase of QCD. Thus at present 
there is a general lack of understanding of how to construct a positive-metric physical subspace in 
the confined phase. Making a virtue of necessity, we are tempted to interpret the difficulty we have 
encountered, of finding the positive-metric physical subspace when BRST is spontaneously broken, 
as a welcome signal that gluons and quarks are not in the physical subspace of the phase we have 
found. This suggests a scenario in which spontaneous breaking of BRST symmetry serves as an order 
parameter of the confined phase. In this interpretation, the two solutions of the equations of motion 
of that we have found, in which BRST is either preserved or spontaneously broken, correspond 
respectively to the deconfined or confined phases of QCD. 

^The i-invariance of Sqz for any value of y is in fact more than is needed for physics. Indeed the theory defined by 
Sqz need not be a gauge theory unless the horizon condition holds, so it would be sufficient for physics if the action were 
BRST-symmetric only for that special value of y determined by the horizon condition. Instead Sqz is BRST-invariant for 
every value of y. 
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3.6 The GZ action and its relation to the Kugo-Ojima confinement cri- 
terion 



3.6.1 Introduction: the Kugo-Ojima criterion 

Some important results concerning a possible origin of confinement, were given in [7, 8]. In these 
papers, confinement was related to the enhancement of the ghost propagator. 

Let us explain this a bit more in detail. The starting point of the analysis of Kugo and Ojima is a 
well defined nilpotent (BRST) symmetry operator Qb and a ghost charge. Their analysis yields two 
results. Firstly they showed that with a well defined nilpotent symmetry, all unphysical states, see 
section 2. 1 .6, form so-called quartets [7] that decouple from the physical spectrum. In this way, only 
physical states, which are closed under the symmetry Qb but not exact, survive. In this way they 
proved that the longitudinal and temporal gauge polarization, the ghost and the antighost fields can be 
excluded from the physical spectrum. In fact this idea is very general, and can be applied whenever a 
system admits a nilpotent symmetry s. 

Secondly, for the Faddeev-Popov action, they also showed the following. Using the equation of mo- 
tion for the gluon field, the conserved global color current can be written as 

J"^=df,F^ + {QB,Dfc''}, (3.1 15) 

and so the color charge, which is the integrated zero component of 7°, is given by 

= I (diFo^i + {Qb,D"o'A) ■ (3.116) 

Thus there are two criteria which need to be satisfied in order to have color confinement, namely 
Q" = for all physical states. The first is that the gluon propagator should not have massless poles, 
so the first term of the last expression vanishes because it is the integral of a derivative^. The second 
criterion is that {Qb,Dq'c''} should be well defined, which is the case when 

m(0) = -1, (3.117) 

with u{p^) defined through the following Green function^, 

I d'xd"^ (Dfc\x)Dl^c%Q))^ ={{^^- ^) u{p^) -^Y"'. (3.118) 

(0)pp stands for the expectation value taken with the Faddeev-Popov action. If the two criteria are 
met, then is well defined and we have (cp| | V)phys ~ ^' color confinement is guaranteed. 

The second criterion can be connected to the ghost propagator. Indeed, in [8], it was shown that 
one can parameterize the ghost propagator, defined as follows 

(c%-p)c\p)) =5«^G(p2), (3.119) 

^If the gluon propagator has massless poles, the first term is ill-defined. 

^Strictly speaking, the KO analysis is done in Minkowski space. We shall however, as any functional or lattice approach, 
consider the corresponding operator in EucUdean space. 
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in terms of^ 



G{p^) = ^7, r4^ • (3-120) 

This relation was also discussed in [22, 20, 23, 24, 25]. It is usually assumed that^ w{p^) = 0, so that 
M = — 1 implies an enhanced ghost propagator. Notice that this scenario is exactly predicted by the 
GZ framework. 



3.6.2 Remarks on the KO criteria 

Two comments should be made concerning the argument of Kugo and Ojima. First of all, in the KO 
framework [7,8], the existence of a globally well-defined BRST charge is assumed. This assumption 
has now been substantiated by the construction of a BRST operator in lattice gauge theory which does 
have non-perturbative validity [126]. On the other hand, the derivation of the two criteria was done 
by employing the usual Faddeev-Popov gauge-fixed action, which does not take into account Gribov 
copies. 

Moreover the KO argument does not really hold for the GZ action, as (1) the GZ action breaks the 
usual BRST symmetry, either explicitly, see equation (3.31), or spontaneously, as in the Maggiore- 
Schaden construction discussed above, and (2) the GZ action does take into account Gribov copies. 

Another point which should be mentioned is that the Kugo and Ojima criterion (3.117) is sometimes 
useful as a starting point in looking for a solution to the Dyson-Schwinger or functional renormalization- 
group equations, see e.g. [26]. 



3.6.3 Imposing u(0) = -1 as a boundary condition in the Faddeev-Popov action 

In [27], the following was tried. The authors imposed the constraint m(0) = — 1 directly into the 
Faddeev-Popov theory, by appropriately modifying the measure one starts from. By using thermo- 
dynamic arguments, as iat the end of sec. 2.5.2, it was concluded that one precisely finds the GZ 
action! The boundary condition m(0) = — 1 is exactly encoded in the horizon condition for the y, see 
expression (2.208). The details can be found in [27, 1 17]. 



3.6.4 Conclusion 

We can thus conclude: 

• Imposing a boundary condition into a theory can have serious consequences. Here for example, 
we have seen that imposing the boundary condition m(0) = — 1 in the Faddeev-Popov measure 
with BRST symmetry, leads us to the GZ action without BRST symmetry. Therefore, imposing 
a boundary condition can change the symmetry content of a theory. 

• The relevance of the KO criterion in the GZ formalism is not clear because the KO analysis was 
done in the Faddeev-Popov formalism. Due to the breaking of the BRST symmetry in the GZ 
action, one cannot simply redo the KO analysis in the GZ framework. 

^Actually, in [8], another notation v(p^) has been used instead of w(p^), the relation being w{p^) = p^v{p^). 
^w{p^) = has been checked up to two loops, see [115]. 
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3.7 The relation of the GZ action to the lattice data 



3.7.1 The lattice data 

It would be interesting to verify the GZ theory. As we are working with a quenched Yang-Mills theory, 
we cannot use experimental data, but we can compare our analytical results with lattice calculations. 
Two particular quantities have been tested in great detail, namely the ghost and the gluon propagator, 
as they are believed to play an important role in confinement scenarios. The first calculations of the 
gluon propagator in 4 dimensions in the Landau gauge were already carried out in 1987, see e.g. [127] 
on very small lattices of the order of 4^ x 8, and of the ghost propagator around 1996, see [128] on 
lattices of the order of 32^. Ever since, many papers appeared on the subject. Let us review them 
briefly. 

The 4d and 3d studies 

Let us start with the gluon and the ghost propagator in 4 and 3 dimensions. There have been many 
investigations of them, see [132, 133, 134, 135, 136, 137, 129, 138, 139, 140, 130, 131, 50, 141, 142, 
143, 144, 145, 146, 29, 33, 28, 147, 42, 34, 37, 148, 149] for a selection of some papers. A useful 
overview of the lattice results can be found in [37], where the results are examined in the Ught of the 
GZ framework. 

Concerning the gluon propagator^^ in 4d in the infrared, it quickly became clear that perturbative 
behavior was not seen, and that the gluon propagator is not enhanced at small momentum. Instead a 
massive behavior was seen [137, 138, 130, 141, 143]. Further studies confirmed this behavior; how- 
ever, in order to reach the deep infrared, larger lattices were required [152, 146]. Therefore, one can 
say that until 2007 not so much was known about the gluon propagator in 4d, although the hope was 
that it would vanish at zero momentum in order to agree with the GZ framework. We must say, many 
studies already reported a finite gluon propagator at zero momentum [130, 144] although the volumes 
were not large enough to reach the far infrared [145]. Finally, in 2007, the papers [29, 33] appeared, 
with simulations on huge lattices^ ^ From these papers, it was found that the gluon propagator in 
4d did not vanish at zero momentum! After these papers, more papers appeared [28, 34, 148, 42], 
which arrived at the same conclusion, i.e. the gluon propagator seems to attain a finite value at zero 
momentum. Some papers were still in agreement with a vanishing gluon propagator, see e.g. [147]. 

Fewer studies of the gluon propagator have been performed in 3d. The first one was in 1999 by Cuc- 
chieri [129], where a gluon propagator was reported which is infrared suppressed. The deep infrared 
was not accessible at that time. Later papers confirmed this suppressed behavior, [139, 140, 153]. It 
was however not clear what happened at zero momentum, but at least in 3d the gluon propagator was 
not in contradiction with a vanishing gluon propagator [140]. Positivity violation was akeady clearly 
seen, see [50]. The hope was that at larger lattices, a turnover would be seen so the gluon propagator 
would vanish at zero momentum. Then the paper [29] appeared, with simulations on huge lattices, 
where a finite gluon propagator a.tp = 0m3d was reported. Paper [28] supported this result. The up- 
shot in 3d is seen in [37] where on lattices of size 320^ the gluon propagator exhibits a clear turnover 
in the infrared namely, D{p) has a maximum below which it decreases as p decreases, and approaches 

^"The propagators were investigated both in SU{2) as in SU (3); for a comparison between these two settings, see [150, 
151]. 

'^It should be mentioned that the lattice spacing is also very large, in the range of 0.18 — 0.22 fm. 
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a finite value D{0) / from above. The decrease of D{p) for /? — 0, has no other explanation than 
the proximity of the Gribov horizon in infrared directions, in accordance with the Gribov scenario, 
but the finite value of D(0) / contradicts perturbative calculations with the GZ action. By contrast, 
in Ad on 128"^ lattices [37] there is something like a shoulder at very low momenta as the finite value 
of Z)(0) is approached. 

The ghost propagator was more complicated to study on the lattice because one has to access the 
inverse of the Faddeev-Popov operator. In Ad, papers [128, 92, 131, 154, 142] showed a behavior 
between the zeroth-order perturbative behavior X/p^ and the \/p^ singularity predicted by the GZ 
framework. Many problems however needed to be overcome [37], including finite-volume effects, 
discretization effects and Gribov-copy effects, which could influence the results obtained to far. The 
hope was that at larger volumes indeed the enhancement of the ghost propagator would be more clearly 
seen. The breakthrough came in 2007, whereby in the same papers [29, 33, 34] as for the gluon propa- 
gator, the ghost propagator was investigated at large lattices up to 128^. The result was very surprising, 
as the enhancement completely disappeared and the \/p^ behavior was recovered in the deep infrared. 

The ghost propagator in 3J was even less studied [153, 29, 30]. The main conclusion is however 

again the same. The larger the volumes, the less enhancement of the ghost propagator is seen. Finally 
a consensus was reached, and one finds again the \/p^ behavior at very large lattices. 

From these lattice results one may conclude that something is missing in the current GZ framework. 
Much research has been done on this problem, see e.g. papers [38, 40, 39, 41, 155], where a dynami- 
cally refined GZ action was proposed in 3 and 4 dimensions. 

The 2d studies 

2d appeared to be different from the 3 J and Ad case. Not so many studies were carried out in 2d, but 
they are all consistent [31, 28, 156]. The gluon propagator goes to zero at zero momentum; this was 
even checked with lattices up to 2560^ [156], while the ghost propagator displays enhanced behavior. 
Therefore, the lattice data in 2d are at least qualitatively in agreement with the GZ framework. How- 
ever, two comments are in order. Firstly, the ghost propagator does not exactly display enhancement 
like 1 / p^, as the following enhanced behavior was reported: consistent with Mm^^Q Q{p^) ~ 
see [31], and in agreement with Schwinger-Dyson calculations [46]. Secondly, a qualitative compari- 
son has not been done so far, and would be needed before drawing any further conclusions. 

Let us also mention that there has been done some work on the GZ action in two dimensions, whereby 
it was shown that the dynamical refinement is not possible [41], therefore meaning that the GZ frame- 
work would still be valid in 2d. This qualitatively agrees with the lattice data so far. 

3.7.2 Reflection positivity 

It was pointed in [158] that if the gluon propagator D{k) vanishes at A: = 0, then reflection positivity 
would be violated. In fact it would be maximally violated, as we now explain. 

Reflection positivity and the other axioms of Euclidean quantum field theory imply that D{k) 
satisfies the Kallen-Lehmann representation, 




(3.121) 
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where p(m^) > is a positive weight. This imphes that the gluon propagator in x-space, 

D(x)= [ dm^ p(m^)DM), (3.122) 
Jo 

is positive for all x, D{x) > 0, because the free propagator of a particle of mass m is positive for all x, 
as one sees from 

/»oo 

DJx)=2~^ dailnay^/^expll-^m^a + x^a-^MyO. (3.123) 
Jo 

If the zero-momentum propagator vanishes, 

D{0) = J d''xD{x) = 0, (3.124) 

and if D{x) is positive (or zero) for all x, then the gluon propagator vanishes identically, D{x) = 
for all X, which is false. Thus, if the gluon propagator vanishes at A: = 0, then reflection positivity is 
violated. In fact it is maximally violated in the sense that if D{Q) = 0, then the gluon correlator D(x) 
is positive and negative in equal measure. The Gribov propagator D{k) = /[{k^Y + m^\ maximally 
violates positivity, because Z)(0) = 0, and the Kallen-Lehmann representation is violated by having 
poles at imaginary mass-square, k'^ = -izint^. 

As discussed in the preceding section, numerical investigation indicates that D{Q) vanishes in Eu- 
clidean dimension 2, whereas it is finite, D(0) > 0, in Euchdean dimension 3 and 4, so reflection 
positivity is maximally violated for d = 2, but it is not maximally violated for <i = 3 and d = 4. 
Nevertheless, numerical investigation does show clear violation of reflection positivity in Euclidean 
dimension 3 [50] and in dimension 4, for both quenched and unquenched cases [51]. This evidence of 
violation of reflection positivity in dimension d = 2,3 and 4 stands in contrast to perturbation theory, 
where BRST symmetry assures that the perturbative gluon spectrum is physical and unitary. Viola- 
tion of reflection positivity is a signal that the gluon is not a physical particle, in accordance with the 
confinement scenario. 



3.8 Color confinement 

In this section we present an exact bound that results from the restriction of the functional integral to 
the Gribov region which implies that the color degree of freedom is confined [158, 159]. 

3.8.1 Exact bounds on free energy 

Let the free energy W{J) be defined by 

exp[W{J)]= f dAp{A)exp{J,A), (3.125) 
Ja 

where only transverse configurations, 9^A^ = 0, that lie inside the Gribov region Q., are integrated 
over, and the source term is given by {J, A) = f d^x J^{x)A'^{x). Since A is identically transverse, only 
the transverse part of the source J contributes, and we impose that it is also identically transverse, 
= 0. The result is quite general, for we shall suppose only that p(A) is a non-negative weight, 
p(A) > 0, that it is normahzed, f dA p{A) = I, with support restricted to the Gribov region, p(A) = 
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for A ^ n. Since the fundamental modular region A is contained in H, the result holds in particular if 
the integral is restricted to A. 



We specialize the source to a single fourier component, 

J^{x)=H^cos{kxi), (3.126) 

where we have aligned the 1-axis along k, and transversality of J/j reads Hf = 0, so the source term is 
given by 

{J,A) = J d'^xH^cos{kxi)Af{x). (3.127) 

where / = 2, ....D. It may be interpreted by analogy with magnetic spin systems, where A^{x) is the 
analog of a spin variable, but with a color index a, and is the analog of an external magnetic field, 
also with a color index, that is modulated by a plane wave cos{kxi). (The external color-magnetic field 
should not be confused with the Yang-Mills color-magnetic field tensor F-j.) In the limit A: — > 0, 
the source becomes a constant color-magnetic field. The free energy W{J) now depends only on the 
parameters k and H", and we write 

W{J) = Wk{H). (3.128) 



This speciahzed free energy is sufficient to calculate the gluon propagator, 

^"%m = 2^^^,, (3.129) 



I J 

where k is aUgned along the 1-axis, and i,j = 2, ...D. Here Wk{H) is the free energy per unit volume, 

w,(//) = ^, (3.130) 

where V is the Euclidean volume. 

An exact bound for Wk{H) on a finite lattice was given in [158, 159] which holds for any numeri- 
cal gauge fixing to Landau gauge, with support restricted to the Gribov region Q., which is to say, with 
support on relative or absolute minima of the gauge fixing functional. In the continuum Umit and at 
large Euclidean volume V, the lattice bound implies the bound in D Euclidean dimensions, 

Wk{H) < {2DYI\\HI (3.131) 

where |//p = i:^_,(//^^)2. 

More recently, a stricter continuum bound for Wk{H) at finite H was obtained [160, 161], that also 
holds for any numerical gauge fixing with support inside the Gribov region H, 

Wk{H) < 2-i/2ittr[(//"//")i/2]. (3.132) 

Here W^H" is the matrix with elements Y.a^t^'j- It positive eigenvalues, and the positive square 
root is understood. This bound is in fact optimal for a probability distribution p(A) of which it is 
known only that its support Ues inside the Gribov region. Either bound yields in the zero-momentum 
limit 

woiH) = XxmwkiH) = 0. (3.133) 
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These bounds imply that the magnetization vanishes in the limit ^ — 0. To show this we write = 
he"^, where h represents the magnitude of the external magnetic field in the direction of the fixed unit 
vector e^, with Ejua(^P^ = 1' ^nd we write Wk{H) = Wk{h). The magnetization is given by 



mk{h) = = (^~'l ^''^ cos(fcci)e? Af{x) )h. (3.134) 



dh 

It vanishes at zero magnetic field, mk{0) = 0, if there is no spontaneous magnetization, as we suppose, 
and increases monotonically, '^"^^ > as follows from 

^ = ^ = (l/2)D(/.)>0, (3.135) 

which is positive because the gluon propagator D{h) is positive at any h. The free energy may be 
expressed in terms of the magnetization and the bound reads, 

^ fo'dh'Mh') ^ tr[(.V) V^] (3.136) 

h h 

by (3.132). Because mk{h) increases monotonically, this bound implies that mk{h) approaches a finite 
limit, ]imh-^oomk{h) = mk{°°) which satisfies 

m,(oo) = hm fodh'mkih') ^ 2-1/2^ tr[(gV)i/2]. (3.137) 

Thus the magnetization at finite h satisfies the bounds, 

< mk{h) < MkH < ti-[(eV)i/2]. (3.138) 

We arrive at the remarkable conclusion that in the limit of constant external magnetic field, ^ — 0, the 
color magnetization per unit volume vanishes, no matter how strong the external magnetic field h, 

limntkih) = {V-^ f d^x < A?(x) exp[ h f d^x e^ A?(x) ] ) = 0. (3.139) 

k^Q Jv Jv 

Thus the system does not respond to a constant external color-magnetic field. In this precise sense, 
the color degree of freedom may be said to be confined. This results from the proximity of the Gribov 
horizon in infrared directions. 



If Wk{H) were analytic in H in the limit ^ — > 0, the bound, \mik^()Wk{H) = 0, would imply that 
all derivatives of the generating function wq{H) vanish, including in particular the gluon propagator 
D{k) (3.129) at A; = 0. In [158, 159] it was assumed that Wk{H) is analytic in 77 at A: = 0, and it was 
concluded that the gluon propagator at A: = vanishes, D(0) = 0. However, as reported in the preced- 
ing section, although this agrees with lattice data in Eucldean dimension 2, it disagrees with recent 
lattice data which indicate a finite value, D{0) / 0, in Euclidean dimensions 3 and 4. If this is true, 
then Wk{H) must become non-analytic in H in the limit A: — > 0. Two models were recently exhibited 
[160, 162] for which the free energy Wk{H) vanishes in the limit \imk-i.oWk{H) = for all H — and 
thus exhibit confinement of color — but for which the gluon propagator D{k) is nevertheless finite at 
A: = 0, in agreement with the numerical data in 3 and 4 Euclidean dimensions. 
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3.8.2 Simple model 

We exhibit here a similar model, defined by 



r / k^H'^H''\ 1/2 

Wk,moi{H) = tr [ [g\k) + -^—) - gik) 



(3.140) 



where g{k) > is a function at our disposal, and H'^H" is the matrix with elements HfH", for 
i,j = 2, ... D. This model possesses the following features of the exact theory [160]: (i) It satisfies 
W/t.mod(0) = 0, which is necessary at // = for a normalized probability distribution /^/A p(A) = 1. 
(ii) It satisfies the bound (3.132), 



Wk,r.oa{H)<2-'l^kir[{H-H'^Yl\ 



(iii) The magnetization 



dh 



tr 



(3.141) 



(3.142) 



where Hf = he", satisfies the bound (3.138). 
(iv) The gluon propagator 



{\/2)D{k,h)=\x 



'g {k) + ^ ^ 



>0. 



(3.143) 



is positive at all h, as required for it to be a covariance. This assures that the magnetization is mono- 
tonically increasing. 



(l/2)D(^,/.) = ^^^^^f^>0. 

ah 



(3.144) 



In this model, the free energy Wk^raod{H) and magnetization mk^raod{H) both vanish in the limit A: — > 
for all h, in agreement with the exact bounds given above, and consistent with confinement of color. 

If Wk.mod(H) were analytic in H at k = 0, then the gluon propagator (1 /2)D{k,h) = ^ '^^^2''^^'' would 
vanish at A: = for all h. However from (3.143) we find for the gluon propagator, D{k) = D{k,h = 0), 




(3.145) 



Whether or not it vanishes at A = depends on the behavior at A = of the function at our disposal, 
g{k). If g{k) = np-k^ at low momentum, where m is some mass, then the model yields at low momen- 
tum 

limD(A) = \, (3.146) 

which is finite, and agrees with numerical results in 3 and 4 Euclidean dimensions. On the other 
hand, if g{k) does not vanish as rapidly as A^ at A: = 0, then we get \\mk^QD{k) = 0, which agrees 
with the numerical data in 2 Euclidean dimensions. Thus we have a model for which the free energy 
Wk{h) vanishes at A = for all h, in agreement with the exact bound provided by the Gribov horizon. 
Nevertheless, the gluon propagator D{k) may remain finite at /i = 0. 

^^The first two terms in the expansion of (3.140) in powers of H are given by wi^jj^„^{H) — ^^^^HfHf + 

g{k)0[k'^ H'^ / g'^ {k)\. For g{k) = m^k^, which makes the first term finite at A: = 0, the coefficient of the H''' term is of 
order 1 /k^, which diverges ask-^Q, as do all higher order coefficients. This exhibits the non-analyticity of w/t.modl^) in 
Hatfe = 0. 
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Chapter 4 

Overview of various approaches to the 
Gribov problem 

In this section we present a critique of the present approach and some related analytic approaches, 
particularly stochastic quantization, with a view toward lessons learned, open problems, and direc- 
tions for future research. 



4.1 GZ action 

In the present review we have gathered in one place the derivation and properties of the GZ action, 
eq. (2.213), which otherwise are widely scattered in the literature. This action, in Landau gauge, 
incorporates a cut-off of the functional integral outside the Gribov region Q. which is done to avoid 
counting Gribov copies that are related by "large" gauge transformations. It must be said at the outset 
that, as explained in sect. 2.2.2, there is an approximation involved in the derivation of the GZ action, 
in which the fundamental modular region A is replaced by the Gribov region H (although it has been 
conjectured that they may give the same expectation values [108, 110]). The GZ action is of interest 
nevertheless because it is local and renormaUzable even though the cut-off is non-perturbative. Its 
non-perturbative character is manifested by the appearance in it of the Gribov mass 

If one calculates with the GZ action instead of the Faddeev-Popov action, and expands in powers 
of the coupling constant g, one obtains an alternative, perturbatively renormaUzable series, in which 
the zeroth-order gluon propagator is the Gribov propagator, DQ{k) = jj^^^. Its poles occur at the 

unphysical locations = ±iy^. This does not correspond to a physical particle but is appropriate for a 
propagator of gluons that are confined. This would clearly be a non-perturbative phenomenon in con- 
ventional Faddeev-Popov perturbation theory, but it appears already in zeroth order in the alternative 
perturbative expansion provided by the GZ action. In this alternative perturbation series, the gluon 
propagator D{k) vanishes atk — 0, D{0) = 0. As discussed in sect. 3.7.2, this corresponds to maximal 
violation of reflection positivity, so the Kallen-Lehmann representation is maximally violated, and 
this is manifested concretely in the unphysical poles of the gluon propagator. Lattice data also show 

'yis not an independent parameter; its value is fixed in terms of Aqcd by the horizon condition ^ = 0, where F is the 
quantum effective action. 
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clear violation of reflection positivity [50, 51], and we take unphysical singularities of correlators of 
gauge-non-invariant fields to be a correct description of confinement in continuum QCD. 

Limitations of the GZ approach in perturbation calculations are mentioned in the Introduction. They 
may possibly be remedied by non-perturbative calculations. 

4.2 Refined GZ action 

In response to this challenge, allowable condensates which preserve renormalizability, such as A^, 
have been treated in a mean-field approach, by introducing them into the GZ-action. This results in 

what is known as the refined GZ action (RGZ) that has been studied extensively [38, 39, 40, 41, 155]. 
By this method, satisfactory agreement with lattice data for the gluon and ghost propagators has been 
achieved, and the glue-ball spectrum has been calculated. A separate review article would be required 
to provide an overview of these developments. 



4.3 Coulomb gauge 

The present review has been devoted to non-perturbative QCD dynamics in the Euclidean Landau 
gauge. However the Coulomb gauge also requires a cut-off to avoid Gribov copies, as Gribov pointed 
out in his original article [10]. The cut-off in the Coulomb gauge has also been implemented by a 
local action, along the hnes of the GZ action in Landau gauge [68]. QCD calculations in the Coulomb 
gauge have also been done by variational calculation in the Hamiltonian formulation [163, 164]. An 
exact bound exists in the Coulomb gauge, that goes by the slogan "no confinement without Coulomb 
confinement." Thus if the (gauge-invariant) Wilson potential V{r) is confining, then the temporal 
gluon propagator in Coulomb gauge, Doo(r,f), has an instantaneous part, ycoui('")5(0' that is also con- 
fining, [165]. The interest of this bound is that the Wilson potential, obtained from gauge-invariant 
Wilson loop P e\p{i f Afjdx'^) , involves «-point functions of all orders n, whereas the gluon propagator 
Doo{r,t) is a 2-point function. A proof of renormalizability in the Coulomb gauge remains on open 
challenge. 



4.4 Dyson-Scliwinger equation 

The Dyson-Schwinger (DS) equation is perhaps the most highly explored analytic approach to non- 
perturbatve calculations in QCD [26, 45, 46, 47, 48, 49, 44]. It is similar in spirit to the functional 
renormalization-group equation [26, 48], and has also been developed in the pinch technique formal- 
ism [43, 24]. Most DS calculations have been done with the Faddeev-Popov action and, because they 
are non-perturbative, the question arises whether Gribov copies are treated correctly. To address this 
question, we must review the derivation of the DS equations. The functional DS equation, which 
expresses the tower of the DS equations for all correlators, results from the identity, 

f 5 

= j^dA^det[M{A)]exp[-SYM{A)/h+iJ,A)/fi], (4.1) 



95 



which states that the integral of a derivative vanishes. Here M(A) = —d^D^{A) is the Faddeev-Popov 
operator and, for Landau gauge, the integral extends over transverse configurations 3 • A = 0, and we 
have introduced h to keep track of the loop order in a diagram. This identity is correct, provided 
that the integrand vanishes on the boundary dR of the region of integration R in (transverse) A-space. 
The region of integration R is usually assumed to be all of A-space. However the identity also holds 
if R is the Gribov region, R = Q., because the Faddeev-Popov determinant, det[M(A)] vanishes on its 
boundary dQ., which is the first Gribov horizon. Indeed, because the Faddeev-Popov determinant is the 
infinite product of eigenvalues det[M(A)] = On ^nl^), the n-th Gribov horizon, defined by ?i„(A) = 0, 
is a nodal surface of the integrand for every integer n, and the identity holds if the boundary dR 
coincides with any one of these nodal surfaces. From the last identity the functional DS equation for 
the generating functional of correlators, 

Zr{J) = [dA det[M(A)] exp[-5YM(A)/n + (/,A)/n] , (4.2) 

JR 

follows, 



Zr{J)=0, (4.3) 



where 5'eff(A) = Sfp{A) - hlndetM{A)} With Zr{J) = exp[WR(7) /h], where Wr{J) is the generating 
functional of connected correlators, the last equation reads 

'^^'^ = 6a^)[''^^^¥j)- ^'-'^ 

Next we change variable from Wr{J) to the quantum effective action Tr{A) by the Legendre transfor- 
mation, 

A(7) = ^; r«(A) = (/,A)-W«(7), /(A) = ^. (4.5) 
This yields the functional DS equation for the quantum effective action rs(A), 

^x;A + n®Ay (4.6) 



8A^(x) 5A^(x) V 5A 

where h'D^{x,y; A) is the gluon propagator in the presence of the source A, and iD(J?) is expressed 
in terms of rR(A) by 



s2p 



The DS equations for the one-particle irreducible correlators are obtained by expanding (A) in pow- 
ers of A, and equating like coefficients. 



Observe now that the functional DS equation depends only on ^eff and is independent of which region 
R, bounded by a nodal surface dR, is chosen, or if all of space is chosen in (4.1). Consequently, for 
each region R that is bounded by a nodal surface — and there are an infinite number of nodal surfaces 
— there is a different solution Zr{J) of the DS equation for the generating functional Z{J), and more- 
over, because the equation for Z(/) is linear, a linear combination of solutions Zr is also a solution for 
Z{J). Each different solution Zr{J) gives a different solution Tr{A), because rR(A) is obtained from 

^In practice Faddeev-Popov ghosts are introduced, but the argument still stands. 
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Zr{J) by an invertible change of variable. Thus there is an enormous ambiguity in the solutions to 
the DS equation in Faddeev-Popov theory [46]. One of these solutions belongs to the choice R = Q., 
which corresponds to a cut-off at the Gribov horizon dQ.. The good news is that the DS equations for 
Faddeev-Popov theory do have a solution that corresponds to a cut-off at the Gribov horizon.^ Taking 
the optimistic view, if we impose enough additional conditions, which result for example from the 
positivity of the integrand detMexp(— 5) in the Gribov region Q., we may resolve the ambiguity and 
select this particular solution. The bad news is that, in practice, we know little about what additional 
conditions would be necessary to resolve this ambiguity completely. In view of this, it is perhaps not 
surprising that two different solutions have been found to the DS equations in Faddeev-Popov theory, 
namely "scaling" and "decoupling" solutions [26]. 

The DS equations for the propagators have also been derived from the GZ action, and its infrared 
behavior was investigated, with two solutions being found [100]. Subsequently one of these solutions 
could be eliminated [101]. Interestingly, the surviving solution has precisely the same infrared critical 
exponents as in Faddeev-Popov theory in the Landau gauge [46, 47]. 

4.5 Stochastic quantization 

An avenue for non-perturbative calculations which is insufficiently explored is stochastic quantization 
[102] with stochastic gauge fixing [103]. This method elegantly by-passes the problem of Gribov 
copies, and is geometrically unobjectionable, whereas the GZ action relies on an approximate solu- 
tion of the Gribov problem in which the fundamental modular region A is replaced by the Gribov 
region Q., as discussed in sect. 2.2. We discuss briefly the advantages and challenges presented by 
stochastic quantization, beginning with a brief introduction to this approach in case it is not familiar 
to the reader. 

4.5.1 Background 

Stochastic quantization is conveniently defined by a quantum field theoretic version of the time- 
dependent Fokker-Planck equation 



This equation determines the probability distribution P{x,t) of Brownian particles subject to a drift 
force Ki{x). It results from the local flow equation P = —V • j, where the current, j = —hVP + KP, 
consists of a diffusion or fluctuation term, where h is the diffusion constant, and a drift term. If the 
drift force is conservative, K = — W(x), the probability distribution P{x,t) relaxes to the Boltzmann 
distribution limt^ocP{x,t) = exp[— V (x) / h], for any initial distribution, provided that exp[— V (x)//i] 
is normalizable. 



^Alternatively, it has been argued [98] that integrating over all of A-space gives the correct answer by systematic cancel- 
lation of Gribov copies because the signed Faddeev-Popov determinent gives the signed intersection number of the gauge 
orbit with the gauge-fixing surface, which is a topological invariant. The question remains: how to select this particular 
solution? 




(4.8) 
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Consider now the functional version of this equation, 



where P{A,t) is a time-dependent probability distribution in Euchdean A-space, t is an artificial "fifth 
time" that corresponds to machine time in Monte Carlo simulations of lattice gauge theory, K^{x; A) 
is a "drift force" in A-space, and X/j is a Euclidean position-vector. If the drift force is conservative, 
K^{x; A) = where the Euclidean action 5(A) is the analog of the potential energy V(x), then 

PiA,t) relaxes to the Euclidean probability distribution, Peq{A) = \imt^coP{A,t) = A^exp[— 5'(A)/^], 
for any initial probability distribution, provided exp[—S{A)/h] is normalizable. Euclidean quantum 
field theory is recovered by calculating expectation-values from the equilibrium distribution 

{0) = j dAO{A)P,^{A). (4.10) 

In a gauge theory, the (would-be) equihbrium distribution exp[— 5'ym('4)//i], where Sym{A) is the 
Euchdean Yang-Mills action, is not in fact normahzable, and it appears that we are back with the 
original problem."* However in gauge theory, we are interested only in gauge-invariant observables. 
For these one may modify the drift force by introducing a gauge-fixing "force" that has the form of 
an infinitesimal gauge transformation, 

because this has no effect on the expectation-value of gauge-invariant observables [103]. The only 
scalar with the correct engineering dimension and color dependence is given by CO^ = a~^d-A^, and 
we take for the drift force, 

<(x;A) = - ^ + a-iDf ^.A^ (4.13) 

where a is a gauge parameter. The second term is a globally restoring force for a > 0, in the sense that 
it always points back toward the origin in A-space, (A, A'gf) < 0, for we have 

a^^{A^,D^d-A) = -a^^{d-A,d-A) < 0. (4.14) 

Expectation- values are calculated from the equilibrium distribution Peq{A) = ^mf^^P{A,t) deter- 
mined by the Fokker-Planck equation with drift force (4.13). 

This method of gauge fixing, that rehes on a gauge-fixing force that is an infinitesimal gauge trans- 
formation, is a geometrically correct procedure that bypasses the problem of Gribov copies. It is not 
available in an action formalism because the gauge-fixing force is non-conservative, and cannot be 
written as the gradient of some gauge-fixing action, 

(4.15) 



^In the original version of stochastic quantization [102], this was dealt with by taking the limit of the expectation- values 
of gauge-invariant observables, 

(O) = lim f dA 0{A)P{A,t) (4.11) 



98 



Note that this is a property of non-Abelian gauge groups. For the U{1) gauge group, we have 
Dfj{A) = d/j, and the gauge-fixing force is conservative, K^f jj = a^^d/jd - A = where Sgf = 

{2a)-^fd'^x{d-A)\ 

4.5.2 Time-dependent approach 

Calculations may be done using either the time-dependent or time-independent Fokker-Planck equa- 
tion. In the time-dependent approach, the solution P{A,t) is expressed as a path integral over paths 
A^{x,t), analogous to an integral over Brownian motion paths with a drift force. Because of the re- 
semblance of the Fokker-Planck to the Schrodinger equation, this path integral formally resembles 
the Feynman path integral in quantum mechanics, but with a local action that is 5-dimensional (and 
should not be confused with the 4-dimensional Euclidean action). It has been shown by BRST meth- 
ods that this method is renormalizable [104, 105, 106, 107]. However calculations are challenging in 
this approach. For example, in a perturbative expansion, the zeroth-order gluon propagator is given by 
[kl + {k^)^]~^, where = 1., , ,4. To our knowledge the time-dependent Fokker-Planck has not been 
used in non-perturbative calculations. 



4.5.3 Time-independent approach 

The equilibrium probability distribution P^q = Mmt^„P{A,t), which is used to calculate expectation 
values of gauge-invariant observables, satisfies the time-independent Fokker-Planck equation, 



A) 



Peq(A), 



(4.16) 



Because the force is not conservative, we cannot give a closed-form solution for the equilibrium 
distribution Peq{A). The basis of the time-independent approach is that this equation may be con- 
verted by a change of variables to a useful functional equation for the correlators, that serves as a 
substitute to the functional DS equation in a theory based on a local action [109]. Before doing so, 
we note that the solution Peq(A) to this equation (if it exists) has no nodal surfaces, and is unique and 
positive up to normalization. This is shown in Appendix C. In contrast, the Faddeev-Popov weight 
det[M(A)] exp[— 5ym(A)/^] has many nodal surfaces, and this leads to a great ambiguity in the solu- 
tion of the functional DS equation in Faddeev-Popov theory, as we have seen. 

To express the time-independent Fokker-Planck equation as an equation for the correlators, we multi- 
ply by exp[(7, A) / h] and integrate over all A. After integration by parts we obtain 



where 



= j d'^x J^{x) J^{x)+K^ i^x; h— 
Z{J) = j dA exp[(y, A)/h] Peq(A) 



Z(7), 



(4.17) 



(4.18) 
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is the generating functional of the correlators. [If we replace J by //, then Z{iJ) = Peq{J) is the 
functional fourier transform of Psq{A).] With the change of variable Z{J) = exp\W{J)/h], where 
W{J) is the generating functional of connected correlators, this reads 



= Jd'^x 



m ^5 



(4.19) 



Finally we make the change of variable from W{J) to the quantum effective action T{A) by the Leg- 
endre transformation (4.5), which expresses the time-independent Fokker-Planck equation as a func- 
tional equation for the quantum effective action, 







5r 



(4.20) 



where h1f^{x^y\ A) is the gluon propagator in the presence of the source A, and is expressed in terms 
of r(A) in (4.7), and = + %J. The term [x; A + h(D^ occurs in the familiar 

functional DS equation in Faddeev-Popov theory (4.6), and we shall not write it out here, while the 
second term, which results from stochastic gauge fixing, is simpler and is given by 

K^il[x-A + h'D^ =a-'Dl^{A)d-A\x)+a-'hgf'''d^^-Dli{xj;A)\y=,. (4.21) 



Although the time-independent Fokker-Planck equation in the form (4.20) appears to have a quite 
different structure from the DS equation (4.6) in Faddeev-Popov theory, it may be used, like the DS 
equation in Faddeev-Popov theory, to generate a perturbation series for the correlators or, more am- 
bitiously, to obtain a non-perturbative solution [109]. In a perturbative expansion, the parameter h 
counts the number of independent loops. This method has the advantage that the solution to equation 
(4.20) is unique. This is true because the solution to the time-independent Fokker-Planck equation 
(4.16), is unique (up to normalization) as shown in Appendix C, and because the change of variable 
from (4.16) to (4.20) is invertible.^ In contrast, the DS equation in Faddeev-Popov theory has many 
different solutions because of the existence of an infinite number of nodal surfaces of the Faddeev- 
Popov determinant detM(A), as explained above. 



In [46] the Landau-gauge limit, a — 0, of this equation was written as the DS equation in Faddeev- 
Popov theory plus a correction term. With neglect of the correction term, that is, from the DS equa- 
tion in Faddeev-Popov theory, the infrared critical exponent of the gluon propagator was found to 
be D{k) ~ {k^f/^, {k^f-^'^^^ and (it2)0 i906 -^^ Euchdean dimension d = 2,3, and 4 respectively [46]. 
The result for li = 4 agrees with the independent calculation of [47]. An improved calculation, based 
directly on the time-independent Fokker-Planck equation (4.20), yielded instead D{k) ~ (^^)" "'*^ in 
Euclidean dimension d = 4 [109], which vanishes at A: = 0, but very weakly. The smallness of the 
exponent is notable because the vanishing of D{k) atk = was imposed by hand (rightly or wrongly) 
as a boundary condition. The result, D{k) (A^)" *''^'' is in reasonable agreement with the (subse- 
quent) lattice calculations on very large lattices, [29] and other references cited above, which favor 
D{k) (^•^)° foT d = 3 and 4, taking into account that there is an unestimated truncation error in the 
analytic calculation from (4.20), and numerical uncertainty in the lattice calculation. 

^The uniqueness of the solution of the exact equation (4.20) does not necessarily imply the uniqueness of the solution 
by a truncation of (4.20). 



100 



In summary, use of the time-independent Fokker-Planck equation (4.20) has two great advantages 
of principle. First the Gribov problem is by-passed by gauge fixing with the gauge-fixing force that is 
an infinitesimal gauge transformation. Second, the solution to (4.20) is unique. The price paid is that 
the time-independent stochastic dynamics involves a non-conservative (but local) drift force, and this 
exceeds the bounds of conventional quantum field theory that is formulated in terms of a local action. 
So if we wish to pursue this approach, we must abandon the comfortable terra firma of conventional 
quantum field theory, and embark on an unexplored sea. New techniques would have to be developed 
to establish renormaUzabiUty, Ward identities, and other needed properties. 

Assuming these results can be established, let us speculate on the characteristics of the time-independent 
stochastic dynamics defined by eq. (4.20). First we should expect that the singularities of correlators 
of gauge-non-invariant fields will be unphysical, as happens with the GZ action. This is practically as- 
sured because (1) the probability distribution gets concentrated inside the Gribov region in the Landau 
gauge limit, a — ^ 0, of (4.20) [46], as it does with the GZ action, and (2) violation of reflection positiv- 
ity in the gluon propagator has been confirmed by numerical simulations in lattice gauge theory. There 
is also a simple a priori reason to expect unphysical singularities to result from (4.20): because of the 
non-conservative drift force, the rules of calculation deviate from the conventional Feynman rules, 
so the singularities that they determine should also be different. If the singularities are unphysical, 
as seems to be assured, then the corresponding particles are confined. If so, one needs a mechanism 
which assures that the singularities of the correlators of gauge-invariant objects, such as glue-balls or 
hadrons, should remain physical. For this too there is a simple argument. The correlators of gauge- 
invariant objects are independent of the gauge parameter oT^, and thus exist in the limit ^ 0. But 
is the coefficient of the gauge-fixing drift force, so in the limit — > 0, the non-conservative 
drift force weakens, and the drift force approaches the physical Urriit, limQ-i_^o^^ = Thus, 
for the correlators of gauge-invariant objects, we expect to recover the physical singularity structure 
associated with a local action. Note also that the gauge-fixing force K^f provides a restoring force that 
is tangent to gauge orbits and, as this force weakens in the limit ^ 0, the probability escapes to in- 
finity along gauge-orbit directions, and we may expect the correlators of gauge-non-invariant objects 
to vanish in this limit, as happens in (non-gauge-fixed) lattice gauge theory. In this scenario in the ver- 
sion of continuum gauge theory provided by time-independent stochastic quantization, confinement 
of gauge non-invariant fields is an almost kinematic consequence of non-Abelian gauge invariance. 
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Appendix A 

Some formulae and extra calculations 



A.l Gaussian integrals 
A.l.l Gaussian integral for scalar variables 

I{AJ) = j [d(p]exp d!^xd!^y(Sf{x)A{x,y)(Sf{y) + j x (Sf{x)J {x 

= C{detA)-^/^exp^Jd''xd'^yJ{x)A-\x,y)J(y), (A.l) 
with C an infinite constant, which, in practice, can always be omitted. 

A.1.2 Gaussian integral for complex conjugated scalar variables 

I{A,J) = J [dcp][d(p]exp - J d^xd^'j cp(x)A(x,j)cp();) + J d'^x{^ix)J^{x) +(!?{x)J^{x)) 

= C{detA)-^exp Jd''xd'^yJ^{x)A-\x,y)J^{y) , (A.l) 
again with C an infinite constant. 

A.1.3 Gaussian integral for Grassmann variables 

/(A,'n,fi) = y"[de][de]exp Jd''xd''ye{x)A{x,y)e{y) + Jd''x{f\{x)e{x) + e{x)r\{x)) 

= CdetAexp- J d''xd''y\\{x)A-\x,y)ii\{y), (A.3) 
again with C an infinite constant. 

A.2 D^(A)co = is a gauge invariant equation. 

To prove that D^(A)to = is a gauge invariant equation, we can write from expression (1.18), 

Djj(o = 3^00 - igA^ro + igcoA^ • (A.4) 
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Now performing a SU{N) transformation, we know that D^a is in the adjoint representation by defi- 
nition, 

D/jO) = ^ UD^iSiU^ = , (A.5) 

so working out this equation we find, 

UD^(SiU'' = Udf.au'' -igUAf^U^UaU^ +igUGiU'^UAf,u'' 

= a^co' - {d^U)(0U^ - U(0{d^,U^) - ig{A'^ + -d^UU^)(o' + ig(o'{A'^ + -d^,UU^) 

= (A.6) 
whereby we made use of equation (1.15) and the simple formula 

f/f/t = i ^ d^UU^ + Ud^U^ = 0. (A.7) 
We have thus indeed proven that D^(A)co = is a gauge invariant equation. 

A.3 o decreases with increasing k^. 

We shall prove that the following function 

decreases with increasing k^. Let us prove this in 2 dimensions for simplicity. We assume k = {kx,ky) 
to be oriented along the x axis. Using polar coordinates, we obtain 

4^ r ^/ 2n /■^''.q l-cos^e 

= Jo (2^-'^^'^ Jo "' k^W-hcose 

whereby we have used the following Poisson-like 6-integral which can be easily calculated using a 
contour integration, 

f^^ l-cos^e _ / ^ if k^<q^ 

Jo '^^k^ + q^-2qkcose ~ \ % ifq^<k^ - ^^"^^ 
Now deriving f{k,A) w.r.t. k^ and using the property of the 9 function: ^Q{x — y) = 5{x — y), we find 

and thus f{k,A) is a decreasing function for increasing A:^. 
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A.4 Determinant of K, 



We calculate the determinant of 



1 2 Ng^ 



VdN^-l ""k^ 




a 



(A. 11) 



-1/2 



We can write 

(det<w) 
Therefore, we need to determine 

TrlnK^ = (A^^ _ i)Trln ^5^k + (^^ 

= (N^-l) Trln(^5^v(^+^'))+Trln(^5^ + 



(A.12) 



1 ) ^K^V 



1 I ^/J^V 



(A^'-l) 



(A.13) 



whereby we used ln(l +x) = x — y + . . .. We can now take the trace of the diagonal elements of the 

2 

second term, and again use x — + . . . = \n{l + x). We obtain. 



Trln*:^^ 



(A^'-l) 



k k 
in 

k 



1 



1 



J^ln 



1 ^ 

it2+ a 



(A.14) 



By working out the last term, we see that it is proportional to a. 



k 



r d'^k f k^ n-\ r d'^k k^ n- 
V / TTT^ln ^+iV^ +y / T^^ln ^-iV^ 



^ a 



J/4 



(A. 15) 



whereby / d^^ln^72 is zero in dimensional regularization. Therefore, in the Umit a 0, becomes 
zero. In conclusion, we find 



(det<(fc)) 



-1/2 



exp 
exp 



(d-l) 



id -I) 



2 J (Iny k^ 
r d'^k 

J (271)^ 



hi U + 



2 , 1 2 (3jVg2 1 



(A.16) 
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Appendix B 

Cohomologies and the doublet theorem 



B.l Cohomology 

Suppose 8 is a nilpotent operator, 8^ = 0. The cohomology of 8 is given by the solutions of the 
equation 

8A = 0, (B.l) 

which cannot be written in the form 

A = Sa . (B.2) 

A quantity A obeying equation (B. 1) is called closed, while a quantity of the form (B.2) is called exact. 
The cohomology of A is thus identified by quantities which are closed but not exact. More precisely, 
a non trivial quantity A is always defined up to the addition of an arbitrary exact part, i.e. one speaks 
of cohomology classes. In fact, take now two closed quantities Ai and A2. These quantities belong to 
the same cohomology class if 

A1-A2 = 8(...), (B.3) 
i.e. when Ai and A2 differ by an exact part. 

In this way one can always write A obeying (B.l) as a sum of a trivial part and a non trivial part. 

A = A„.ttiv + 80^, (B.4) 

whereby An.triv does not contain parts that can be written as 8(. . .). In quantum field theory, these non 
trivial parts shall be the most interesting parts, as they will be related to the renormalization of the 
physical parameters of the theory. 

B.2 Doublet theorem 

Now there is a very important theorem which shall be very useful later on [1 16]. Suppose our theory 
contains a pair of fields, sources or parameters (m,-, v,) which form a doublet: 

8m,- = Vi 8vi = , (B.5) 
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whereby the subscript / is a certain index (e.g. color). We assume Uj to be commuting, while v, is an 
anticommuting quantity. Then we can prove that m,- and v, shall never enter the non trivial part of the 
cohomology of 8. 

The proof is as follows. We introduce two operators P and A 

P = dx[ui^+Vi^ 

J \ OUi ov, 



A = Jdx (^Ui^^ , (B.6) 
Functionally, we write for the nilpotent operator 6 

6 = v,-^, (B.7) 



OUi 



so we obtain 



A = / dx Ui^ UiVjz—zr — 

J \ OUi OVi OUj ) 



and thus 



Analogously we also have 



{6,A} = A (B.8) 



Ph = / dx M/Vj^— 5 \-Vi^ 

J \ OUi OUj OUi OVi OUj 

„ . f f d 3 3 3 

8 P = dx[vj^+ VjUi^^ + VjVi^ 

J \ OUj OUj OUi OUj 



3 d_ 

3v; 



and thus 

[8,P]=0. (B.9) 

As P is a counting operator for the total number of m,- and v,-, we can expand^ A, see expression (B.4), 
in eigenvectors of P, 

A=£A„, (B.IO) 

n>0 

whereby PA„ = nA„ and n represents the total number of m,- and v, in A„. Now from the cohomology 
condition (B.l) and the commutation relation (B.9), we find that 

0=£n8A„ ^n8A„ = 0. (B.ll) 

n>0 n>l 

Looking at expression (B.5), we easily obtain that 

8A„ = Vn>l. (B.12) 



^We assume A to be a polynomial in and Vj. 
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Finally, using this property and invoking expression (B.8), we obtain 

A = Ao+Y-PAn 

«>i" 

= Ao+V-5AA„ 

= Ao + 5(...). (B.13) 

In conclusion, as 6^ = 

5A = 5Ao, (B.14) 

whereby Aq is independent of the doublet (m,-, v,). The quantities m,- and v, shall thus never enter the 
non trivial part of the cohomology. 
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Appendix C 

Positivity and uniqueness of solution of 
time-independent Fokker-Planck 
equation 



C.l Positivity 

The proof is given in a finite number of dimensions, but that should be a good guide. We use discrete 
notation and write A,- instead of A^{x), so the time-independent Fokker-Planck equation reads, 

= di[-di + Ki{A)]F{A), (C.l) 

where 3, = J^. Suppose that this equation possesses a solution F{A) that is positive somewhere and 

negative elsewhere. Call R the region where F(A) is negative, so F{A) < for A in the interior of 
R, it is positive outside, F{A) > for A ^ /?, and the boundary dR of the region 7? is a nodal surface, 
F{A) = for A G dR. Now integrate the equation over the region R, 

0= fdA di[-di + Ki{A)]F(A), (C.2) 
Jr 

where dA is the volume element in A-space, and use Gauss's law to obtain 

= [ dLi[-di + Ki{A)]F{A), (C.3) 

JdR 

where dUi is the surface element oriented along the outward normal n,- on dR. By hypothesis F{A) 
vanishes on the boundary dR, which gives 

= / d-LidiFiA). (C.4) 

JdR 

Moreover F{A) is negative inside R and positive outside, so F{A) increases as one leaves the region 
R, and the gradient along JZ, is positive, dLi 3,F (A) > 0, which is a contradiction. Consequently the 
solution F{A) is either positive everywhere or negative everywhere, and F{A) caimot have any nodal 
surface. We may normalize F{A) so it is positive everywhere F{A) > 0, with / dA F{A) = 1. 
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C.2 Uniqueness 



Now suppose that there are two Hnearly independent solutions to the time-independent Fokker-PIanck 
equation, Fi{A) and F2{A), and that each is normalized / dA Fi{A) = f dA F2{A) = I. Because the 
time-independent Fokker-PIanck equation is linear, the difference F{A) = Fi{A) — F2(A) is also a 
solution, and it satisfies / dA F{A) = 0. However we have just proven that F{A) is either positive 
everywhere or negative everywhere, or zero, and we conclude F{A) = Fi{A) — F2{A) = 0. Thus if 
the time-independent Fokker-PIanck equation possesses a non-zero solution Fi{A), this solution is 
positive and unique up to normalization. This result in no way depends upon whether the drift force 
Ki{A) is conservative. However when Ki{A) is non-conservative, the equilibrium current does not 
vanish jgq ; = — ^i^eq + ^!^eq / 0^ for if jgq i = then the drift force is given by Ki = dilnP^^, and is 
conservative. 
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